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The Belle experiment, running at the KEKB e + e asymmetric energy collider during the 
first decade of the century, achieved its original objective of measuring precisely differences 
between particles and anti-particles in the B system. After collecting 1000 fb _1 of data at 
various T resonances, Belle also obtained the many other physics results described in this 
article. 
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§1. Introduction 

In the sections that follow, we describe the physics accomplishments of the Belle 
experiment, which ran at the KEKB 1 ) e + e _ asymmetric energy collider in Tsukuba, 
Japan between 1999 and 2010. KEKB broke all records for integrated and instanta- 
neous luminosity for a high energy accelerator. As a result Belle was able to integrate 
over 1000 fb _1 or one inverse attobarn of data. 

Belle was designed and optimized for the observation oi CP violation in the B 
meson system. In 2001, Belle (along with BaBar, a competing and similar experiment 
located in Stanford, California) was indeed able to observe large CP asymmetries 
in B decays, which were expected and consistent with the theoretical proposal of 
Kobayashi and Maskawa. This experimental result was explicitly recognized in the 
2008 Physics Nobel Prize citation. 

Nevertheless, the Belle spectrometer was a general purpose device with reason- 
able solid coverage as well as high quality vertexing with silicon strip detectors, 
charged particle tracking with a central drift chamber, and excellent electromagnetic 
calorimetry as well as muon and Kl detection. These detector capabilities allowed 
Belle to not only cover most of the important topics in B physics (in addition to 
the CP violation measurements) but also to make important discoveries in charm 
physics, tau lepton physics, hadron spectroscopy, and two-photon physics. 

Most of the Belle luminosity was recorded on or near the T(4S) resonance, which 
is the optimal center of mass (CM) energy for the production of BB pairs used in 
B physics analysis. However, KEKB has some flexibility in energy and Belle also 
recorded a series of unique data sets at the T(IS), T(2S), and T(55) resonances. 
The latter data set is of special interest in hadron spectroscopy as a large number of 
new and some exotic states were found in analyses of this sample. 

§2. The Belle detector and its data samples 

2.1. Overview 

The Belle detector is located at the interaction region of an asymmetric energy 
e + e~ collider, called KEKB. 1 ) Belle is optimized to measure time-dependent CP 
violation in B-meson decay. Therefore, the detector has good vertex resolution 
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and good particle identification capabilities for leptons and hadrons. The detector 
material is minimized to reduce multiple scattering for charged particles and to 
maintain high efficiency and good resolution for low energy photons. The acceptance 
is asymmetric (covering the polar angle region from 17° to 150°) to match the boost 
from the asymmetric 8 on 3.5 GeV energy collisions. Belle is a general purpose 4tt 
detector, which can accommodate various physics programs, including studies of r 
pairs, two-photon physics, and qq continuum processes. 

Figure 1 shows the Belle detector configuration. The detector is built around a 
1.5 Tesla superconducting solenoid and iron structure. The beam crossing angle is 
±11 mrad. i?-meson decay vertices are measured by a double-sided silicon vertex 
detector (SVD) situated around a cylindrical beryllium beam pipe. There are two 
inner detector configurations: SVD1 (three layers before the summer of 2003) and 
SVD2 (four layers). Charged particle tracking is provided by a central drift chamber 
(CDC). Particle identification is provided by dE/dx measurements in the CDC, aero- 
gel Cerenkov counters (ACC), and time-of-flight counters (TOF) situated radially 
outside of the CDC. Electromagnetic showers are detected by an array of CsI(Tl) 
crystals (ECL) located inside the solenoid coil. Muons and Kl mesons are identified 
by arrays of resistive plate counters (KLM) interspersed in the iron yoke. An array 
of bismuth germanate oxide (BGO) crystals called the extreme forward calorimeter 
(EFC) is located on the surface of the cryostats of the compensation solenoid mag- 
nets in the forward and backward directions. The EFC is used as an active shield 
against beam background and also measures the online luminosity. Each subdetec- 
tor is briefly described in the following subsections; more detailed information is 
available in Ref. 2 ) 

2.2. The beam pipe 

The Belle detector beam pipe 3 ) is connected to the KEKB accelerator beam 
pipe. The pipe is a double-wall beryllium structure with liquid paraffin cooling to 
remove the heat generated by the beams. The inner diameter is only 30 mm (40 mm) 
for the SVD2 (SVD1) inner detector configuration to optimize vertex resolution. A 
10 //m-thick layer of gold is sputtered inside the beryllium wall to prevent synchrotron 
radiation photons from entering the detector. 

2.3. The SVD detector 

The SVD consists of four layers in a barrel-only design (SVD2 4 )). Each layer is 
independently constructed and consists of ladders. Each ladder consists of double- 
sided silicon strip detectors (DSSDs) reinforced by support ribs. The design uses 
two types of DSSDs. For the second version, the DSSDs of the 4th layer are shorter 
and wider than those of the other layers. The readout chain for the DSSDs is 
based on a VA1 (Viking architecture) integrated circuit. The back-end electronics 
is a system of flash analog-to-digital converters (FADCs) and field programmable 
gate arrays (FPGAs), which perform online common- mode noise subtraction, data 
sparsification, and data formatting. 

Before the summer of 2003, there were three DSSD layers with slightly less 
angular coverage (SVD1). In addition, the beam pipe diameter was larger (40 mm 




Fig. 1. Side view of the Belle detector. 

versus 30 mm for the final data taking configuration). 

2.4. The CDC detector 

The structure of the CDC is asymmetric in the z direction in order to optimize 
the angular coverage. The longest wires are 2400 mm long. The inner radius extends 
to 80 mm without any walls in order to obtain good tracking efficiency for low 
momentum tracks with minimal intervening material. The outer radius is 880 mm. 
The forward and backward smaller-r regions have conical shapes in order to clear 
the accelerator components while maximizing the acceptance. A low-Z gas, a He- 
-C<iH% (50/50) gas mixture, is used in order to minimize multiple scattering. The 
chamber has 50 cylindrical layers, each containing between three and six axial or 
small-angle stereo layers, and three cathode strip layers. The CDC has a total of 
8400 drift cells. We chose three layers for the two innermost stereo super-layers and 
four layers for the three outer stereo super-layers in order to provide a highly efficient 
and fast z-trigger, which is combined with the information from the cathode strips. 

During the summer of 2003, the cathode part of the CDC was replaced by a 
compact small cell type drift chamber in order to make enough space for the SVD2 
vertex detector. The cell sizes are only 5 mm in both the radial and azimuthal 
directions to accommodate two layers (128 cells per layer) in a limited space. The 
maximum drift time is rather small (~100 nsec); this feature can provide the first 
trigger signal for the SVD2 readout latch. 
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2.5. The ACC subsystem 

The ACC consists of 960 counter modules segmented into 60 cells in the (ft 
direction for the barrel part and 228 modules arranged in 5 concentric layers for 
the forward end-cap part of the detector. All the counters are arranged in a semi- 
projective geometry, pointing to the interaction point (IP). In order to obtain good 
pion/kaon separation to cover the entire kinematical range of two-body B decays, the 
refractive indices of the aerogel blocks vary between 1.01 and 1.03, depending on their 
polar angle region. Five aerogel tiles are stacked in a thin (0.2 mm thick) aluminum 
box of approximate dimensions 12 x 12 x xl2 cm 3 . In order to detect Cerenkov light 
effectively, one or two fine mesh-type photomultiplier tubes (FM-PMTs), which are 
operated in a 1.5 T magnetic field, are attached directly to the aerogel on the sides 
of the box. We use Hamamatsu Photonics PMTs of three different diameters: 3, 

2.5. and 2 inches, depending on the refractive index of the aerogel block, in order to 
obtain a uniform response for relativistic particles. 

2.6. The TOF subsystem 

The TOF system consists of 128 TOF counters and 64 thin trigger scintillation 
counters (TSC). Two trapezoidal shaped TOF counters and one TSC counter, with a 
1.5 cm intervening radial gap, form a single module. In total, 64 TOF/TSC modules 
located at a radius of 1.2 m from the IP cover a polar angle range from 34° to 120°. 
The thicknesses of the scintillators (BC408, Bicron) are 4 cm and 0.5 cm for the 
TOF and TSC counters, respectively. The fine mesh PMTs operating inside the 1.5 
T magnetic field, with a 2-inch diameter and 24 stages, were attached to both ends 
of the TOF counter with an air gap of 0.1 mm. For the TSCs, the tubes were glued 
to the light guides at the backward ends of the counters. 

2.7. The ECL detector subsystem 

A highly segmented array of CsI(Tl) crystals with silicon photodiode readout 
were selected for the ECL. 5 ) Each crystal has a tower-like shape and is arranged so 
that it nearly points to the IP. The calorimeter covers the full Belle angular region. 
A small gap between the barrel and end-cap crystals provides a pathway for the 
cables and room for supporting members of the inner detectors. The entire system 
contains 8736 counters. The size of each crystal is typically 55 mm x 55 mm (front 
face) and 65 mm x 65 mm (rear face). The 30 cm length (16 radiation lengths) is 
chosen to avoid deterioration of the energy resolution for high energy gammas due 
to fluctuations in the shower leakage out the rear of the counter. Each counter is 
read out by an independent pair of silicon PIN photodiodes and charge sensitive 
preamplifiers attached at the end of the crystal. 

2.8. The KLM detector 

The KLM consists of alternating layers of charged particle detectors and 4.7 cm- 
thick iron plates, which are the magnetic flux return in the barrel and endcap re- 
gions. 6 ) There are 15 detector layers and 14 iron layers in the octagonal barrel region 
and 14 detector layers in each of the forward and backward end-caps. The iron plates 
provide a total of 3.9 interaction lengths of material for a particle traveling normal to 



Physics achievements from the Belle experiment 



7 



the detector planes. The detection of charged particles is provided by glass-electrode 
resistive plate counters (RPCs). The resistive plate counters have two parallel plate 
electrodes of 2.4 mm-thick commercially available float glass. The bulk resistivity of 
the glass is 10 12 — 10 13 ,!? cm at room temperature. To distribute the high voltage 
on the glass, the outer surface was coated with carbon ink, which achieves a surface 
resistivity of 10 6 — 10 7 ,!?/square. The discharge signal can then be obtained from 
external pickup strips. The readout of 38K pickup strips is accomplished with the 
use of custom-made VME-based discriminator /time multiplexing boards. 

2.9. Trigger and data acquisition 

The Belle trigger system consists of a Level- 1 hardware trigger and a Level-3 
software trigger. 7 ) The latter is implemented in the online computer farm. The 
Level-1 trigger system consists of a subdetector trigger system and a central trigger 
system called the global decision logic (GDL). The subdetector trigger systems are 
based on two categories: track triggers and energy triggers. The CDC and TOF 
are used to yield trigger signals for charged particles. The ECL trigger system 8 ) 
provides triggers based on total energy deposit and cluster counting of crystal hits. 
These two categories have sufficient redundancy. The KLM trigger gives additional 
information on muons. The EFC triggers are used for tagging two-photon events as 
well as Bhabha events. 

The Belle data acquisition system used one type of multi-hit TDC modules for 
all subsystems except for the SVD. The signal pulse height is recorded as timing 
information using a charge to time conversion chip (Q-to-T chip). Precise timing 
information in the TOF is recorded by commercial TDC modules with special time 
expansion modules. The TDC modules did not have a pipe-line readout scheme. 
Therefore, the readout deadtime is large (around 30 ^sec). There were several elec- 
tronics upgrades in order to reduce deadtime carried out during the latter parts of 
the experiment. The TDC modules were gradually replaced with pipe- lined TDCs 
(2.8 /usee) for most of the subdetectors in the 2007-2009 running period. 9 ) It was 
carefully checked that these electronics upgrades did not affect the data quality. 

Belle turned off the detector high voltage during beam injection, as do other 
experiments. The KEKB injection time was slightly longer than at PEP-II (the 
collider hosting the BaBar experiment) and the average efficiency was lower. In 
order to reduce such losses, a continuous injection scheme 1 -* was implemented in 
January 2004. The detector high voltage was kept on and the trigger signals were 
vetoed for a 3.5 msec interval just after each beam injection. This scheme leads to 
3.5% deadtime only in the case of a 10 Hz injection rate. After adopting continuous 
injection, the KEKB machine became more stable and the peak luminosity improved 
due to the leveling of the beam currents. 

2.10. Detector performance 

The charged track reconstruction mainly uses the CDC. Good momentum res- 
olution is obtained by combining CDC tracks together with SVD hit information, 
especially for low momentum tracks, thanks to the limited amount of intervening 
material. The following expression gives the momentum resolution for a charged 
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track as a function of its transverse momentum: 

a pt / Pt = 0.0019 p t [GeV/c] 0.0030/^ 



(2-1) 



The typical mass resolution of D° mesons is 5 MeV in hadronic events. The z-vertex 
resolution is 61 /iin in the J/tp — >■ ^ mode. A similar resolution is also obtained 
in the r — <p plane. The energy resolution of the ECL is 1.7% for Bhabha events. A 
7T° mass resolution of 4.8 MeV is obtained for low momentum photons in hadronic 
events. 
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Fig. 2. Kaon identification efficiency and fake 
rate as a function of momentum. 
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Fig. 3. Muon identification efficiency and fake 
rate as a function of momentum. 



Pion/kaon/proton separation is obtained by combining ACC, TOF, and CDC 
dE/dx information. The kaon efficiency and the fake rate are shown in Fig. 2. The 
typical electron identification efficiency is 90% with a small fake rate (0.3%). Muons 
are also identified with 90% efficiency (2% fake rate) for charged tracks with momenta 
larger than 0.8 GeV (Fig. 3). More detailed information is available in Refs. 10 ^ 

2.11. Luminosity 

Belle started data taking on 1 June 1999. After that, data runs were taken 
for 6-9 months every year until the final shutdown on 30 June 2010. The total 
integrated luminosity reached 1040 fb _1 , as shown in Fig. 4. Belle took most of its 
data at the energy of the T(4S) resonance in order to study £?-meson decay. Off- 
resonance data were collected 60 MeV below the resonance peak energy for 10% of the 
running time about every two months in order to determine the non-BB background. 



Physics achievements from the Belle experiment 



9 



Belle logged luminosity 



The first non-T(4S') data were taken at the 
energy of the T(55) resonance for just three 
days in 2005. In the same year, T(3S) res- 
onance data were taken to search for invis- 
ible decay modes of the T(1S) resonance. 
The last T(4S) resonance data were taken 
in June 2008. During the last two years of 
operation, T(1S), T(2S), and T(5S) res- 
onance data samples were taken as well as 
energy scans between the T(4S) and T(6S) 
resonances. The integrated luminosity col- 
lected by Belle for each CM energy is listed 
in Table I and is calculated using barrel 
Bhabha events after removing bad runs, 
which could not be used in physics anal- 
ysis due to serious detector problems. The 
systematic error in the luminosity measurement is about 1.4%; the statistical error is 
usually small compared with the systematic error. Integrated luminosities for T(4S) 
data are shown separately for the SVD1 and SVD2 data sets, which were taken 
with different inner detector hardware configurations as described in the previous 
subsection. Other resonance and scan data were taken in the SVD2 configuration. 
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Fig. 4. Integrated 
Belle. 
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Table I. Summary of the luminosity integrated by Belle, broken down by CM energy. 



Resonance 


On-peak 


Off-peak 


Number of resonances 




luminosity (fb _1 ) 


luminosity (fb _1 ) 




T(15) 


5.7 


1.8 


102 xl0 b 


r(2S*) 


24.9 


1.7 


158 xlO 6 


Y{3S) 


2.9 


0.25 


11 xlO 6 


r(4S) SVDl 


140.0 


15.6 


152 xlO 6 BB 


r{4S) SVD2 


571.0 


73.8 


620 xlO 6 BB 


T(5S) 


121.4 


1.7 


7.1 xlO 6 B s Bs 


Scan 




27.6 





§3. CKM angle measurements 

3.1. The Kobayashi-Maskawa model and unitarity triangle 

The phenomenon of CP violation was one of the major unresolved issues in 
elementary particle physics after its discovery in 1964 in neutral kaon decay. 11 ) In 
1973, M. Kobayashi and T. Maskawa proposed a model in which a quark-mixing 
matrix among six quark flavors includes a single irreducible complex phase that 
causes CP violation. 12 ) Conventionally the quark-mixing matrix is written as: 13 ) 

V ud V us V ub \ ( l-A 2 /2 A AX 3 (p-ir])\ 

V cd V cs V cb \ = \ -A 1 - A 2 /2 ^A 2 + o(#)l) 

V td Vu V tb J \AX\\-p-iri) -AX 2 1 / 
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where the nontrivial complex phases are assigned to V ub and Vtd- Due to the unitarity 
of this matrix, the following relation is expected to hold, in particular for the terms 
involving the 6-quark: 



v td vi + v cd v; b + v ud v. 



'ub 



0. 



(3-2) 



This expression can be visualized as a closed triangle in the complex plane as shown in 



Fig. 

Here the phase of Vtd plays a fundamen- 
tal role and induces time-dependent CP 
asymmetries via interference with ampli- 
tudes containing V cb and V u b- Measure- 
ments of the relevant time-dependent CP 
violation parameters are used to determine 
the CP- violating angles, 4>i an d 02 , 14 ^ that 
are described in Sects. 3.4 and 3.5. In con- 
trast, the angle 03 is determined by the di- 
rect CP asymmetries in B — > DK^ decays 
and is discussed in Sect. 3.6. 



3.2. 



CP violation and B°-B° mixing 



Neutral B mesons, B° and -B u , can 
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v ud v ub 




(0,0) Vcd v cb * (1.0) 



Fig. 5. The unitarity triangle relevant to 
B decays. The CP violation param- 
eters are defined as the angles tf>\, <f>2 
and d>3. 



transform or mix into their antiparticles 

through box diagrams as shown in Fig. 6. The frequency of the mixing transition 
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Fig. 6. Box diagrams that contribute to B -B mixing. 



(oscillation) is Am d = (0.507±0.004) ps" 1 , 15 ) while the lifetime (t b o) is 1.519±0.007 



ps. 



15) 



A.I. Sanda, A.R. Carter, and I.I. Bigi showed that a sizable CP violation can 
appear in P-meson decays if B°-B° mixing is large. In neutral B decays to 
CP eigenstate (fcp), both B° and B can decay to the same final state. Because 
of B°-B° mixing, the decay proceeds through two paths; one from direct decay, 
B° — > fcp, and the other through B°~B° mixing, B° B° ->■ f C p. These two 
amplitudes have a phase difference of 4> m ix — %4>D where m i x is the weak phase of 
B°-B° mixing, axg{V t dV* b /V^ d Vtb) , and 4>d is the weak phase of the B° — > fcp decay. 
In the Wolfenstein representation, (p m i x = 2<pi and the phase difference is given as 
2(0i — 4>d)- The interference term for the two amplitudes has opposite signs for B° 
and B° decays and leads to CP violation effects proportional to sin2(0i — 4>d)- 
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3.3. Experimental approach at a B-factory 

At a P-factory, pairs of neutral B mesons in a coherent state with C = — 1 are 
produced by T(4S) — > B°B° decays In a decay in which one B meson decays to fcp 
and the other B meson decays to a flavor specific final state, / tag , the decay rate is 
given as 

e -\At\/r B0 r i 
V(At,q;S,A) = — <l + q- S Am d At) + A cos( Am d At) \. (3-3) 



At b o 

Here At = top — hag is the difference between the proper decay times of fcp and 
/tag) Q = ±1 is the flavor of /tag (+1 for B° — > /tag)- The quantities S and A are 
CP violation parameters that are dependent on the decay mode. The parameter S 
describes mixing-induced CP violation and is given by S = —ijcp sin 2((pi — 4>d), 
where rjcp is the CP eigenvalue of fcp- The other parameter, A, corresponds to 
direct CP violation (i.e. no CP violation in the B° <-> B° transition rates). It should 
be noted that, depending on the weak phase of the decay, CP violation measurements 
give information on the various angles of the unitarity triangle. The asymmetry in 
the rate of B° and B° decays is given by 

A{At)s l\^ +V j^-l\ A y l '- S ' A l = S - Am d At + A cos Am d At (3-4) 
v ' V(At,+l;S,A)+V(At,-l;S,A) v 1 

An experimental measurement of time-dependent CP violation at a l?-factory 
includes the following steps: 

1. Reconstruct one B decaying to fcp- 

2. Determine q using all available information on the B — > / tag decay. 

3. Reconstruct vertices for fcp and / tag and determine At from the distance 
between the two B vertices. 

4. Obtain S and A by fitting the At distribution of reconstructed signal candi- 
dates. 

Each step is described in more detail below. 

3.4. Measurement of <fii 

At the quark level neutral B meson decays into (cc)K° are induced by a b — > ccs 
transition. Since both leading and sub-leading order diagrams of this process contain 
neither V u b nor Vtd, there is no complex phase in the decay amplitude. Thus <po is 
zero and the mixing-induced CP violation parameter S is directly related to one of 
the CP- violating angles, (f>i. In the SM, 

<S = — rjcp ■ sin 2(f)\ and A ~ (3'5) 

are expected. 

3.4.1. B° —7- (cc)K° reconstruction 

We reconstruct J/ipKg, J/tpK®, i[)(2S)Kg, and xaKg as the fcp in neutral B 
meson decays to (cc)K°. J/ip mesons are reconstructed via their decay into oppo- 
sitely charged lepton pairs (e + e~ or [i + n~) while ip(2S) mesons are reconstructed 
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by lepton pairs as well as J/ifm + Tr~ final states. We reconstruct Xcl mesons in the 
mode and K$ mesons in the tt + tt~ final state. 

For B° — > J/ipKg, ip(2S)Kg, and XciKg candidates, the B signal is identified us- 
ing two kinematic variables calculated in the T(45) CM: the energy difference AE = 
E B —E^ egjn and the beam-energy constrained mass Mb c = v^(^bcam) 2 ~~ (p*b) 2 ' wnere 
-^beam * s the beam energy in the CM of the T(4S) resonance, and E* B and p* B are the 
CM energy and momentum of the reconstructed B candidate, respectively. In the 
B° — > J/ipK9 case, candidate K 9 mesons are selected using information recorded in 
the ECL and/or the KLM. Since the K L energy cannot be measured, we determine 
only its direction. Thus B° — > J/ipK L candidates are identified by the value of p* B 
calculated using a two-body decay kinematic assumption. 

The Mb c distribution for signal candidates with a stringent AE requirement 
(\AE\ < 40 MeV for J/i/jK^, \AE\ < 30 MeV for ip(2S)K%, and \AE\ < 25 MeV 
for XciKg) as well as the p* B distribution for J/ipK® candidates are shown in Fig. 7. 
The signal yields and purities are estimated for each fcp mode and given in Table II. 




Fig. 7. M bc distribution within the AE signal region for B° -s> J/ipK% (black), ip(2S)K° s (blue), 
and Xci Kg (magenta); the superimposed curve (red) shows the fit result for all modes combined 
(left) and the p* B distribution for B° — > J/ipK% candidates with the results of the fit separately 
shown as signal (open histogram), background with a real J ftp and a real K ( [ (yellow), back- 
ground with a real J/ip but without a real K% (green), and background without a real J/ijj 
(blue) (right). 



Table II. Signal yield (iVaig), CP eigenvalue (rjcp), and purity for each B° — >• fcp mode. 



B decay mode 


Vcp 




Purity (%) 


J/^K S 


-1 


12649 ± 114 


97 


i>{2S){l + t-)K% 


-1 


904 ± 31 


92 


iP{2S){J/iP-k + k-)K s 


-1 


1067 ± 33 


90 




-1 


940 ± 33 


86 




+1 


10040 ± 154 


63 
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3.4.2. Flavor tagging 

For the events in which we reconstructed B° — > fcp candidates, the neutral 
B flavor is identified from the decay products of the accompanying B meson. The 
available information is obtained from leptons, kaons, A baryons, and pions. Leptons 
directly coming from B decay and secondary leptons and strange particles in the 
cascade decays carry the mother 6-flavor information. Low momentum tagging pions 
may come from D* ± decays. In addition, there are high momentum pions originating 
from B° —7- D**) + 7r~ or Z)W p~ decays. Both types of tagging pions give some 
information about 6-flavor. The information from all the decay products is handled 
by a multi-dimensional likelihood approach with corresponding look-up tables. 17 ) 

To calibrate w, we select a flavor specific final state of neutral B meson decays 
such as semileptonic B° — > D* + £~i> decays and hadronic B° — > D(* + ^ir~ and D + * p~ 
decays. We then determine the wrong tag fraction w by measuring the time evolution 
of the opposite-sign flavor asymmetry, as it exhibits a At dependence proportional to 
(1 — 2w) cos (Am^At). We also determine Aw, which is the difference in w between 
q = +1 and —1 events. For B° — > J/ipK® decay, we obtain the effective tagging 
efficiency, e e d = e(l — 2w) 2 = (30.1 ± 0.4)%, where e is the tagging efficiency. 

3.4.3. At determination and its resolution 

In energy-asymmetric e + e~ collisions at KEKB, the T(4S) is produced with a 
Lorentz boost of j3j = 0.425 nearly along the z-axis, which is defined as the direction 
anti-parallel to the e + beam at Belle. Since B mesons are approximately at rest with 
respect to the T(AS), we can measure At by measuring the displacement between 
the two B meson decay vertices in the z direction, Az, 

& -IT-- ( 3 - 6 ) 

he 

The B meson decay vertex is reconstructed by a Lagrange multiplier approach, 
which minimizes the x 2 calculated from the decay vertex position and the daughter 
particle tracks. 18 ** We call this procedure a "vertex fit". The vertex fit is carried 
out using daughter tracks with a sufficient (minimal) number of SVD hits and a 
constraint on the interaction-region profile in the plane perpendicular to the beam 
axis. 

Because of the negligible flight length of J /ip or ip(2S) mesons, the vertex re- 
constructed from their daughter lepton tracks can represent the B° — >■ fcp decay 
vertex; its resolution is found to be approximately 75 /xm. On the other hand, the 
B° — > /tag vertex is obtained with well-reconstructed tracks that are not assigned to 
fcp- Here, high momentum leptons are always retained because they usually come 
directly from semileptonic B meson decays. Since /t ag may contain long-lived par- 
ticles such as D + , D°, Kg, and so on, the vertex reconstructed using the daughter 
tracks coming from these intermediate particles can deviate from the true B° — > / tag 
vertex. This effect is minimized by removing tracks that are identified by a large 
contribution to the vertex fit \ 2 - The / tag vertex position resolution is found to be 
approximately 165 ^im. 

In the Belle experiment, the contributions to At measurement error are divided 
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into three categories: detector measurement error, the effect of secondary particles 
in /tag vertex reconstruction, and the kinematical approximation, At ~ Az/((3jc). 
These three effects are convoluted on an event-by-event basis to obtain the At res- 
olution function, which is used in a maximum likelihood fit to extract S and A as 
discussed in the next section. 

3.4.4. Extracting CP violation parameters 

We determine sin20i and A from a maximum likelihood fit using At and q 
information obtained on an event-by-event basis from signal candidates. By taking 
the effect of incorrect flavor assignment into account, the probability density function 
(PDF) expected for the signal distribution is given by 



r sig (At) 

e -\At\/r B0 



1 - qAwi + q(l - 2wi) 



4t b o 

(-r) C p)sm2(f) 1 sm(Am d At) + Acos{Am d At) \. (3-7) 



The distribution is convoluted with the At resolution function R s \ g (At), which takes 
into account the finite vertex resolution as described in Sect. 3.4.3. The background 
PDF Vb)ig{At) is determined by the events found in a sideband region well away from 
the signal region in M^-AE space as well as Monte Carlo (MC) events. A small 
component of broad outliers in the Az distribution, caused by misreconstruction, 
is represented by a Gaussian function V \(At) with a ps 30 ps. We determine the 
following likelihood value for each event indexed by i: 

Vi(At i ,q i ;sm2(/) 1 ,A) 

/oo 
V^At^R^AU - At')d(At') (3-8) 
-oo 

+(1 - foiifbkgPbkgiAU) + UPoi(AU), 

where f a \ is the outlier fraction, / s i g and /bkg are the signal and background probabil- 
ities calculated as functions of AE and M^ c . The CP violation parameters, sin20i 
and A, are determined by maximizing the likelihood function 

L(sin 2#i,.A) =Y[ri(At i ,q i ;sm2(f) 1 ,A), (3-9) 

i 

where the product runs over all events. A fit to the candidate events results in the 
CP violation parameters, 19 ) 

sin 2^1 = 0.667 ± 0.023(stat) ± 0.012(syst), (3-10) 
A = 0.006 ± 0.016(stat) ± 0.012(syst). 

The background-subtracted At distribution for q = +1 and q = — 1 events and the 
asymmetry for events with good quality tags are shown in Fig. 8. The world average 
of sin 20i is now 0.68 ± 0.02, which is a firm SM reference. 
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Fig. 8. The background-subtracted At distribution for q = +1 (red) and q = — 1 (blue) events and 
asymmetry for events with good quality tags in (cc)K% (left) and J/%j)K% (right) decays. 



3.4.5. Search for new physics using CP violation measurements in b — > s penguin 
modes. 

B meson decays involving penguin diagrams are thought to be a sensitive probe 
for new physics (NP) beyond the SM because of the one-loop nature of penguins. 
NP could appear as deviations of CP violation parameters from the SM expectation. 
In this section, some highlight results for penguin modes are reviewed. 

In SM b — > sqq hadronic B decays, the relevant coupling is V* b Vt s and the weak 
phase is the same as in the b — > ccs transition, e.g. B° — > (cc)K° decay. Therefore, 
the main point is to check whether the penguin CP violation results deviate from the 
SM expectation, S = — rjcp sin 20i and A = 0. In this context, the time-dependent 
CP-violating parameters are denoted as sin 2cj)f and A. The modes B° -)■ 4>K°, 
n'K°, and K°K°K° that involve only b — > sss processes are of special interest, since 
the SM theoretical uncertainty for CP violation is small for these decay processes. 

In the Belle experiment, attempts to perform measurements of time-dependent 
CP violation in b — >■ sqq induced decays with B° — > n'K^ and 4>K^ modes were 
made from the earliest stage of data taking, starting in 2002. In 2003, using a 
152 x 10 6 BB data sample, the value of S in B° -)• 4>K° S flipped sign and exhibited 
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a 3.5a deviation from the S parameter measured in B° — > (cc)K° modes. 20 ) This 
was very striking and suggestive of an NP effect. In 2006, with a larger statistics 
data sample corresponding to 535 x 10 6 BB, updated measurements were reported. 
These measurements added B° — > rj'K® and (j)K® decays to the B° — > tj'K and 
<f>K° sample. 21 ^ The results are summarized in Table III. In B° — > rj'K decay, CP 

Table III. Measurements of CP violation parameters, s'm2<f>1 a and A, in B° — > rj'K , (f>K°, and 
K%K%K% modes with a 535 x 10 6 BB data sample. The first and second errors are statistical 
and systematic errors, respectively. 



B decay mode 


sin 20°!" 


A 




+0.64 ±0.10 ±0.04 


-0.01 ±0.07 ±0.05 


4>K° 


±0.50 ±0.21 ±0.06 


±0.07 ±0.15 ±0.05 


K%K%K% 


±0.30 ± 0.32 ± 0.08 


-0.31 ±0.20 ±0.07 



violation is observed with a statistical significance of 5.6<t. In all these three B decay 
modes, the large deviation from B° — > (cc)K° has disappeared. 

In spite of the small theoretical uncertainty, experimentally, several contribu- 
tions overlap in B° — > (pK° because of the relatively wide natural widths of the 
resonances that contribute in the K + K~ final state. In order to resolve these inter- 
fering contributions, Belle fits the time-dependent Dalitz distribution by expressing 
each contribution at the amplitude level for the B° — > K + K~Kg candidate events. 
With this technique, the extracted parameter is not sin2</>| ff but rather the angle 
4>l s itself and A. Therefore the result does not have a two-fold ambiguity between 
<f>l s and 7r/2 — </>f ff . In B° — > K + K~K$ decays, we find four solutions related to 
resonant amplitude interference. The preferred one is identified using external in- 
formation related to /o(980) and fx (assumed to be /o(1500)) branching fractions. 
The obtained CP violation parameters are summarized in Table IV. 22 ) These are 

Table IV. CP violation parameters in B° — > K + K~ K% time-dependent Dalitz analysis, 0f ff and 
A. The first, second, and third errors are statistical, experimental systematic, and Dalitz model 
uncertainties, respectively. 



B decay mode 


4>f 


A 


<PK» S 


(32.2 ± 9.0 ± 2.6 ±1.4)° 


±0.04 ±0.20 ±0.10 ±0.02 




(31.3 ±9.0 ±3.4 ±4.0)° 


-0.30 ±0.29 ±0.11 ±0.09 



consistent with the CP violation in B° — > ccK° decays at the la level. 

Including other b — > s mediated B decays, the precision of sin 2<pf is still sta- 
tistically limited, typically 0.1 ~ 0.2. Obtaining O(10~ 2 ) sensitivity requires an 
integrated luminosity of O(10 ab -1 ), and a Super B-factory experiment. 

3.5. Measurement of 4>2 

After the first observation of CP violation in B meson decays, which gave a 
measurement of (f>i, a precise measurement of 4>2 became the next target of CP 
violation measurements for the validation of the Kobayashi-Maskawa model. The 
first Belle measurement of CP asymmetry parameters in B° — > tt + tt~ decay 23 ) was 
reported in March 2002, representing the second decay mode (after B — > ccK°) with 
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a time-dependent CP violation measurement. 

The decay modes used for measurements are those proceeding via b — >■ u 
transition, such as B° — > it + it~, B° — > p + 7r~, B° — > p + p~ ■ The b — > u transition is 
shown in Fig. 9(left) and includes the Cabibbo-Kobayashi-Maskawa (CKM) element, 
V u h\ it can be shown that the time dependent CP asymmetry is then given as 
S = sin 2^2 and A 2± 0. However, an additional amplitude, a "penguin diagram" 
(Fig. 9(right)), contributes and has a phase that is different from the tree diagram 
(V t d instead of V u b). This causes a deviation of S from sin 202 and a non-zero A. 




The first measurement was attempted using the B — > tt + tt~ decay mode. 
This decay has the simplest two-body topology and was one of the first well estab- 
lished charmless B decays. The reconstruction of the decay is straightforward: a 
pair of oppositely charged pions with an invariant mass consistent with the l?-meson 
mass (Mb c = mg) is selected; the B meson energy in CM is required to be consistent 
with the beam energy {AE = 0). However, the selected sample suffers from a very 
large background from the e + e~ — > qq (q = u,d,s,c) continuum process since the 
same kinematic properties can easily be faked by two oppositely charged pions frag- 
mented from primary quarks and carrying about half of their momentum. Another 
significant background is from B° — > K + it~ decay, where the kaon is misidentified 
as a pion. The branching fraction for the former decay mode is about four times 
higher than that of B° — > 7r + 7r~. In this case, the reconstructed AE is shifted by 
—40 MeV and good K/ir separation and good momentum resolution are important 
to reduce this background. 

The continuum background is suppressed utilizing a difference in the global event 
topology for the two classes of events; continuum events have a two-jet like shape 
while BB events have an isotropic shape as the two B mesons are produced almost 
at rest in the CM. To quantify the event shape, we use a Fisher discriminant 24 ) com- 
bining modified Fox-Wolfram moments. 25 ) We form a likelihood C s for signal 
(continuum background) using the Fisher discriminant and the angle between the 
flight direction of the B candidate and the beam direction in the CM, cos 0b- The 
likelihood ratio TZ = C s / (C s + £5) is used as the final continuum suppression param- 
eter. In the early analyses, 23 ) we imposed a tight requirement on TZ by optimizing 
S/ \/ S + B, where S and B are the expected number of signal and background events, 
respectively. In a later analysis, 26 ) we optimized the TZ requirement depending on 
the flavor tagging quality. 

The AE distribution of B° — > it + it~ candidates is shown in Fig. 10. Background 
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events due to three-body decays populate the negative AE region but they do not 
contribute in the B° — > ir + ir~ signal region (\AE\ < 0.064 GeV). 
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§200 
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71+71" 
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qq 

3-body 



</>300 



■0.2 0.2 0.4 
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Fig. 10. AE distribution of B° ->■ tt+tt" can- 
didates. In order to enhance the signal, re- 
quirements are imposed on the two other 
variables, Mbc and 72.. 




Fig. 11. At distribution for B° and B° tagged 
B — > ir + n~ events (top) and the CP asym- 
metry together with the fit result (bottom). 



The vertex reconstruction and the flavor tagging are performed in the same 
way as for the sin2</>i measurements. The CP violation parameters are extracted 
from a fit to the At distribution for the events in the signal region in AE and Mb c 
([5.271, 5.287] GeV/c 2 ). The PDFs include the signal, continuum background, and 
B° — > K + ir~ background. The first result was reported using 48 x 10 6 BB pairs: 23 ) 



Stttc — — 1.21^27(stat)^Qj3(syst) 
A** = +0.941°;^ (stat) ± 0.09(syst) 



(3-11) 



In the latest analysis using 535 x 10 6 BB pairs, 26 ' a stringent selection on K/tt par- 
ticle identification is not imposed and instead the B° — > K + ir~ decays are included 
as a component in the fit to extract the CP violation parameters. This increases 
the signal detection efficiency by 23% and improves the measurement errors by 10%. 
The results are 26 ) 



S n7T = -0.61 ± O.lO(stat) ± 0.04(syst) 
Attt = +0.55 ± 0.08(stat) ± 0.05(syst). 



(3-12) 



The At distribution and the asymmetry together with fit results are shown in 
Fig. 11. A clear non-zero A^m i*^* & clear direct CP violation, is seen (the asymmetry 
exhibits a significant cosine term). As shown above, the first measurement already 
indicated CP violation in decays with a significance of 2.9a. The first evidence 
of direct CP violation in a B decay mode was reported with 3.2cr significance in 
January 2004 using a sample of 152 x 10 6 BB pairs , 27 ) Although this claim was not 
widely accepted at that time because the result of the BaBar collaboration showed 
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Fig. 12. 1 — C.L. for a range of <f>2 values as ob- ^ 2 1 9 ' 

tained with an isospin analysis of B — > nn 

decays. The solid and dashed lines indicate Fig. 13. 1-C.L. as a function of <f) 2 from the 
C.L. = 68.3% and 95%, respectively. average of the Belle and BaBar results of 

B° -¥ n+Tv-, B° -> pn, B° -»■ pp. 



a rather small Ann value, the latest world average, A n7T = 0.38 ± 0.06, 28 ) establishes 
CP violation in B° — > Tr + ir~ decays with a significance above ha. 

The large direct CP violation indicates that the contribution of the penguin 
diagram is sizable and the deviation of S nn from sin 2^2 may be significant. The 
angle 02 can be extracted using the isospin relation among branching fractions and 
CP asymmetries of B° — > tt + it~, B — > 7r°7r°, and B + — > tt + tt° decays; this was 
first proposed by M. Gronau and D. London. 29 ) The result, shown in Fig. 12, is 
obtained using the results for S nn and given above and the world average values 
of branching fractions of the three B — > tttt modes and direct CP asymmetry in 
B° —7- 7r°7r°. Using this method, there are multiple discrete ambiguities for the angle 
4>2- The solution that is closest to the global fit result 30 - 1 gives (p2 = (97 ± 11)°. 

The final state in B° — > p + ir~ decay is not a CP eigenstate, but the decay 
proceeds through the same quark diagrams as B° — > ir + ir~. Since B° and B° can 
decay to p + vr~, time-dependent CP violation can occur and provide information on 
4>2- Here the final state is 7r + 7r~7r° and the decay B° — > 7r + 7r _ 7r° contains three 
intermediate states; B° — > p + vr~, p _ vr + , and p°ir°. These three amplitudes interfere 
and their magnitudes and relative strong phases can be extracted from a Dalitz plot 
amplitude analysis. Knowing the hadronic phases of these amplitudes in the Dalitz 
plane, a time-dependent Dalitz plot analysis allows the determination of 4>2- This 
method provides 4>2 without ambiguities (assuming large signal statistics) except for 

4>2 — > 4>2 + 7T- 

The reconstruction and continuum suppression are similar to the B° — > tt + tt~ 
analysis with an additional 7r° reconstructed in the tt° — > 77 decay mode. CP 
violation parameters are obtained from a three-dimensional fit to the distribution of 
At and two Dalitz distribution parameters, M 2 + and M 2 _ . Belle performed the 
analysis using 449 xlO 6 BB pairs. ' The amplitudes include p(770) and higher mass 



20 



J. Brodzicka et al. 



resonances, p(1450) and p(1700). The time-dependent Dalitz plot distribution is 
parameterized with 27 real parameters describing the components that have different 
time- and Dalitz plot behaviors. CP violation parameters for B° — > p^vr^, B° 
7r°7r° decays and 02 are extracted from these parameters. We obtain a 68° < 02 < 95° 
at a 68.3% confidence level (C.L.) interval for the solution consistent with the global 
fit result. A large region (0° < 2 < 5°, 23° < 2 < 34°, and 109° < 2 < 180°) 
also remains. With a larger data sample, a more restrictive constraint without 
ambiguities is expected from this measurement. 

In the B° — > p + p~ mode a pseudoscalar decays into two vector particles and 
the final state is a mixture of CP-even and CP-odd amplitudes. In order to extract 
the fraction of each CP component, an angular analysis is required. Fortunately, 
the fraction of the longitudinal polarization turns out to be close to 100%, 33 )~ 35 ) 
simplifying the measurement. The signal candidates are reconstructed in p ± — > tt^tt 
decays. Because of two 7r°s in the final state, the combinatorial background due to 
fake 7T° candidates is very large. The results using 535 x 10 6 BB pairs are: 36 ) 

A p+p - = +0.16 ± 0.21(stat) ± 0.07(syst) 

S p+p - = +0.19 ± 0.30(stat) ± 0.07(syst) (3-13) 

In this mode, 02 can be obtained using an isospin relation similar to that in the 
B° —7- 7r + 7r~ case. Because the branching fraction for B° — > p°p° is much smaller 
than those of B° — > p + p~ and B + — > p + p°, the deviation of 02 from the measured 
value is small and some ambiguities are degenerate. So far only an upper limit on 
B(B° — > p°p°) has been obtained; this is used in the isospin analysis. The isospin 
analysis gives 62° < 2 < 106° at the 68.3% C.L. 

All of the above results and results from the BaBar collaboration can be com- 
bined to obtain the 02 constraint shown in Fig. 13. 30 ) 02 = (89.0^42)° is obtained 
at a 68.3% C.L. 

3.6. Measurement of 03 

The angle 0i has been now measured with high precision (Sect. 3.4). Mea- 
surement of the angle 2 is more difficult due to theoretical uncertainties from the 
contributions of penguin diagrams (Sect. 3.5). Precise determination of the third 
angle of the unitarity triangle, 03, is possible using B ± — > DK ± decays. However, 
it requires much more data than determinations of the other angles. The determi- 
nation of 03 is theoretically clean due to the absence of loop contributions; 03 can 
be determined using tree-level processes only, exploiting the interference between 
b — > cud and b — > ucd transitions that occurs when a process involves a neutral D 
meson reconstructed in a final state accessible to both D° and D decays. Therefore, 
03 provides an SM benchmark, and its precise measurement is crucial in order to 
disentangle non-SM contributions to other processes, via global CKM fits. 

Several different D decays have been studied in order to maximize the sensitivity 
to 03. The archetype is the use of D decays to CP eigenstates, a method proposed 
by M. Gronau, D. London, and D. Wyler (and called the GLW method). 37 ) Belle 
makes use of CP-even modes (Di), such as K + K~, and CP-odd modes (-D2), such 
as Kg7T°. To extract 03 using the GLW method, the following observables sensitive 
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to CP violation are used: the asymmetries 

B(B~ -»• D li2 K~) - B{B+ -»• D h2 K + ) 2r B smS B sin, 



^1.2 = 



6(5- -»• Di 2 if-) + -> £>i 2^+) 1 + r\ ± 2r B cos «J B cos 03 

(3-14) 

and the ratios 

B(B-^DOK~)+B(B+^D°K+) = 1 + ± 2r B cob cos (3-15) 

where re is the ratio of the magnitudes of the two tree diagrams shown in Fig. 14 
and 5 B is their strong-phase difference. The value of r B is given by the product of the 
ratio of the CKM matrix elements | V^ 6 V^ S |/| V^V^ S | ~ 0.38 and the color suppression 
factor, which altogether results in a value of around 0.1. In the expressions above, 
mixing and CP violation in the neutral D meson system are neglected. Among these 




u b fc . » u 

K" "4 D' 

s % c 

b p. » C „ S 

B . . D° . K 

u ^ u u ^ u 

Fig. 14. Feynman diagrams for B~ — > D°K~ and B~ — > D K~ . 

four observables, 1Zi$ and Aip, only three are independent (since AilZi = —A-fit'i)- 
Recently, Belle updated their GLW 

analysis using their final data sample of Table y Resultg of the GLW anal ig for 

772xlO b BB pairs. ' The analysis uses DK^ mode. 

D° decays to K + K~ and tt + tt^ as CP- 1 03 ± 07 ± 03 

even modes (Fig. 15), Kgir and Kgi] as TZ 2 1.13 ± 0.09 ± 0.05 

CP-odd modes. From Eqs. 3-14-3-15, Ai +0.29 ± 0.06 ± 0.02 

the signs of the A\ and A2 asymmetries -0.12 ± 0.06 ± 0.01 

should be opposite, which is confirmed 

by experiment (Table V). 

The difficulties in the application of the GLW methods arise primarily due to the 
small magnitude of the CP asymmetry of the B^ — > DcpK^ decay, which may lead 
to significant systematic uncertainties in the observation of the CP violation. An 
alternative approach was proposed by D. Atwood, I. Dunietz, and A. Soni. 39 ) Instead 
of using D° decays to CP eigenstates, the ADS method uses Cabibbo-favored and 
doubly Cabibbo-suppressed decays: D — > K~ir + and D — > K~n + . In the decays 
B + [K~ir + ]£)K + and B~ — > [K + tt~]dK~ , the suppressed B decay is followed by 
a Cabibbo-allowed D° decay, and vice versa. Therefore, the interfering amplitudes 
are of similar magnitude, and one can expect a large CP asymmetry. Unfortunately, 
the branching ratios of the decays mentioned above are small. The observable that 
is measured in the ADS method is the ratio of the suppressed and allowed branching 
fractions: 

B(B ± -> [K^it^dK*) 22 , 
■Kads = B ^ B± — ^ k ±^ dK ±^ = r B + r D + 2r B r D cos <fe cos 8, (3-16) 
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A E (GeV) A E (GeV) 

Fig. 15. Signals for B — > D±K ± decays. The left (right) figure is for B~ (B + ) decays. The plotted 
variable, AE, peaks at zero for signal decays, while background from B -+ Drr* appears as a 
satellite peak at positive values. 



and 

A B ( B ~ ~t [K + v-]pK-) - B{B+ -» [K-ir+] D K+) . . 

AADS = B(B- ^ [K+v-] D K-) + B{B+ -> [K-«+] D K+) = sm^ sm^ADS, 

(3-17) 

where rjj is the ratio of the doubly Cabibbo-suppressed and Cabibbo-allowed D° 
decay amplitudes and 5 is the sum of strong phase differences in B and D decays: 
S = S B + S D . The ADS analysis 40 ) using the full T(4S) data sample was reported 
by the Belle collaboration (Fig. 16). The analysis uses B ± -> DK± decays with D° 
decaying to K + tt~ and K~tt + (and their charge-conjugated partners). The signal 
yield obtained is 56^ J| events, which corresponds to the first evidence of an ADS 
signal (with a significance of 4.1c); the ratio of the suppressed and allowed modes 
is summarized in Table VI. Although the analyses with B^ — > DK^ decays give 
the most precise results, different B decays have also been studied. The use of two 
additional decay modes, D* — > Dtt° and D* — > D7, provides an extra handle on the 
extraction of 4>3 from B^ — > D*K^, which is becoming visible in the most recent 
results. 38 ^ 



Table VI. Results of the Belle ADS analyses. 



Mode 




72-ads 


-4ads 


B -> DK 

B -> D*K, D* - 

B -> D*K, D* - 


+ Z?7r° 
-> D-y 


n niRq+°- U044 + u - U0U7 

U.U100_q .0041 -0.0013 

n nin+ 008 + - 001 

U.U1U_ .007-0.002 

0.036^:^2 ± °-° 02 


n OQ+U.I26 4-U.U4 
U.O»_ .28-0.03 
, n 4 +l.l+0.2 
+ u -^-0. 7-0.1 

-0.5llo;29 ±0.08 



A Dalitz plot analysis of a three-body D meson final state allows one to obtain 
all the information required for determination of ^3 in a single decay mode. Three- 
body final states such as K®ir + 7r~ have been suggested as promising modes 41 ) for 
the extraction of 03 . As in the GLW and ADS methods, the two amplitudes interfere 
if the D° and D mesons decay into the same final state Kg7r + 7r~ . Assuming no CP 
asymmetry in neutral D decays, the amplitude for B + — > D[Kstt + tt^]K + decay as 
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Fig. 16. Signal for B — > DK ± decays from Belle ADS analysis. In these AE and NB (contin- 
uum suppression variable) distributions, [K + -k~]dK~ components are shown by thicker dashed 
curves (red). 

a function of Dalitz plot variables m 2 , = m 2 n , and m 2 = m 2 a is 

f B + = fD(ml,ml)+r B e^^ B f D (ml,m 2 + ) (3-18) 

where /n(ml,m?.) is the amplitude of the D° — > Kgir + 7r~ decay. Similarly, the 
amplitude for B~ — > D[KsiT + ir^]K^ decay is 

f B - = f D {ml,m\) +r B e- i ^ +iSB f D (ml,m 2 __). (3-19) 

The W K^tt+tt- decay amplitude fr> can be determined from a large sample 
of flavor-tagged 5° — > Kgir+ir^ decays produced in continuum e + e~ annihilation. 
Once fr> is known, a simultaneous fit to B + and B~ data allows the contributions 
of r B , 03 and 5b to be separated. The method has only two-fold ambiguity: (03, 8b) 
and (03 + 180°, 5 B + 180°) solutions cannot be distinguished. To test the consistency 
of the fit, the same procedure was applied to the B^ — > D^n^ control samples 
and the B^ — > D^'K^ signal. A combined unbinned maximum likelihood fit to the 
B + and B~ samples with free parameters r^, 03, and 5b yields the values given 
in Table VII. Combining B ± -> DK± and B ± Z?*iv" ± , we obtain 42 ) the value 
03 = (78+}i ± 4 ± 9)°, where the sources of uncertainties are statistical, systematic, 
and due to imperfect knowledge of the amplitude model that describes D — > Kgir + ir~ 
decays. The last source of uncertainty can be eliminated by binning the Dalitz plot 

Table VII. Results of Belle Dalitz plot analyses. 

Mode cj> 3 (°) 5 B Q r B 

B -> DK ±5 137+ig ± 4 0.16 ± 0.04 ± 0.01 

B -> D*K 74+^ ± 4 342±^ ± 3 0.20 ± 0.07 ± 0.01 

(Refs. 41 )' 43 -*), using information on the average strong phase difference between D° 
and D° decays in each bin that can be determined using quantum correlated ^(3770) 
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data. Results have been published recently by CLEO-c. 44 -* The measured strong 
phase difference is used to obtain a model-independent result: 45 ) 

03 = (77 ± 15±4±4)°, (3-20) 

where the last uncertainty is due to the statistical precision of the CLEO-c results. 

§4. Measurement of \V c b\ and semileptonic, and leptonic B decays 

4.1. Introduction 

The Cabibbo-Kobayashi-Maskawa (CKM) 
matrix elements \V c b\ and \V u b\ are deter- 
mined from semileptonic B — > Xlv (£ = 
e, n) decays to charmed and charmless fi- 
nal states, respectively (Fig. 17). These 
decays are chosen because semileptonic de- 
cays proceed via leading-order weak in- 
teractions and thus are free of possible 
non-Standard Model contributions. Their 
branching fractions are sizable compared to 

purely leptonic B -»• Iv decays, and have Fi S- 17 - Illustration of the semileptonic 
, , , ■ , ■ , i n B meson decay B —¥ Xlv. 

hadromc uncertainties that are well con- 
trolled by various theoretical techniques. 

In this section, we also discuss purely leptonic and semileptonic B decays involv- 
ing a heavy r lepton. At present these decays are not relevant for the determination 
of \V c b\ and \V u b\ but are studied because of their sensitivity to the charged Higgs 
boson and other manifestations of new physics. 

There are two orthogonal approaches to measuring semileptonic decays and de- 
termining | V c b | and \V u b\: Analyses can either be exclusive, i.e., these reconstruct 
only a specific semileptonic final state, such as D*£v, ttIu, .... Alternatively, the 
analysis can be inclusive, which means that it is sensitive to all semileptonic final 
states, X c lv or X u £v, in a given region of phase space, where X c and X u refer to a 
hadronic system with charm or without charm, respectively. 

Exclusive and inclusive analyses are affected by different experimental uncer- 
tainties. In addition, different and largely independent theoretical approaches are 
used to describe the QCD contributions in exclusive and inclusive decays. Since both 
approaches rely on different experimental techniques and involve different theoret- 
ical approximations, they complement each other and provide largely independent 
determinations of comparable accuracy for \V c b\ and \V u b\. This in turn provides a 
crucial cross check of the methods and our understanding of semileptonic B decays 
in general. 

4.2. \V cb \ 

4.2.1. \V c b\ from exclusive semileptonic decays 

The determination of \V c b\ from exclusive decays is based on the B — > D*lv or 
B —7- Dlv decay modes. Experimentally, one has to measure the differential decay 
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rate as a function of the velocity transfer w, denned as 

w = = , (4-1) 

m B m D (*) 2m B m DM 

where tub and m D {,) are the masses of the B and the charmed mesons, Pg and 
P D (*) are their four-momenta, and q 2 = (Pi + P v ) 2 . The point w = 1 is referred to 
as zero recoil, because there the charmed meson is at rest in the B meson frame. To 
determine \V c b\, the experimental analyses extrapolate the decay rate to w = 1, as 
theory can determine the decay form factors with greater accuracy at this kinematic 
point. When neglecting the lepton mass, i.e., considering only electrons and muons, 
the differential decay rate of B — > D*lv as a function of w is given by 46 ) 

S = " m D ,fV^lx{w)T\ W )\V ch \ 2 . (4-2) 

Here, G F is Fermi's constant equal to (1.16637 ± 0.00001) x 10~ 5 GeV" 2 and X (w) 
is a known phase space factor, 

„2" 



x( w ) = ( w + iy 



w 1 — 2wr + r 

1 + 4- 



w + 1 (1 — r 



,2 



(4-3) 



where r = mo* / 'tub- The dynamics of the decay are contained in the form factor 
J-(w), which can be parameterized by the normalization ^(1), the slope p D *, and 
the amplitude ratios i?i(l) and i?2(l) in the framework of the heavy quark effective 
theory (HQET). 47 ) 

A similar expression can be derived for the differential rate of the decay B — > 
D£v, 

d £ = ^? (mB + mD)2{w2 - i) 37 ^ 2 wi^i 2 • ( 44 ) 

As the D meson is a pseudoscalar, the form factor Q(w) of this decay is simpler than 
F(w) and can be parameterized by the normalization £7(1) and the slope p 2 D only. 47 ) 

In the limit of infinite quark masses, the form factors F(w) and Q(w) coincide 
with the Isgur-Wise function, 48 ) which is normalized to unity at zero recoil, w = 1. 
Corrections to the heavy quark limit have been calculated in the framework of lattice 
QCD (LQCD). In LQCD, the QCD action is discretized on a Euclidean spacetime 
lattice and calculations are performed numerically on computers using Monte Carlo 
methods. Physical results are then recovered in the limit of zero lattice spacing. 
Because lattice results are obtained from QCD first principles, they can be improved 
to arbitrary precision, given sufficient computing resources. 

The form factor values at w = 1 are the main theoretical input needed for the 
determination of |Vc&| from exclusive decays and also the main source of theoretical 
uncertainty. The current LQCD value of ^"(1), describing the decays B — > D*£v, 
is 49 ) 

= 0.908 ± 0.017 . (4-5) 
The LQCD B — > Div form factor at zero recoil is calculated to be 50 ) 

0(1) = 1.074 ± 0.024 . (4-6) 
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The Belle measurement of B° ->■ D*^£ + u^ is based on 772 x 10 6 BB events, 
resulting in about 120000 reconstructed decays. In this analysis the decay chain 
D*~ —7- D°it~ followed by D° — > K + ir~ is reconstructed and D* candidates are 
combined with a charged lepton I (I = e, fi) with momentum between 0.8 GeV 
and 2.4 GeV. As the analysis makes no requirement on the second B meson in the 
event, the direction of the neutrino is not precisely known. However, using the 
cos 9 by variable, 

2E B E Y - m 2 B - ml 

cose BY = ^ ' ^) 

with Y = D*l, the B momentum vector is constrained to a cone centered on the 
D*£ direction. By averaging over the possible B directions one can approximate 
the neutrino momentum and calculate the kinematic variables of the decay (w and 
three decay angles) . The parameters of the form factor F{w) are obtained by fitting 
these four kinematic distributions. The very large data sample led to much reduced 
statistical and systematic uncertainties. 

The result of the Belle analysis (after rescaling input parameters to their most 
recent values 15 )) is 

F(l)\V cb \ = (34.7 ± 0.2(stat) ± l.O(syst)) x 10 -3 , (4-8) 

where the dominant systematic uncertainties stem from charged track reconstruction. 
Assuming the form factor normalization of Eq. 4-5, we obtain 

\V cb \ = (38.2 ± l.l(exp) ± 0.7(th)) x 10~ 3 . (4-9) 

The experimental uncertainty is at the level of 3.0% while the theoretical uncertainty 
from lattice QCD amounts to 1.9%. 

In addition, the decay B — > Dlv has been studied at Belle using 10.8 x 10 6 
BB events. 52 ^ For the determination of \V c b\, the decay B — > D*lv is preferred over 
B — > Dlv for both theoretical and experimental reasons: On the theory side, the rate 
at zero recoil is lower for B — > Dlv than for B — > D*£u due to the factor (w 2 — l) 3 / 2 
in the expression for the width (instead of Vw 2 — 1 in the D* case). Experimentally, 
due to the presence of the slow pion in the decay D* — > Dir, the D* signal is cleaner 
than the D signal and backgrounds in the analysis of B — > Dlv are typically the 
limiting factor. 

The result of the Belle analysis (after rescaling input parameters to their most 
recent values 15 )) is 

Q(\)\V ch \ = (40.8 ± 4.4(stat) ± 5.2(syst)) x 10" 3 , (4-10) 

with the dominant systematic uncertainties from background estimation. Assuming 
the G(l) value from Eq. 4-6, we obtain 

\V cb \ = (38.0 ± 6.3(exp) ± 0.8(th)) x 10~ 3 . (4-11) 

This determination of \V c b\ is consistent with the B — > D*£u value but has a signifi- 
cantly larger uncertainty. 
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4.2.2. \V c b\ from inclusive semileptonic decays 

The theoretical tool for calculating the inclusive semileptonic decay width r(B — > 
X c tv) of the B meson is the operator product expansion (OPE). In this framework, 
a simplified form reads 53 )' 54 ) 

r( B -> X c t„) - ^4^(1 + «M <?><"> + «<"> «*><"> + 0(^)) . (442) 
1927r d v m% ml ' 

where the expansion parameter is the 6-quark mass nib- At leading order in 1/mb, the 
OPE result coincides with the parton model, i.e., with the decay width of a (hypo- 
thetical) free 6-quark. Corrections to the free 6-quark decay arise at order 1/m?: the 
term (Os)( / u) denotes the expectation values of local dimension 5 operators, which 
depend on the renormalization scale fi. A detailed analysis shows that only two op- 
erators appear at 0(l/m?): the kinetic operator, related to the kinetic energy of the 
6-quark inside the B hadron, and the chromomagnetic operator, related to the B*-B 
hyperfine mass splitting. At C(l/m 3 ), new operators appear. These expectation 
values of local operators describe basic hadronic properties of the B meson and do 
not depend on the observable (here r(B — > X c tv)) calculated using the OPE. As 
they contain soft hadronic physics, they cannot be calculated by perturbative QCD. 

These matrix elements are multiplied by the Wilson coefficients c$, c$, . . . , which 
encode the short-distance QCD contributions to the process and thus can be calcu- 
lated in perturbation theory as a series in powers of a s . Hence, the OPE factorizes the 
calculable and the non-calculable contributions to the semileptonic width. Even more 
interestingly, the hadronic matrix elements in the non-calculable part also appear 
in similar OPE expressions for other inclusive observables in semileptonic B decays. 
By measuring these additional observables, one can determine the non-perturbative 
OPE parameters, substitute them into the expression of the semileptonic width, and 
measure \V c b\ with a total precision of about 1-2%. This is the basic idea underlying 
the global fit analysis of \V c b\ discussed in this section. 

These other observables are the (truncated) moments of the lepton energy Eg 
(in the B rest frame) and the hadronic mass squared m\ spectra in B — > Xiv. The 
quantity m 2 x is the invariant mass squared of the hadronic system X c accompanying 
the lepton-neutrino pair. The lepton energy moments are defined as 

\^e)E cut - ^-t^ — r , (4-13) 

where E cut is the lower lepton energy threshold and 

R n (E cut ) = [ Ef^dEi . (4-14) 
JE e >E cut a^e 

Here, dr/dEg is the partial semileptonic width as a function of the lepton energy. 
The hadronic mass moments are 

/ 2n\ S n (E cut ) 

(m x )E cut = q f „ ; , (4-15) 
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with 

Sn(E cut ) = [ m f^dm 2 x . (4-16) 

JE e >E cut am x 

Here the integration over the B — > X c lv phase space is restricted by the require- 
ment Eg > E cut . These observables can be expanded in OPEs similar to Eq. 4-12, 
containing the same non-perturbative parameters. 

In practice, the semileptonic width and moments in B — > X c tv have been calcu- 
lated in two theoretical frameworks, referred to by the name of the renormalization 
scheme used for the quark masses (though this is not the only difference in the 
calculations): The calculations in the kinetic scheme are now available at next- 
to-next-to-leading (NNLO) order in a s . 53 )' 55 ) At leading order in the OPE, the 
non-perturbative parameters are the quark masses rrn, and m c . At 0(1/ ml) the 
parameters are p% and /j,q, and at 0(l/m\) the parameters p\, and p\ s appear. 
Independent expressions have been obtained in the IS scheme. 54 ^ Here, the long- 
distance parameters are mj, at leading order, Ai and A2 at 0(1/171%) and pi, 7~i_3 at 
0(l/mf). Note that the numerical values of the quark masses in the two schemes 
cannot be compared directly due to their different definitions. 

Belle has measured moments of inclusive observables in B — > X c tv decays. 56 -*' 57 ) 
The lepton energy Eg and hadronic mass squared m x spectra in B — > X c £u are 
based on 152 x 10 6 T(4S) — > BB events. These analyses first fully reconstruct 
the decay of one B meson (-Btag) in the event in a hadronic mode (or a hadronic 
tag). The tracks and clusters associated with -B tag are removed from the event. The 
semileptonic decay of the second B meson in the event (-B s ig) is then identified by 
searching for a charged lepton among the remaining particles in the event. In the 
lepton energy analysis, 56 ) the electron momentum spectrum p* in the B meson rest 
frame is measured down to 0.4 GeV/c. In the hadronic mass study, 57 ) all remaining 
particles in the event, after excluding the charged lepton (either an electron or muon), 
are combined to reconstruct the hadronic X system. The m x spectrum is measured 
for lepton energies above 0.7 GeV in the B meson rest frame. 

The observed spectra are distorted by resolution and acceptance effects and 
cannot be used directly to obtain the moments. In the Belle analyses, acceptance 
and finite resolution effects are corrected by unfolding the observed spectra using the 
singular value decomposition (SVD) algorithm. 58 ) Belle measures the lepton energy 
moments (Eg) for k = 0,1,2,3,4 and minimum lepton energies ranging from 0.4 
to 2.0 GeV. Moments of the hadronic mass (m k x ) are measured for k = 2,4 and 
minimum lepton energies between 0.7 and 1.9 GeV. 

To determine \ V c b\, Belle performs fits 59 ) to 14 lepton energy moments, 7 hadronic 
mass moments and 4 moments of the photon energy spectrum in B — > X s ~j based 
on OPE expressions derived in the kinetic 53 ^' 55 )' 60 ) and IS schemes 54 ) (Fig. 18). 
Both theoretical frameworks are considered independently and yield very consistent 
results: The fit to the Belle data in the kinetic scheme yields 

\V cb \ = (41.58 ± 0.90) x 10~ 3 , (4-17) 

while in the IS scheme we obtain 

\V cb \ = (41.56 ± 0.68) x 10~ 3 . (4-18) 
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While the result in the IS scheme is more precise (1.6% uncertainty compared to 
2.2% in the kinetic scheme), it should be noted that the assumptions on the dominant 
theory error are significantly different. 
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Fig. 18. Global fit of moments in B — > X c lv decays measured at Belle to theoretical expressions 
obtained in the kinetic scheme. The error bars show the experimental uncertainties. The error 
bands represent the theory error. Filled circles are data points used in the fit, and open circles 
are unused measurements. 



4.3. | F ub | 

4.3.1. \V u b\ from exclusive B — > X u tv decays 

The absolute value of V u b, one of the least known CKM elements, can be deter- 
mined from rate measurements of exclusive charmless semileptonic decays, such as 
B — > tt£v, B — > ptv and B — > ujiv. Of these, B° — > tv~£ + i/ decay has been the most 
extensively studied both theoretically and experimentally. The decay rates and \ V u b\ 
are related as 

dr(B° -> n-t+v) _ Gp 2 3 , 2 

aW 2 ~ 24^ Kbl PM+iQ )] ' ( } 

where Gf is the Fermi coupling constant, and f+{q 2 ) is the B — > it transition form 
factor, which is calculated in lattice QCD and by QCD sum rules. Compared to the 
inclusive measurements, described below, the exclusive measurements are relatively 
straightforward experimentally, but suffer from large theoretical uncertainties in the 
form factors, which must be determined from non-perturbative QCD calculations. 

The Belle collaboration has measured B — > n/ p/uilv decays. 61 )~ 63 ) The most 
recent measurement of the B° — > Tt~~l + u decay 63 - 1 uses a data sample containing 
657 x 10 6 BB pairs, and has the best precision for the \V u b\ determination. In 
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this analysis, signals are reconstructed by combining an oppositely charged pion 
and lepton (either electron or muon), originating from a common vertex. For the 
reconstruction of the undetected neutrino, the missing energy and momentum in the 
cm. frame are defined as -E m i ss = 2E\ )Cam — EiEi and p m i ss = —Uipl, respectively, 
where -Ebeam is the beam energy in the cm. frame, and the sums include all charged 
and neutral particle candidates in the event. We require E'miss > GeV, and the 
neutrino 4-momentum is taken to be p v = ( |p m iss | , Pmiss), since the determination 
of i?miss is more accurate than that of the missing energy. As in the analysis of 
B° —7- D*~£ + is, using the variable cosOby (Eq. 4-7), the B momentum vector is 
constrained to lie on a cone centered on the tt~£ + direction; signals can then be 
selected by requiring |cos#by| < 1. Background from continuum e + e~ — > qq(q = 
u, d, s, c) jets are reduced using an event topology requirement based on the second 
Fox-Wolfram moment. Signals are extracted, for each of 13 q 2 bins ranging from 
to 26 GeV 2 /c 2 , by fitting the two-dime msional distribution of the beam energy 

constrained mass Mt, c = \J ^ cam — + Pi + P~l\ 2 and the energy difference AE = 

-E-beam — (E w + Eg + E v ). Figure 19 shows the obtained q 2 distribution. The total 
branching fraction, integrated over the entire q 2 region, is 

B(B° -»• 7r~e + v) = [1.49 ± 0.04(stat) ± 0.07(syst)] x 10~ 4 . (4-20) 

The value of \V u b\ can be determined from the measured differential q 2 dis- 
tribution using Eq.(4T9). Following the procedure proposed by the FNAL/MILC 
collaboration , 64 ) \V u b\ can be extracted from a simultaneous fit to experimental and 
lattice QCD results from the FNAL/MILC collaboration, as shown in Fig. 20. In 
this approach, q 2 is transformed to a dimensionless quantity z, and both the exper- 
imental and lattice QCD distributions are fit to a third-order polynomial with \ V u b\ 
determined as a relative normalization between the lattice QCD and experimental 
results. We find \V ub \ = (3.43 ±0.33) x 10~ 3 , as shown in Table VIII. The table also 
lists | V u b | values determined using only a fraction of the overall phase space, leading 
to less precise but statistically compatible results. The form factor f + (q 2 ) predic- 
tions are based on the light cone sum rule (LCSR) and lattice QCD (LQCD), which 
can be applied in the regions q 2 < 16 GeV 2 /c 4 and q 2 > 16 GeV 2 /c 4 , respectively. 

Table VIII. Summary of \V u b\ results from a recent B° —5- n~ ' tv measurement by Belle. 



Theory 


q< (GeVVc 4 ) 


\V ub \ (xlCT 3 ) 


LCSR 65) 


< 16 


3.64±0.1lt£!o 


HPQCD 66) 
FNAL 50) 


> 16 

> 16 


3.55 ±0.13±g;|? 
3.78 ±0.141°;^ 


FNAL/MILC 64) 


all regions 


3.43 ± 0.33 



4.3.2. \V u b\ from inclusive B — > X u lv decays 

For inclusive B — > X u lv decays, the theoretical description relies on the OPE, 
as in the case of inclusive B — > X c lv decays. However, B — > X u lu decays are 
about 50 times less abundant than B — > X c tv decays, and thus the experimental 
sensitivity to B — > X u tv and \ V U \\ is highest in the region of phase space that is less 
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Fig. 19. Measured q 2 distribution for the 
B° — > it - Iv decay. The curve repre- 
sents a fit to an empirical form factor pa- 
rameterization. The four histograms show 
various form factor predictions (dashed: 
ISGW2; plain: HPQCD; dotted: FNAL; 
dot-dashed: LCSR). 



Fig. 20. \V u b\ extraction from a simultane- 
ous fit to experimental (closed circles) and 
FNAL/MILC lattice QCD results (open 
circles) . 



impacted by the dominant background from B — > X c lv decays. In this phase space 
region, however, non-perturbative corrections are kinematically enhanced, and as a 
result non-perturbative dynamics become an 0(1) effect. Extracting \V u b\ requires 
the use of theoretical parameterizations called shape functions (SF) to describe the 
unmeasured regions of phase space. 

A classical method is to measure the lepton momentum spectrum at the end- 
point of the spectrum (pf 71 > 2.3GeV/c), where the b — > c decay is forbidden. This 
method allows the measurement of \V U (,\ with small data samples, but suffers from a 
large extrapolation error, because only a limited portion of the phase space (~10 % 
of the total) is measured. Belle reported a result using this method in 2005. 67 ) The 
high luminosity data at Belle enable us to also measure kinematic variables such as 
the invariant mass of the X u hadronic system, mj, and the four-momentum transfer 
of the B meson to the X u system, q. This enables us to control the experimental and 
theoretical errors by optimizing the region of phase space for the measurement. Belle 
reported the first measurement using m x — q 2 for the inclusive B — > X u lv decay. 68 ) 

More recently, Belle reported a measurement of the partial branching fraction 
of B — > X u iu decays with a lepton momentum threshold of 1 GeV/c using a multi- 
variate data mining technique, with a data sample containing 657 x 10 e BB pairs. 69 ** 
This method allows us to access ~ 90% of the B — > X u lv phase space and minimizes 
the dependence on an SF. The measurement is made by fully reconstructing one 
B meson (Stag) in hadronic decays, and measuring the semileptonic decay of the 
other B meson (-B s ig) with a high momentum electron or muon. The B — > X u lv 
decays are selected based on a nonlinear multivariate boosted decision tree (BDT), 
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which incorporates a total of 17 discriminating variables, such as the kinematical 
quantities of candidate semileptonic decays, number of kaons in the event, M\, c of 
-Stag) etc. The candidates passing the selection of the BDT classifier are analyzed 
in a two-dimensional fit in the (mx, Q 2 ) plane. The hadronic invariant mass mx 
is calculated from the measured momenta of all charged tracks and neutral clusters 
that are not associated to B tSLg reconstruction or used as a lepton candidate. The 
momentum transfer is calculated as q = Prus) ~ PB ts , g — Px- Figure 21 shows the 
one-dimensional projections of the (mx,q 2 ) distribution with a lepton momentum 
requirement of p^ B > l.OGeV/c, fitted with distributions for the B — > X u £u signal, 
B — > X c lv and other backgrounds mainly from secondary and misidentified leptons. 
The partial branching fraction for p* B > l.OGeV/c is 

AB{B -> XJv\ p*f > l.OGeV/c) = 1.963 x (1 ± 0.088(stat) ± 0.081(syst)) x 1(T 3 . 

(4-21) 

A \ V u b\ value is obtained from the partial branching fraction using \ V u b\ 2 = AB u t v j (tbAR), 
where AR is the predicted B —> X u tv partial rate in the given phase space region, 
and tb is the average B lifetime. Table IX presents \V u b\ results based on differ- 
ent theoretical prescriptions that predict AR. Here the results were obtained by 
the Heavy Flavor Averaging Group (HFAG) using the most recent calculations and 
input parameters. 28 ^ The results are consistent within their stated theoretical un- 
certainties, and have an overall uncertainty of ~ 7%. 

As described above, there is a tension between the \V U \,\ values extracted from 
the exclusive and inclusive methods, which are subject to further clarification with 
improved experimental and theoretical errors in the future. 




M x (GeV/c 2 ) q 2 (GeV 2 /c 2 ) 

Fig. 21. Projections of the mx — q 2 fit in bins of mx (left) and q 2 (right). 

4.4. Purely leptonic B~ — > t~~vi [t = e, fi, or r) decays 

In the SM, B~ — > l~v^ decays to purely leptonic final states {I = e, fi or r) 
occur via annihilation of the two quarks in the initial state, b and u, to a W~ boson 
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Table IX. \V u b\ values obtained using the inclusive B —¥ X u £nu measurement by Belle and input 
parameters (mi, and /xj). The errors quoted on \ V u b\ correspond to experimental and theoretical 
uncertainties, respectively. 



Theory 


m b (GeV) 


l& (GcV x ) 


\V ub \ (xlO 


-*) 


BLNP 70) 


4.588 ± 0.025 


0.189;™ 


4.47 ±0.27^ 


-0.19 
-0.21 


DGE 71) 


4.194 ±0.043 




4.60 ± 0.27j 


-0.11 
-0.13 


GGOU 72 > 


4.560 ± 0.023 


0.453 ± 0.036 


4.54 ± 0.27j 


-0.10 
-0.11 


ADFR 73) 


4.194 ±0.043 




4.48 ± 0.30j 


-0.19 
-0.19 



(Fig. 22). The branching fraction for a B — > t decay is given by 

B(B - -> n t ) = ( 1 - fi\V^r B , (4-22) 

where Gf is the weak interaction coupling constant, and tub are the lepton and 
B + meson masses, respectively, tb is the B~ lifetime, \V u j,\ is the magnitude of a 
CKM matrix element, and fs is the B~ meson decay constant. All these input 
parameters have been directly measured with good precision except for fg. The 
value of fs can be obtained using LQCD calculations. Since LQCD calculations are 
based on first principles of QCD, it is possible to calculate the SM expectation for 
B(B~ —7- l~i>i) with high precision. Therefore, measurement of fs via B~ — > l~i>i 
decays can provide a stringent test of the LQCD, within the framework of the SM . 

On the other hand, particles from 
physics beyond the SM, for example, a 
charged Higgs boson in supersymmetry or a 
generic two-Higgs doublet model, may take 
the place of the W~ in Fig. 22 and modify 
the branching fraction. Moreover, in the 
minimum flavor violation NP scheme, it is 
expected that the relative branching frac- Fig. 22. A Feynman diagram for the SM 
tions of charged lepton modes will remain B~ — > t~v t process, 

the same as those predicted by the SM. Ac- 
cordingly, measuring the branching fractions of B~ — > l~vn {I = e, ji or r) modes 
and their relative ratios can provide a very sensitive probe for NP beyond the SM. 

Due to helicity suppression, the branching fraction (Eq. 4-22) is proportional to 
the square of the charged lepton mass, mf . As a result, the SM branching fractions 
for e~v e and modes are suppressed in comparison to the t~v t mode by factors 

of ~ 10 7 and ~ 200, respectively. At the time of this report, there exists evidence 
for B~ — > t~u t from Belle 74 )' 75 ) and BaBar, 76 ) but no evidence has yet been found 
for the B~ — > e~9 e and B~ — > modes. 

4.4.1. B- -)• t~v t 

While the large mass of the r lepton significantly enhances the branching fraction 
of B~ —7- t~v t compared to other modes, the presence of one or more neutrinos 
from the r decay make it difficult to cleanly detect B~ — > t~v t decays. In the 
process e + e~ — > T(4S) — > BB, signal sensitivity is greatly improved by completely 
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reconstructing or "tagging" one B meson (-B tag ); the signature of the signal is then 
searched for in the other B meson (B S ig)- Experimentally, two different tagging 
methods have been applied to measure B(B~ — > t~i> t ): reconstructing a full decay 
chain of a hadronic final state ("hadronic tagging") or reconstructing all particles 
except for a neutrino in semileptonic -Btag ~~ * D^tv decays ( "semileptonic tagging"). 

4.4.2. Hadronic tagging analysis 

The first evidence for B~ — > t~v t decays was obtained in a hadronic tagging 
analysis by Belle 74 ) using 449 x 10 6 BB events, which obtained B{B~ — > t~v t ) = 
(1.79 49-0 5i) x 10~ 4 - Recently, Belle has updated the hadronic tagging analysis of 
B~ —7- t~u t , analyzing the full Belle data sample containing 772 x 10 6 BB events. 77 -* 

In the most recent analysis, the data sample is fully reprocessed with much 
improved tracking and slightly improved neutral cluster detection. A new hadronic 
tagging algorithm using a Bayesian artificial neural network has been developed and 
applied to the analysis. 78 - > As a result of all these improvements, the statistics of 
the -Btag sample has increased by nearly a factor of three. Figure 23 shows the Mb c 
distribution of -B tag candidate events, in comparison with that from the previous 
analysis. 



x10 3 




Fig. 23. Mbc distributions for the -Btag candi- 
date events. The triangles, open circles, 
and solid circles represent the distributions 
obtained by applying the original tagging 
algorithm 74 - 1 to the previous data set, ap- 
plying improved hadronic tagging to the 
previous data set, and applying improved 
tagging to the latest fully reprocessed data 
set, respectively. The solid and dotted 
curves show the sum and the background 
component, respectively, of the fit to the 
full data sample. 




ML s (GeV 2 /c 4 ) 



Fig. 24. Distributions of -Becl (top) and 
M^ iss (bottom) combined for all the r~ de- 
cays. The Mmi ss distribution is shown for 
the signal region -Eecl < 0.2 GeV. The 
solid circles with error bars are data. The 
solid histograms show the projections of the 
fits. The dashed and dotted histograms 
show the signal and background compo- 
nents, respectively. 



Once the -B t ag candidates are selected, we search for B — > r v T decays using 
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the particles not belonging to 2? tag in these events. The r~ lepton is identified in four 
decay modes: r~ — > e~i> e v T , p~v^v T , tt~v t , and ir~ir°v T . Signal candidate events 
are required to have only one track with charge opposite to B tag . The charged tracks 
are required to be consistent with being either an electron, muon, or pion. For the 
r~ — > -k~it®v t mode, with 7r° — > 77, the invariant mass of the ir~ir° system must 
be within 0.15 GeV/c 2 of the nominal p~ mass. There should be no other detected 
particles after removing the particles from the -Btag and the charged tracks and 7r°'s 
from the B slg . In particular, events containing extra tt° and K\ candidates are 
rejected. 

The signal yield is evaluated by fitting the two-dimensional distribution of Eecl 
and M^ iss , where -Eecl is the sum of the energies of neutral clusters that are not 
associated with either the Bt ag or the 7r° candidate in the r~ — > -k~it®v t decay 
and M^ iss is the missing mass squared defined by M^ iss = (Eqm - Es tag — -EB sig ) 2 - 
l^-Btag +P-B S i g | 2 w ith the energies and momenta measured in the CM frame. To reduce 
background, we require M^ iss > 0.7 (GeV/c 2 ) 2 . Figure 24 shows the projections of 
the result of the fit on Eecl and M^ iss where the four r decay modes are combined. 
The preliminary fitted signal yield is 62 ^22 ^ 6 events and the branching fraction 
Shad (in the hadronic tagging analysis) is 

B had (B- -> t-9 t ) = (0.72 t° f 5 ± 0.11) x 10~ 4 . (4-23) 

The signal significance is 3. Oct including systematic uncertainty. This result is con- 
sistent with the previous measurement considering the overlap of the event samples. 

4.4.3. Semileptonic tagging analysis 

In the semileptonic tagging analysis of B~ — > T~i> T , B tSLg is reconstructed in 
B + — > D*°£ + v and B + — > D £ + u decays, where £ is an electron or muon. Since 
semileptonic tagging imposes fewer constraints on the B slg kinematics, only r~ de- 
cays to l~V(V T {£ = e,fi) and ir~v T are used for B S i g reconstruction. Except for D 
or D , one £ + for B tag , and one £~ or 7r~ for B s [ g , we allow no other charged track 
or neutral particle in the event. 

One of the main variables to suppress background events is the cosine of the 
angle, cos9 B £>(*)£, between the momentum of B tag and that of D and £ + system. 
This variable is defined in the same way as the variable cos 9 by discussed in the 
previous section, but with Y = D^°£ + . Correctly reconstructed B tSug candidates 
populate the physical range — 1 < cos6> BD (»)^ < 1. Signal candidates are selected 
based on P^ m (the lepton momentum of B tSLg in the CM frame), cos B D (*)^, and P s c j™ 
(the CM- frame momentum of the charged track from B s [ g ). The selection criteria 
depend on the r decay mode of B s { g . After all selections, the signal yield (n s ) is 
obtained by fitting the -Eecl distribution. From a combined fit to the three r~ 
decay modes, n s = 143 +35 events is obtained. The signal significance is found to 
be 3.6<t including the systematic uncertainty. The branching fraction £>sl (in the 
semileptonic tagging analysis) is 

bsUb- t-u t ) = (1.54 ±8:H±g:S) >< ^ ■ (±24) 
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4.4.4. The combined result 

The two results, Shad and Ssl, are combined after taking the correlation in 
the systematic uncertainties between the two results into account. 79 ) The signal 
significance for the combined result is 4. Oct and the average branching fraction is 

B(B~ -»• t~v t ) = (0.96 ± 0.22 ± 0.13) x 10~ 4 . (4-25) 

The result is consistent with the SM expectation obtained from other experimental 
constraints. Using this result along with the input values found from the most recent 
world averages, 15 ) we obtain fB\V ub \ = (7.4 ± 0.8 ± 0.5) x 10~ 4 GeV. This result sets 
stringent constraints on the parameters of various models involving charged Higgs 
bosons. 

4.4.5. B -»• tTv t {I = e,n) 

As discussed above, the B~~ — > e~v e and B~ — > decays are suppressed 

compared to B~ — > t~v t due to helicity suppression. On the other hand, these 
decays have a clear experimental signature: the monochromatic energy of the charged 
lepton in the rest frame of the signal B. Two methods have been applied to measure 
these decays: a loose reconstruction analysis and a hadronic tagging analysis. 

In the loose reconstruction analysis, where a data sample containing 277 x 10 6 
BB pairs is used, 80 ) the signal candidates are selected mainly via a tight requirement 
on pf, which is the charged lepton momentum (magnitude) in the signal B rest 
frame. The signal yield is then obtained by fitting the Mt, c distribution, where M^ c 
is calculated by including all detected particles in the event except for the signal 
charged lepton. No significant excess of signal in any mode is found. We set the 
following upper limits on the corresponding branching fractions at the 90% C.L.: 

B(B~ -»• e~v e ) < 0.98 x 10~ 6 , (4-26) 
B(B- -»• n'v^) < 1.7 x 10~ 6 . (4-27) 

The hadronic tagging analysis is based on a method similar to that described in 
Sect. 4.4.2 and uses the full data set of Belle containing 772 x 10 6 BB pairs. After 
selecting signal candidates primarily using the Mt, c and AE variables of the Stag and 
requiring that the B s i g be consistent with B~ — > £~i>e, including a requirement on 
Eecl, the expected background in the signal region, 2.6 < pf < 2.7 GeV/c, is much 
less than one event. The background estimate is determined by examining data and 
MC events in the sideband of pf below the signal region. 

The signal yield is obtained by counting the events in the pf signal region. No 
events are found in any mode and we set 90% C.L. upper limits on the branching 
fractions using the POLE 81 ) program taking the uncertainty in signal efficiency and 
the expected background with its uncertainty into account. The preliminary upper 
limits (at 90% C.L.) for the branching fractions Shad (by hadronic tagging analysis) 
are: 82 ) 



ShadCB~ -> e ~^e) < 3.5 x 10~ 6 , 
Shad (B~ -> < 2.5 x 10- 6 . 



(4-28) 
(4-29) 
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Although the constraints are not as stringent as those obtained in the loose recon- 
struction analysis, the amount of background is much smaller, nearly zero; hence it 
is anticipated that the sensitivity may improve almost linearly with the increase of 
statistics. Therefore, the hadronic tagging analysis will be very interesting in the 
next-generation super U-factory experiments such as Belle II. 

4.5. B -)• D^tv decays 

Compared to ordinary semileptonic decays B° — > D* + l~V£ with £ = e or [i, 
B — > D^t~v t decays, occurring through a quark-level b — > ct~v t process, are 
suppressed because of the large r mass. The predicted branching fractions, based on 
the SM, are approximately 1.4% and 0.7% for B — > D*t~v t and B — > Dt~v t decays, 
respectively. 83 ^ On the other hand, the large r lepton mass makes them sensitive to 
interactions with a charged Higgs, where the H + may replace the virtual W, thereby 
modifying the branching fraction. Therefore, these B — > D^t~v t modes can be a 
very effective probe to search for indirect evidence of charged Higgs or other NP 
hypotheses beyond the SM. Moreover, compared with B~ — > t~9 t , these decay 
modes provide more observables to search for NP, e.g. the polarization of the r 
lepton. On the experimental side, however, it is very difficult to measure these 
modes because of the multiple neutrinos in the final state, the low lepton momenta, 
and the large associated background contamination. 

The first observation of B — > D^*\~v T decays was reported by Belle in the 
B — > D* + t~v t mode using an event sample of 535 x 10 6 BB pairs. 84 ) In contrast 
to the hadronic tagging analysis (see Sect. 4.4.2), a loose reconstruction of the ac- 
companying B (Stag), where all particles not belonging to the signal decay chain are 
included without taking subdecay information into account, was used and tighter 
kinematic constraints were applied for improved background suppression. The sig- 
nal yield was obtained by fitting the distribution of the beam-constrained mass Mb c 
of the Stag- A clear signal excess of 60 *ii events was observed with a significance 
of 5.2a including systematic uncertainties. The measured branching fraction was 
B(B° -> D* + t-u t ) = (2.02 1°; 4 ° ± 0.37)%. 

Belle has also published measurements of other B — > D^t~v t decay modes. 
Analyzing a data sample of 657 x 10 6 BB pairs, using a similar analysis to that 
described above, 446 +|g events of the B~ — > D*®t~v t decay mode are observed 
with a significance of 8.1<r and 146 ^ events of the B~ — > D t~v t decay mode 
are obtained, providing the first evidence of this mode with a significance of 3.5<r. 85 ) 
The branching fractions are B{B~ -»• D*°r"^ T ) = (2.12 ±^±0.29)% and B(B~ -> 
D°t-u t ) = (0.77 ± 0.22 ± 0.12)%. 

A preliminary branching fraction of the B — > D + t~v t mode is measured by 
an analysis that uses a hadronic tagging method similar to the one described in 
Sect. 4.4.2: B(B° -»• D + t~v t ) = (1.01 l^lan ± 0.10)%, 86 ) where the third error 
comes from the branching fraction uncertainty of the normalization mode, B° — > 
D + l~V£. The branching fractions of the other B — > D^t~v t decay modes are also 
obtained in this analysis; the results are consistent with published results. 84 )' 85 ) 

Recently, BaBar has claimed that the branching fractions of B — > D*t~v t and 
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B — > Dt~v t are larger than SM expectations at a combined significance of 3.4<r. 87 ^ 
We note that all the branching fractions of B — > D^t~v t modes measured by 
Belle are also larger than the SM-predicted values. 83 '' It will be interesting to see the 
final Belle results on these modes using improved hadronic tagging and the full data 
sample of 772 x 10 6 B~B pairs. 

§5. Rare B decays 

5.1. Charmless hadronic decays 




< 4 + 



Fig. 25. Feynman diagrams for the quark-level transitions that dominantly contribute to charmless 
hadronic B decays: (left) color allowed and (middle) color suppressed b — > u tree diagrams, and 
(right) b — > (s, d)g penguin diagrams. 

Charmless hadronic B decays give rise to final states with two or more hadrons 
that do not contain any charm quark. These decays are suppressed in the SM, mostly 
proceeding via the CKM-suppressed b — > u tree level transition and b — > (s, d)g 
penguin diagrams, as shown in Fig. 25. Compared to the CKM-favored b — > c 
transition such as B° — > J/ip K°, the golden channel for determining the angle (f>i of 
the unitarity triangle (Sect. 3), their branching fractions are about two to four orders 
of magnitude lower. By virtue of this suppression, charmless decays provide a good 
window to probe new physics beyond the SM. For instance, these rare decays have 
the potential to reveal the contribution of heavy, non-SM virtual particles in penguin 
loops. Branching fraction calculations within the SM — whether they are based on 
QCD factorization, 88 ^ SU(3) flavor symmetry, 89 ) or perturbative QCD 90 ) — suffer 
from large theoretical uncertainties. However, one can examine physics observables in 
which theory errors as well as common experimental systematic uncertainties largely 
cancel out. Such observables include direct CP asymmetries, ratios of branching 
fractions, and longitudinal polarization fractions (in the case in which the decay final 
states consist of two vector particles). In this section, we summarize Belle's results 
on charmless B decays, which have resulted in close to 100 journal publications. 

5.1.1. Experimental methodology 

Before examining various categories of charmless decays, we wish to describe the 
important experimental considerations. As these decays are suppressed in the SM 
one needs to be extremely careful in devising selection algorithms to select candidate 
events and fitting methods used for extracting the final signal yields. We identify B 
mesons using two kinematical variables: the beam-energy constrained mass, M^ c = 
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V^beam ~ \pb\ 2 , and the energy difference AE = E B - E h 

earn) where -Ebcam is the 

beam energy, and Eb and ps are the energy and momentum of B candidates in the 
CM frame, respectively. The dominant background contribution is from e + e~ — > qq 
(q = u, d, s, c) continuum processes. To suppress this background, we use variables 
based on the event topology; these selections rely on the fact that B decays are nearly 
isotropic, in contrast to jet-like continuum events. In some analyses, additional 
discrimination is provided by variables pertaining to the nature of B decay, e.g., 
the flight-length difference along the beam direction between the signal and recoil 
B decay. All this information is combined into either a likelihood ratio or a neural 
network to optimize the sensitivity. B decays proceeding via a CKM-favored b — > 
c transition can have a final state that is either the same as our signal or mis- 
reconstructed. Since branching fractions for b — > c decays are much larger and the 
charm mesons involved are quite narrow, we suppress their contributions by applying 
a veto on the reconstructed invariant mass of daughter particles of the charm meson. 
Backgrounds from other B decays, especially those due to particle misidentification, 
pose a special challenge. The AE and charged-hadron identification variables help 
in discriminating such backgrounds. The final signal yield is extracted by means of 
an unbinned maximum likelihood fit to the discriminating variables, Mb c , AE, and 
continuum suppression variable (in some analyses we have used an optimized and 
tight requirement on the latter using the expected signal significance from Monte 
Carlo simulations as a figure of merit). For decays involving narrow or non-zero 
spin resonances in the final state, we employ the invariant mass and helicity angle 
distributions to further enhance the sensitivity of our results. 

5.1.2. Results on charmless decays 

Charmless hadronic B decays can be roughly divided into "two-body", "quasi- 
two-body (Q2B)", and "three-body" categories. Although two-body decays are easily 
identifiable by the presence of two long-lived final state particles, such as Kit, the 
latter two classes of decays are somewhat intertwined. For instance, when one per- 
forms the Dalitz plot analysis of a three-body final state, one can access information 
on related Q2B decays along with the three-body nonresonant decay. The Dalitz 
plot approach is the most appropriate method when dealing with broad intermedi- 
ate resonances, e.g., p(770). However, if the intermediate resonances are narrow or 
the resonances decaying to the same final state do not interfere, we can use a Q2B 
approach where the interference is accounted for as an additional source of systematic 
error. Since the aforementioned distinction between the two categories is not "black 
and white", we begin the discussion with two-body decays, then Q2B decays are 
described together with intermediate resonance final states from three-body Dalitz 
analyses, and finally we move on to results on three-body nonresonant and inclusive 
final states. 

Two-body decays 

B meson decays to two stable hadrons are kinematically easy to identify, having 
average momenta larger than a typical B decay. The kinematics also provides a 
good handle on the continuum background. A formidable challenge is posed by 
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the feed-across background arising due to particle misidentification (mostly, kaons 
misidentifled as pions). We tackle this issue by performing a simultaneous fit to the 
event samples that can cross feed into each other. 

Table X summarizes the branching fraction and CP asymmetry results for var- 
ious charmless two-body B decays from Belle. The most notable result here has 
been the observation of a non-zero difference of CP violation asymmetry in the 
B° -»• K+-K- and B + -»• K + ir° decays: AA K * = A C p{K + tt°) - A C p(K + ir-) = 
+0.112 ± 0.027 ± 0.007. 91) This discrepancy, also called the AA Kn puzzle, may be 
explained either by a large contribution from the color-suppressed tree diagram 92 ) 
or a new physics contribution in the electroweak penguin. 93 ) Before concluding on 
this issue, we must improve the uncertainties on CP violation results for the de- 
cay B° — > K°7T . This would allow us to precisely test the prediction of an isospin 
sum rule 94 ) given by A CP (K+ir-) + A CP {K^+)0£± - A CP {K+^) 2 ^1] - 

Acp(K°tt°) jy/Pf^f] = 0> wn ere r is the partial width. Belle's latest update 91 ) 
reports a sum of —0.270 ± 0.132 ± 0.060 with 1.9<r significance. 



Table X. Data samples used (N B g), branching fractions (B), 90% confidence-level upper limits on 
B (UL), and direct CP asymmetries (Acp) obtained for various charmless two-body B decays. 
The two uncertainties quoted here and elsewhere are statistical and systematic, respectively. 



Final state 


N BE (10 6 ) 


B (KT 6 ) 


UL (10^ 6 ) 


Acp 


Ref. 


K+K~ 


772 


0.10 ±0.08 ±0.04 


0.20 




91) 


K+K° 


772 


1.11 ±0.19 ±0.05 




±0.014 ±0.168 ±0.002 


91) 


K°K° 


772 


1.26 ±0.19 ±0.05 






91) 


K+TY- 


772 


20.00 ± 0.34 ± 0.60 




-0.069 ±0.014 ±0.007 


91) 


K+tv" 


772 


12.62 ±0.31 ±0.56 




±0.043 ± 0.024 ± 0.002 


91) 


K°7T + 


772 


23.97 ±0.53 ±0.71 




-0.011 ±0.021 ±0.006 


91) 


K°7T° 


772 


9.68 ± 0.46 ± 0.50 






91) 




772 


5.04 ±0.21 ±0.18 






91) 


7T+7T 


772 


5.86 ±0.26 ±0.38 




±0.025 ± 0.043 ± 0.007 


91) 


ir°ir° 


275 


O+0.4+0.2 
z -^-0. 5-0.3 




+0A4+° H ± 0.17 


95) 


PP 


449 




0.41 




96) 


pA 


449 




0.49 




96) 


AA 


449 




0.69 




96) 



Quasi-two-body (Q2B) decays 

Q2B analyses assume that the intermediate resonances decaying to the same 
final state (such as p(770) and /o(980) decaying to ir + ir~) do not interfere. This 
treatment allows us to compare branching fraction results with measurements from 
earlier experiments, in which the effects of interference were treated as a part of the 
systematic error. The extent and nature of the background (as most charmless B 
decays suffer from a low signal-to-background ratio) and of the nonresonant signal 
component strongly influence our analysis strategy. The helicity angle plays an 
important role in such Q2B analyses when the intermediate resonances have a non- 
zero spin; we can use it either as a simple selection criterion or to extract physics 
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observables, such as the fraction of longitudinal polarization fi, directly from the fit. 
The helicity angle Oh for a resonance is defined as the angle between the momentum 
vector of one of its daughter particles and the direction opposite to the 5-meson 
momentum in the resonance rest frame. 97 ) 

In the discussions that follow, the decays have been grouped according to their 
spin. For each spin grouping, the results for branching fractions, direct CP asym- 
metries, and longitudinal polarization fractions (where applicable) are listed in the 
accompanying tables. In Table XI we start with final states comprising at least an 
■n, or rj meson that is reconstructed in the two channels r\ — > 77 and rj — > ir + ir~ir°, 
or rj' — > r]iT + iT~ and 7/ — > p°7, respectively. Among the highlighted results are the 
first observation of B° — > r/K°, and evidence for direct CP violation in the decays 
B + —7- T]K + and B + — > r/7r + with significances of 3.8cr and 3.0a, respectively. 98 - ) The 
latter results call for a large interference between the b — > s penguin process and the 
CKM-suppressed, color-favored b — >■ u tree transition, both of which contribute to 
B+ ^ V h+ (h = K,ir). 



Table XI. Data samples used (N B g), branching fractions (B), 90% confidence-level upper limits on 
B (UL), and direct CP asymmetries {Acp) obtained for various charmless Q2B decays with an 
r) or rj' meson in the final state. 



Final state 


N BS (10 6 ) 


B (ICT 6 ) 


UL (10" 6 ) 


A C p 


Ref. 


7]K+ 


772 


2.12 ±0.23 ±0.11 




-0.38 ±0.11 ±0.01 


98) 


V K° 


772 


1.27±g;|| ± 0.08 






98) 


r]iT + 


772 


4.07 ±0.26 ±0.21 




-0.19 ±0.06 ±0.01 


98) 




152 


1.2 ±0.7 ±0.1 


2.5 




99) 


m 


152 


0.7i°;I±0.1 


2.0 




99) 


V 'K+ 


386 


69.2 ±2.2 ±3.7 




±0.028 ± 0.028 ±0.021 


100) 


T)'K° 


386 


58.9t^'j ± 4.3 






100) 


rj'iT + 


386 


1.76 ±0.671^50.14 




±0.2015',;^ ± 0.04 


100) 


r,V 


386 


2.79 ± 1.021^5-0.34 






100) 


v'v 


535 




4.5 




101) 


vW 


535 




6.5 




101) 


V K*+ 


449 


19.3li g ± 15 




±0.03 ±0.10 ±0.01 


102) 


r,K*° 


449 


15.2 ±1.2 ±1.0 




±0.17 ±0.08 ±0.01 


102) 


VP + 


449 


4.1±^;| ±0.4 


6.5 


-0.041° ± 0.01 


102) 


VP° 


449 


0.84l°;^ ±0.19 


1.9 




102) 


V 'K*+ 


535 




2.9 




101) 


v 'k*° 


535 




2.6 




101) 


v'p + 


535 




5.8 




101) 


v'p° 


535 




1.3 




101) 




535 




2.2 




101) 


n'4> 


535 




0.5 




101) 



The branching fractions and CP asymmetries for other Q2B decays without an 
rj or r\ meson in the final state are summarized in Table XII. Most of the results 
are obtained as a by-product of a three-body Dalitz plot analysis. The systematic 
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uncertainties in the table include the experimental systematic as well as Dalitz-plot 
model dependence, where applicable. We report the first evidence of CP violation 
in the Q2B decay B + — > p°K + exceeding the 3a level. Note that this was the first 
evidence for direct CP violation in a charged meson decay, a phenomenon that was 
already observed in decays of neutral K 112 ^ and 5 91 )> 113 ) mesons, and very recently 
in D° decays. 114 ) 

In Table XIII we present results obtained from the vector-vector final states. 
One naively expects B — > VV decays to be dominated by longitudinal polarization 
amplitudes since /l = 1 — 4my /rag ~ 0.9, 121 ) where my (m^) is the mass of the 
vector (B) meson. Contrary to this expectation, it is found out that the decays 
dominated by the b — > s penguin transition such as B — > <pK* have fi values closer 
to 0.5. However, decays proceeding via the b — > u tree diagram, notably B — > pp, 
follow the expected trend. This so-called polarization puzzle could be explained by 
the presence of new particles in the penguin loop. 122 ) However, large SM corrections 
appear to be a more plausible explanation. 



Table XII. Data samples used (N B g), branching fractions (B), 90% confidence-level upper limits on 
B (UL), and direct CP asymmetries (Acp) obtained for various charmless Q2B decays without 
an rj or rj meson in the final state. 



Final state 


N BS (10 6 ) 


B (10" 6 ) 


UL (10~ 6 ) 


A C p 


Ref. 


/o (980) if 


388 


7.6 ± l-7±?:? 






103) 


/o(980)AT+ 


386 


8.78 ± 0.82±°;®| 




-0.077 ± 0.0651°;!$; 


104) 


/ 2 (1270).K" + 


386 


f.33±0.30l°;3i 




-0.59 ±0.22 ±0.04 


104) 


/ 2 (1270)A: 


388 




2.5 




103) 


7^ (f430) +7T- 


388 


49.7±3.8l|-J 






103) 


J ftT *(1430)°7r + 


386 


51.6±1.7±™ 




±0.076 ± 0.0381^22 


104) 


K-+-K- 


388 


8.4 ± l.llo.^ 




-0.21 ±0.11 ±0.07 


103) 


K*°it+ 


386 


9.67 ± 0.64tg-jg 




-0.149 ±0.064 ±0.022 


104) 


K"%° 


85 


0.4l^ ±0.1 


3.5 




105) 


coK + 


388 


8.1 ±0.6 ±0.6 




+0.05±^^ ± 0.01 


106) 


ujK° 


388 


4.4l°;® ± 0.4 






106) 


LJTT + 


388 


6.9 ±0.6 ±0.5 




-0.02 ±0.09 ±0.01 


106) 


U,7f° 


388 


0.5t°-j±0.1 


2.0 




106) 


(j>K+ 


f52 


9.60 ± 0.92lJ;°| 




±0.01 ±0.12 ±0.05 


107) 


<t>K° 


85 


9.0t?'j ± 0.7 






108) 


<j>1T + 


657 


n no+O.09+0.06 
U.UO_g .08-0.03 


0.33 




109) 




657 


n n7 +0.06+0.04 
u - u ' -0.04-0.08 


0.15 




109) 


0(1680)^+ 


f52 




0.8 




107) 


p~K+ 


85 


ic 1+3.4+2.4 
1J - 1 -3.3-2.6 




, n 99+O.22+O.O6 
-|-U.ZZ_ 23-0 .02 


105) 


p°K+ 


386 


3.89 ± 0.47±21? 




±0.30 ± o.n±g:SJ 


104) 


p a K a 


388 


6.1 ± 1.0±i2 






103) 


p+7V° 


f52 


13.2±2.3ti;| 




±0.06 ±0.171^ 


110) 


P c v+ 


32 


s.ol^o ±0.7 






111) 


A° 


449 


3.0 ±0.5 ±0.7 






32) 


p^7V ± 


449 


22.6 ±1.1 ±4.4 






32) 



Physics achievements from the Belle experiment 



43 



Table XIII. Data samples used (N BB -), branching fractions (B), 90% confidence-level upper limits 
on B (UL), longitudinal polarization fraction (/l), and direct CP asymmetries (Acp) obtained 
for various charmless Q2B decays with vector-vector final states. 



Final state 


N BB (10 6 ) 


b (itr 6 ) 


UL (10~ 6 ) 


h 


A C p 


Ref. 




657 




0.2 






115) 


K*°K*° 


657 


n Ofi+0.33+0.10 
u -^°-0.29-0.08 


0.8 






115) 


K*°p+ 


275 


8.9 ±1.7 ±1.2 




0.43 ±0.1lt°;°l 




116) 


K*°p° 


657 


r, -i +0.8+0.9 
z - 1 -0.7-0.5 


3.4 






117) 


ujK*° 


657 


1.8 ± 0.7±g;| 




0.56 ± 0.29^ 08 




118) 


<j>K* + 


85, 257 


fi 7 +2.1+0.7 
u - ' -1.9-1.0 




0.52 ±0.08 ±0.03 


-0.02 ±0.14 ±0.03 


108), 119) 


4>K*° 


85, 257 


10 n+ 16 + - 7 

lu - u -1.5-0.8 




0.45 ± 0.05 ± 0.02 


±0.02 ±0.09 ±0.02 


108), 119) 


p+p- 


275 


22.8 ± 3.8± 2 2 3 6 




0-94lo;o4 ± 0.03 




34) 


p+p" 


85 


31.7 ± 7.1±|;f 




0.95 ±0.11 ±0.02 


±0.00 ±0.22 ±0.03 


33) 


p°p° 


657 


0.4 ± 0.4±'±J 


1.0 






120) 



Three-body decays 

As was discussed above, the Dalitz-plot method is the most robust analysis 
technique for a three-body decay, especially for B — > 3 pseudoscalars. This method 
has greater complexity but at the same time provides a better understanding of the 
underlying physics. A subtle point here is that one needs a good deal of statistics 
before carrying out a full-fledged Dalitz plot analysis. As the integrated luminosity 
continued to increase at Belle, starting with measurements of inclusive branching 
fractions and charge asymmetries, various rare decay analyses slowly evolved into 
a detailed study of the three-body phase space, e.g., B + — > K+ir+n^ . At times, 
study of Q2B final states served as an intermediate step. The choice of analysis 
technique is mostly dictated by the luminosity, expected signal and background, and 
the understanding of the intermediate resonances involved. Table XIV summarizes 
results on the branching fraction and CP asymmetry for various decays with three- 
body mesonic final states. 

A great deal of effort has also been applied to studying charmless three-body 
baryonic decays. Quite often Belle has reported results before its sister experiment, 
BaBar, in these kind of studies. In Table XV we attempt to summarize the results 
obtained in these baryonic decays. An intriguing feature of the results is the peaking 
of baryon-antibaryon pair mass distributions toward threshold. These enhancements 
have generated much theoretical interest. 130 ) 
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Table XIV. Data samples used (N BB ), branching fractions (£>), 90% confidence- level upper limits 
on B (UL), and direct CP asymmetries (Acp) obtained for various charmless decays with three- 
body mesonic final states. 



Final state 


N BS (10 6 ) 


B (10~ 6 ) 


UL (KT 6 ) 


Acp 


Ref. 


K+K+tt- 


85 




2.4 




123) 


K+K-K+ 


152 


30.6 ± 1.2 ±2.3 






107) 


K + K~K° 


85 


28.3 ±3.3 ±4.0 






123) 


K+K--K+ 


85 


9.3 ±2.3 ±1.1 


13 




123) 


K+K° S K° S 


85 


13.4 ± 1.9 ± 1.5 






123) 


K + ir + ir~ nonres. 


386 


16.9 ± Ot^ 






104) 


K+ty+ty- 


386 


48.8 ± 1.1 ± 3.6 




±0.049 ± 0.026 ± 0.020 


104) 


K + ir~-K° nonres. 


85 


r 7+2.7+0.5 
°- ' -2.5-0.4 


9.4 




105) 




85 


36.614 3 ± 3 -° 




±0.07 ±0.11 ±0.01 


105) 


K--K + TT + 


85 




4.5 




123) 


K°K-tt+ 


85 




18 




123) 


K tv -k nonres. 


388 


19.9 ± 2.5l 2 ;o 






103) 


K°TT + TT- 


388 


47.5 ±2.4 ±3.7 






103) 


K° S K° S K° S 


85 


4.2±i;| ± 0.8 






123) 




85 




3.2 




123) 


u)K + n~ nonres. 


657 


5.1 ±0.7 ±0.7 






118) 


(t>(t>K + 


449 


3.2±° ■% ± 0.3 




+0.01±g;i| ± 0.02 


124) 


<jxj)K 


449 


2.3±o ? ± 0.2 






124) 


P°K+tt- 


657 


2.8 ±0.5 ±0.5 






117) 


p% + 7v- 


657 


5.9t3'j ± 2.7 


12 




120) 


fo(980)K + ir- 


657 


1.4 ± 0.4t°; 4 


2.1 




117) 


fo(980)-K + iv- 


657 


0-3±i;l ±0.9 


3.8 




120) 



Table XV. Data samples used (N BB ), branching fractions (B), 90% confidence-level upper limits 
on B (UL), and direct CP asymmetries (Acp) obtained for various charmless decays with three- 
body baryonic final states. 



Final state 


N BB (10 6 ) 


B (10~ s ) 


UL (10~ 6 ) 


Acp 


Ref. 


ppir + 


449 


1.57to.ls ±0.12 




-0.17 ±0.10 ±0.02 


125) 


ppK + 


449 


5-OOto 1 ;^ ± 0.32 




-0.02 ±0.05 ±0.02 


125) 


ppK° 


535 


2.5lto.29 ± 0.21 






126) 


ppK* + 


535 


3.38t_oll ± 0.39 




-0.01 ±0.19 ±0.02 


126) 


ppK*° 


535 


1.18i°;i±0.11 




-0.08 ±0.20 ±0.02 


126) 


pAn° 


449 


3-OOlg-jg ± 0.33 




±0.01 ±0.17 ±0.04 


127) 


pA-K~ 


449 


3.23lo jg ± 0.29 




-0.02 ±0.10 ±0.03 


127) 


pAK~ 


85 




0.82 




128) 


pE°TY- 


85 


3 - 97 ^o.80 ± 0- 56 






128) 


AAn+ 


152 




2.80 




129) 


AAK+ 


152 


2.9ll° ?° ± 0.38 






129) 
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5.2. Radiative penguin decays 

Decay processes of a 6 quark that emit a pho- 
ton are not allowed at the tree level in the SM, 7 
and require a so-called radiative "penguin" loop ^!lyv_ fJ 

(Fig. 26). The dominant contribution in the SM *V 
is from a loop with a top quark and a weak boson. b * I I * s 

However, these heavy SM particles may be re- ^fc. S 

placed by hypothetical particles such charged ^ 
Higgs boson or supersymmetric particles. In such u, d u,d 

a scenario, the decay rate or other observables „. „„ „ ,. r 

big. 26. ieynman diagram for the 

could be drastically modified. Hence, radiative & -> S7 process 

decays have been extensively studied to search 
for and to constrain physics beyond the SM. 

5.2.1. Inclusive B — > X s j measurement 

At the hadron level, the quark-level b — > 57 transition is represented by a ra- 
diative B meson decay into a high energy photon and an inclusive hadronic final 
state with a unit strangeness denoted by the symbol X s . The clean signature of the 
high energy photon makes it possible to measure the decay rate without reconstruct- 
ing the X s . The SM transition rate is calculated including next-to-next-to- leading 
logarithmic corrections to 7% precision. 131 ) 

A large and dominant background to B — > X s j is from the 7r°s (and to a lesser 
extent rf) in the e + e~ — > qq continuum, which subsequently decay into a pair of 
photons. Although this background is several orders of magnitude larger than the 
inclusive photon signal, backgrounds that are not from a B meson decay can be sta- 
tistically subtracted by using the off-resonance data sample taken at 60 MeV below 
the T(4S) resonance. However, since only 10% of integrated luminosity is taken off- 
resonance, this continuum background remains the main source of the statistical and 
systematic error. The remaining backgrounds are from B meson decays, where the 
photon backgrounds are dominantly (in order of their importance) from the 7T°s, 77s, 
radiative decays of other hadrons, final state radiation and electron bremsstrahlung, 
and mis-reconstructed K^s and (anti-)neutrons. The inclusive ir° and 77 produc- 
tion rate from a B meson is directly measured in data, and used to subtract the 
corresponding background contribution, while other sub-dominant contributions are 
subtracted using MC expectations after correcting for the measured data-MC differ- 
ences. The photon energy spectrum, which is monochromatic if b — >■ 57 is strictly a 
two-body process, is broadened by QCD corrections and the Fermi motion of the b 
quark in the B meson. 132 )' 133 ) The measured spectrum in the T(4S) rest frame is 
further broadened by the small momentum of the B meson and the detector reso- 
lution. The branching fraction has to be integrated over the entire photon energy 
range. It becomes more difficult to do so for lower energies as the signal contribution 
becomes smaller and the background becomes insurmountably large. It is now cus- 
tomary to compare the extrapolated branching fraction in the range E^ > 1.6 GeV 
to theoretical predictions. Experimental efforts to lower this bound have been the 
focus of most past B — > X s ^j measurements. Using 657 x 10 6 BB events, Belle 
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measured B -»• X s ~j with £ 7 > 1.7 GeV. 134 ) This should cover (98.5 ± 0.4)% of the 
spectrum above 1.6 GeV. 133 ) The spectrum is shown in Fig. 27 and the branching 
fraction was measured to be 



B(B -»• X sT , £ 7 > 1.7 GeV) = (3.45 ± 0.15(stat) ± 0.40(syst)) x 10" 



statistical 



and 



(5-1) 
systematic. 



30000 




Fig. 27. Photon 
B -> X sl . 



2.6 2.8 
[GeV] 

energy spectrum from 



where the errors are 

Together with BaBar's measurement, the 
world average 28 ^ extrapolated for E 1 > 
1.6 GeV is B(B -> X s7 ) = (3.55 ± 
0.24(exp) ± 0.09(model) x 10" 4 , where the 
first error is a combined experimental (sta- 
tistical and systematic) uncertainty and the 
second is the model error in the extrap- 
olation. This can be compared with the 
theory prediction 131 ) of B(B ->■ X s ^) = 
(3.15 ± 0.23) x 10~ 4 . The results are con- 
sistent, and have been used to constrain 
new physics scenarios. For example, the 
charged Higgs mass is bounded to be above 
295 GeV. 

5.2.2. Exclusive radiative B decays with 
b — > «7 

Exclusive radiative B meson decay 
modes, such as B ->■ K* (892)7, 135 ) have 

been more precisely measured, since one can fully constrain and effectively sup- 
press the background of the decay kinematics using the beam-energy constrained 
mass (Mb c ) and the energy difference (AE). However, theoretical predictions suf- 
fer from large uncertainties in the exclusive form factors, which cannot be reliably 
determined. 136 ) -138 ) 

The B — > K* (892)7 constitutes about 15% of the total B ->■ X s ^ branching 
fraction. Since the B meson has spin zero and the photon has spin one and is 
longitudinally polarized, the X s system cannot be a single kaon (with spin zero), a 
resonance, or an S-wave Kir system. Of the higher kaonic resonances, only B — > 
-fC|(1430)7 139 ) an d B — > Ki(1270)7 140 ) have been measured. In particular, higher 
kaonic resonances around 1.4 GeV have a complicated structure, and among these the 
(1270) contribution was found to be dominant. 141 ** In the multi-body final states, 
many modes have been measured: B -»• Kttttj, 13 ^ B -»• i^7, 142 ) B -»• ifr/7, 143 ) 
B -»• Kpj, 139 ^ B -»• K07, 144 ) and B -»• Apj. 14 ^ 

One way to reduce the theoretical uncertainty is to take ratios or asymmetries. 
In particular, the time-dependent CP asymmetry for a radiative decay into a self- 
conjugate final state has a unique feature. In the SM, the final state, e.g. i^^7r°7, is 

not a CP eigenstate since the photon is dominantly left-handed from B° (with a b 
quark) decay and thus does not mix with the decay from B° (with a b quark) with 
a right-handed photon. The spin flip is suppressed by the quark mass ratio 2m s /mb 
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and hence the time-dependent CP asymmetry is also suppressed in the SM to a few 
per cent. 146 ) Therefore, this asymmetry in the b — > 57 process is sensitive to non-SM 
right-handed currents. 

In 5 ^ i^*(892)°7, the rate to the K%(-> 7r + 7r")7r°7 final state is only 1/9 
of that for K + tt~j. The time-dependent asymmetry is measured by extrapolating 
the Kg momentum from the Kg decay vertex to the interaction region. Therefore, 
the detection efficiency and statistics of the final signal sample are not large. The 
coefficient to the sine term is measured with 535 x 10 6 BB to be 147 ) 

S K , 0l = -0.32+°;! (stat) ± 0.05(syst). (5-2) 

This study can be extended to B° — > P°Q°'j, where P° and Q° are any pseu- 
doscalars, 148 ) or to the P°V°j state if the spin parity of the P°V° system is deter- 
mined. Time-dependent asymmetries have been measured for K s tt°j, K s p°~f, and 
Kg<jry states, although none of them is yet able to constrain the right-handed cur- 
rent. This study is one of the promising modes in the search for physics beyond the 
SM with the high statistics data samples expected at Belle II. 

5.2.3. Radiative B decays with b — > dj 

The b — > dj penguin loop is suppressed with respect to b — > 57 by |Vtd/V^ s | , 
and therefore is sensitive to this ratio. It is particularly interesting because a more 
precise determination of |14d/Vt s | was not available until the B s mixing rate was 
measured 149 - 1 and even after that it provided an independent test of this ratio of 
CKM parameters. 

Since the dominant diagram is suppressed, there are more contributions from 
subleading diagrams. These could lead to a large direct CP violation or large isospin 
asymmetry, although they also modify the determination of |V^/Vi S |. On the other 
hand, time-dependent asymmetry is expected to be even smaller, since the phase 
from V t( i in mixing and b — > dj transition cancel. 146 )' 15 °) Contributions from non- 
SM physics can therefore be relatively enhanced and may be more clearly visible 
than in the b — > S7 case. 

Because of the similarity of the kinematics, the large b — > 57 process is a severe 
background to the suppressed b — > dj process. In the reconstruction of an exclusive 
decay mode, particle identification devices are crucial to separate the kaon in b — >■ S7 
from the pion in b — > dj. Exclusive 6—7-^7 decay modes such as B — > pj and 
B W7 have been searched for since the start of Belle, and were finally observed 
with 386 x 10 6 BB pairs in a combined measurement. 151 ) Charged and neutral 
modes are combined assuming isospin symmetry B{B + — > p + j) = 2B(B° — > p°j) 
and B(B° ->• p ^) = B(B° ->• ^7). The latest result with 657 x 10 6 BB pairs is 
shown in Fig. 28 and the combined branching fraction is measured to be 152 ) 

B(B -»• (p,w)7) = (1.14 ± 0.20(stat)+°;}°(syst)) x 1(T 6 . (5-3) 

As B — > p7 is suppressed compared to B — > K*^ by |T4d/^4s| 2 > known kine- 
matic corrections, and less-known form factor ratios and corrections for subleading 
diagrams, the result is combined with a corresponding analysis on B — > K*^j to 
constrain |V^/Vt s |. The result is 




\Vtd/V ts \ = 0.195^° (exp) ± 0.015(theo) (5-4) 

where the first error is a combined statistical and systematic uncertainty and the 
second error is the theory uncertainty on the ratio. 

The B — y signal is found to be stronger than B + — > p + 7- This corre- 
sponds to a large isospin asymmetry, which is defined as A(pj) = 2 ^ B ° + B(B + — > 

p + j)/B(B° — > j0°7) — 1. The isospin asymmetry is calculated as 

A(py) = -0.48±g;fJ(stat)±g;^(syBt). (5-5) 

BaBar also measures this ratio and finds the same tendency; the combined isospin 
asymmetry is ~ 3a away from the SM expectation, which could be at most ~10%. 
As the statistical error is still large, the high statistics expected at Belle II will be 
necessary to clarify this tension. 

5.3. Electroweak penguin decays 



i- i- 



h 




u.d u.d u.d u,d 

Fig. 29. Feynman diagrams for the b —¥ sl + l~ process. 

The b — > s(d)l + l~ transitions proceed at lowest order in the SM via Z/j pen- 
guin diagrams or a W box diagram (Fig. 29). The b — > s{d)viVi transitions also 
proceed through similar diagrams except for the 7 penguin diagram. NP mediated 
by SUSY particles or a possible fourth generation may contribute to the penguin loop 
or box diagram and as a result branching fractions and other properties could be 
modified. 153 ^ Such NP contributions may change the Wilson coefficients that param- 
eterize the strength of the short distance interactions. This is similar to b — > s(d)j, 
but has a richer structure. The decay B — > K*l + l~ is of particular interest since 
its large branching fraction facilitates the examination of various observables that 
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are sensitive to NP. For instance, the lepton forward-backward asymmetry (Afb), 
the K* polarization (.Fl), and the K*l + l~ isospin asymmetry (Aj) as functions of 
dilepton mass squared (q 2 ) differ from the SM expectations in various NP models. 154 ) 
The neutral pure leptonic decays B° — > l + l~ and B° — > v\pi proceed mainly 
through the box and Z boson mediated annihilation diagrams, which are equiva- 
lent to the diagrams for b — > dl + l~ and b — > dviyi- In the SM these decays are 
also helicity suppressed and, compared to the charged purely leptonic decays (see 
Sect. 4.3), the branching fractions are about three orders of magnitude smaller for 
the corresponding generation. 155 ) The SM branching fraction of B° — > v{Vi is at 
the level of icr 20 . 156 ) The lepton-flavor-violating decay B° — > e ± fi T is not an elec- 
troweak penguin decay and is forbidden in the SM, but can occur in the Pati-Salam 
model 157 -* or supersymmetric models, 158 ) and can be searched for simultaneously. 
A positive signal for any of these decay modes with the current Belle data sample 
would demonstrate NP in the loop. 

5.3.1. Exclusive 6 — » s(d) l + l~ decays 

The study of the decay B — > K^l + l~ started at the beginning of Belle and was 
updated several times. We reported the first observations of B — > Kl + l~ 159 ^ and 
B -»• K*l + l~ 16 ^ with 31.3 x 10 6 and 152 x 10 6 BB pairs, respectively. In 2006, 
Belle published the first measurements of the forward-backward asymmetry and the 
ratios of Wilson coefficients A 9 /A 7 and A w /A 7 using 386 x 10 6 BB pairs. 161 ) An 
unbinned maximum likelihood fit to q 2 and cos#; was used to extract the ratios of 
the Wilson coefficients, where 9i is the angle between the momenta of a negative 
(positive) lepton and the B (B) meson in the dilepton rest frame. 

The latest analysis in 2008 162) used 657 x 10 6 BB pairs; more observables were 
measured. Candidate B — > K^l + l~ decays were reconstructed in 10 channels: 
K + n~, KgTT + , K + ir° for K*, K + , and Kg for K, with e + e~ and lepton pairs. 

The dilepton mass of each candidate was required to be outside of the J/ip and 
tp(2S) mass regions to avoid the large charmonium background, and above the 7T° 
mass for e + e~ pairs to avoid the 7r° Dalitz decay, photon conversion, and the pole at 
q 2 = 0. Two major backgrounds were considered: the continuum and BB events in 
which both B mesons decay semileptonically. These backgrounds were suppressed 
by imposing requirements on the signal-continuum and signal-Si? likelihood ratios. 

After requiring the candidate AE to lie in the signal region, the signal yields in 
each q 2 bin were extracted from an unbinned likelihood fit to Mb c and the Ktt mass 
(Mk-k) for the K*l + l~ mode and Mb c only for the Kl + l~ mode. The corresponding 
branching fractions were thus obtained. The Fl and Afb parameters were extracted 
from fits to cos 9k* and cos 6*; in the signal region, where 9k* is the angle between 
the kaon direction and the direction opposite to the B meson in the K* rest frame. 
The signal PDFs for the cos 9k* and cos 9\ variables are the product of the following 
two functions, 



[If l cos 2 9 K * + \{l-F L ){l- cos 2 9 K *)] x e(cos 9 K *) 
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and 

\-F L (l- cos 2 6 l ) + l(l-F L )(l + cos 2 9i) + Afb cos 0,] x e(cos 0,), 
4 o 

where e(cos0R-*) and e(cos0;) are the reconstruction efficiencies. For the B — > Kl + l~ 
modes, Fl is set to 1. Furthermore, this analysis also reported the isospin asymmetry 
defined as 

_ (t b +/t b o) x B(K^°l+l-) - BjK^l+r) 
1 ~ (t b+ /t b0 ) x B(K(*)oi+l-) +B(K(*)±l+l-y 

where t b +/t b o is the ratio of B + to B 
lifetimes. These observables were measured 
for the first time in six q 2 bins as shown 
in Fig. 30. Although the uncertainties in 
the Afb values are still large, the posi- 
tive central values in all q 2 bins suggested 
a non-zero Afb(q 2 )- This phenomenon 
would have been an undeniable signature of 
NP, but unfortunately did not persist with 
larger data samples at the LHC hadron col- 
lider. 163 ) Two more observables, the direct 
CP- violating asymmetry and the lepton fla- 
vor ratio of the muon to electron modes, 
were also measured. The latter is sensitive 
to Higgs emission and could be larger than 
the SM expectation in the two Higgs dou- 
blet model at large tan/3. 164 ) 

The observed values are Acp(K*l + l~) = 
-0.10 ± 0.10 ± 0.01 and A C p{Kl + l-) = 
0.04 db 0.10 ± 0.02, consistent with no asym- 
metry, and R K * = 0.83 ± 0.18 ± 0.08 and R K = 1.03 ± 0.19 ± 0.06, similar to the 
SM values. The measurements of so many observables demonstrate the richness and 
potential of the B -> K^l + l~ decay. 

A search for the exclusive b — > dl + l~ process, B — > irl + l~ (ir = ir + or tt°), was 
performed using 657 x 10 6 i?-B pairs. 165 ) No obvious signal was observed and upper 
limits on the branching fractions at the 90% C.L. were obtained: B(B + — > ir + l + l~) < 
4.9 x 10~ 8 and B(B° — > ir°l + l~) < 15.4 x 10~ 8 . These limits are approaching the 
SM expectations, which are O(10 -8 ). 

5.3.2. Inclusive B — > X s l + l~ decay 

The inclusive measurement of the b — > sl + l~ process is experimentally chal- 
lenging, but can be compared with theoretically clean predictions. The standard 
technique is to analyze B — > X s l + l~ events with a semi-inclusive approach, where 
the X s is reconstructed in 18 different combinations of either a K + or Kg combined 
with to 4 pions, of which up to one ir° is allowed. This set of final states covers 
around 62% of X s decay states. The missing states were taken into account in the 
signal efficiency obtained from MC simulations. 




10 12 

q 2 (GeV 2 /c 2 ) 

Fig. 30. The q 2 dependence of Fl (top), 
Afb (middle), and Ai (bottom) in 
six bins. The results for B — > K*l + l~ 
are the filled circles and those for B 
Kl + l~ are open circles (Ai only). 
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The first observation of B — > X s l + l was reported by Belle using 65.4 x 10 6 BB 
pairs in 2003. 166) The latest Belle results in 2009 used 657 x 10 6 BB pairs. 167) As 
in the exclusive analysis, signal candidates were selected with AE and then the M^ c 
distribution is used to extract the signal yield. The dilepton mass was required to 
be outside of the J/ip and ip(2S) regions and the low mass region below 0.2 GeV/c 2 . 
The signal yields were extracted from an unbinned maximum likelihood fit to the 




0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 u 5 10 15 20 25 



Mass(Xs) (GeV/c 2 ) q 2 (GeV/c 2 ) 2 

Fig. 31. The top two plots show the M^ c distributions with the fit curves superimposed for the 
entire sample (left) and for 1.0 GeV/c 2 < Mx s < 2.0 GeV/c 2 (right). Points with error bars are 
data; the dominant background, peaking background, and self-cross-feed components are the 
yellow, green, and blue solid shaded areas, respectively. The bottom plots show the dB/dMx s 
(left) and dB/dq 2 (right) distributions for data (points) and the SM expectation (histograms). 

Afb c distribution. In addition to the dominant backgrounds from BB pairs and 
continuum, an effort was made to investigate the peaking background and include it 
in the fit. Two kinds of peaking background were considered: charmonium peaking 
background and hadronic peaking background. The former includes the residual of 
B — > J/ipX s ,B — > ip(2S)X s events after the J/ip and tp(2S) vetoes, and a possible 
contribution from higher ip resonances, such as the ^(3770), ^(4140), and ^(4160). 
The hadronic peaking background contains B — > X s hh and B — > X s hlv events in 
which one or two hadrons are misidentified as leptons. The peaking backgrounds were 
estimated directly from data or simulations and the corresponding yields were fixed 
in the fit. Finally, the last component considered is the self-cross-feed background, 
in which the B daughter particles are not correctly selected. Its probability density 
function was modeled as a histogram with the ratio of the normalization of the 
self-cross-feed background to signal fixed to the MC simulation value in the fit. 
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The probability density function for the dominant background was modeled by an 
ARGUS function with the parameters floated in the fit. A simultaneous fit to the 
X s l + l~ and Xge^ 1 ^ samples was performed with the same ARGUS parameters for 
the dominant background. 

The branching fractions of B — > X s l + l~ were reported as a function of Mx s and 
q 2 separately. Fit results for the total sample and a subset with Mx s > 1.0 GeV/c 2 
are shown in the top two plots of Fig. 31, and the differential branching fractions 
as functions of Mx 3 and q 2 are shown in the bottom two plots. The differential 
distributions are compared with the SM expectation, 168 ) and found to be in good 
agreement. The branching fraction of B — > X s l + l~ in the entire Mx s range was 
obtained by summing the branching fraction in each Mx 3 region and correcting 
for the X s l+l~ fraction in the J/ip, i/)(2S), and M Xs > 2.0 GeV/c 2 regions. The 
branching fraction with the dilepton mass above 0.2 GeV/c 2 is thus measured to be 
B(B -> X s l + l~) = (3.33±0.80lg;^) x 10~ 6 . We also reported the branching fractions 
separately for the electron and muon modes using the same analysis procedure, 
B(B -»• X s e + e~) = (4.56 ± 1.15±g;^) x 10~ 6 and B(B -»• X s n+iT) = (1-91 ± 
1.02±g;l|) x lO" 6 . 

5.3.3. Searches for B° decays to invisible final states 

Searches for B° decays to invisible final states are rather challenging. The 
same strategy used in the B + — > t + v t analysis was applied to identify the sig- 
nal (Sect. 4.4.2). Candidate events were selected by fully reconstructing a B° meson 
and requiring no additional charged, 7r°, or K\ particles in the rest of the event. 
The signal can be identified by requiring no or very little extra calorimeter energy 
(-E'ecl) in the event. Furthermore, two variables were used to distinguish the signal 
and the continuum background: cos 9b and cos 9t, where the latter is the cosine of 
the angle of the Stag thrust axis with respect to the beam axis in the CM frame. 
The continuum was suppressed by making requirements on cos 9t and cos 9b- 

The signal yield was extracted from an unbinned extended likelihood fit to -Eecl 
and cos 9b- Candidate events in the fit were categorized as signal, BB, and non-S 
backgrounds, where the latter includes the continuum and a small e + e~ — > t + t~ 
background. Using a sample of 657 x 10 6 BB pairs, the signal yield obtained in the fit 
was 8.9l5;5 events. Since no significant signal was observed, we provide the branching 
fraction upper limit including systematics at the 90% C.L. of B(B° — > invisible) < 
1.3 x 10~ 4 . 174) The expected upper limit from the MC study is 1.1 x 10~ 4 . 

5.3.4. Search for B° -> l+l~ 

The results of searches for the decays B° — > e + e~, and e ± / u T (collectively 

denoted by B° l + l~) were reported at the beginning of Belle using only 85 x 10 6 
BB pairs. 175 ) Since the background for the two energetic leptons is relatively small, 
the Belle analysis was able to suppress the background effectively while maintaining 
a high reconstruction efficiency. After all the selection criteria, no events were found 
in any of the three modes. 175 ) The upper limits are: B(B° — > e + e~) < 1.9 x 10~ 7 , 
B(B° -»• < 1.6 x 10" 7 , and B{B° -»• e ± /i =F ) < 1.7 x 10~ 7 . Furthermore, 

a lower bound on the mass of the Pati-Salam leptoquark model of 46 TeV/c 2 was 
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obtained using the upper limit for the e\x mode. 

§6. Tau physics 

The tau lepton is an extremely convenient probe to search for NP beyond the 
SM because of the well-understood mechanisms that govern its production and decay 
in electroweak interactions. With its large mass, it is the only lepton that can 
decay into hadrons, thus providing a clean environment to study QCD effects in 
the 1 GeV energy region. Tau physics at Belle is categorized by two themes; NP 
searches and SM precision measurements. To probe NP, we search for lepton-flavor 
violating (LFV) decays, CPV in the charged lepton sector, and the electric dipole 
moment (EDM) of the tau lepton. For SM precision measurements, we measure the 
r lepton mass, the branching fractions of various hadronic decay modes, and their 
invariant mass distributions. In this section, we summarize the results obtained 
from the world's largest data sample (about 10 9 t+t~ pairs) accumulated at the 
Belle experiment. 

6.1. New physics searches 

6.1.1. Tau lepton flavor violation 

An observation of LFV would be a clear signature of NP since LFV in charged 
leptons has a negligibly small probability in the SM, O(10 -54 ) — O(10 -52 ), even 
if neutrino oscillations are taken into account. 176 -* Since the r is the most massive 
charged lepton, it has many possible LFV decay modes. Belle has examined as many 
decay modes as possible in the LFV searches, since the specific mechanisms of NP 
are unknown. 

Models including supersymmetry (SUSY), which is the most popular scenario 
beyond the SM, can naturally induce LFV at one loop. In many SUSY models, 
including see-saw extensions and grand unified theories, r — > fx'y is expected to have 
the largest branching fraction of all the possible r LFV decays. In some cases, 
however, such as the Higgs-mediated scenario, r decay into \xr\ or \i\x[x can become 
more probable. By measuring the branching fractions for various r LFV decays, 
one may be able to determine the NP model favored by nature. Among the various 
modes studied in Belle, we focus here on three possibilities, r — > £j, ££'£", and £P°, 
where £ stands for e or \i and P° is 7r°, 77, or rf . 

In an LFV analysis, in order to evaluate the signal yield, two independent vari- 
ables are used: the reconstructed mass of the signal and the difference between the 
sum of energies of the signal r daughters and the beam energy (AE) in the CM 
frame. In the r — > fij case, these variables are defined as 

M M7 = ^E^-P^ (6-1) 

AE = E^ - -Ebe^m, (6-2) 

where E m and P m are the sum of the energies and the magnitude of the vector sum 
of the momenta for the /j and the 7, respectively. The superscript CM indicates 
that the variable is defined in the CM frame, e.g. E™ m ^ s * ne beam energy in the 
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CM frame. For signal, M^~ f and AE should be in the vicinity of M„~ ~ m T and 
AE ~ (GeV), while for the background, M^ 7 and AE will smoothly vary without 
any special peaking structure. Taking into account the resolution of the detector and 
the correlation between M^ 7 and AE, we use an elliptical signal region. To avoid 
bias, we perform a blind analysis: the data in the signal region are blinded when 
determining the selection criteria and the systematic uncertainties. After fixing these 
quantities, we open the blind and evaluate the number of signal events in the signal 
region. 

r — > £7 



We have searched for r — > £7 with a data set corresponding to produced 4.9 x 10 s 
t + t~ pairs. 177 ) The main background (BG) is from r — > £v(i> T + extra 7 events 
and radiative di-muon (for ^7) or Bhabha (for &y) events. The observed M^—AE 
distributions are shown in Figs. 32(a) and (b) for r — > [vy and r — > &y, respectively. 
The signal yield is evaluated from an extended unbinned maximum-likelihood fit to 
the M^-AE distribution. We found no excess in the signal region. We thus obtain 
an upper limit on the branching fraction for t — >■ ^7 (e^y) of 4.5 x 10~ 8 (1.2 x 1CP 7 ) 
at 90% C.L. 



-/►■■■>■ ■ ■ 




A/ mv (GeV/c 2 ) 



1.7 1.75 l.i 
M. (GeV/c 2 ) 



Fig. 32. M^-AE distributions in the search for (a) r — > /17 and (b) r — > ej. 177 ' The black dots 
and shaded boxes show the data and signal MC, respectively, and the ellipse is the 2a signal 
region. 



r ->• we" 

The decays r — > ££'£" have been searched for with nearly the entire data sample 
of 7.2 x 10 8 t + t~ pairs obtained by Belle. 178 ) Figures 33(a) and (b) show the three- 
lepton invariant mass versus AE {M^-AE) distributions for the t~ — > e~e + e~ 
and t~ — > [i~ candidates after selection, respectively. No events in the signal 
region have been found in any of the six modes; the 90% C.L. upper limits on the 
branching fractions in units of 10 -8 are given in Table XVI. The obtained upper 
limits are two or three times more restrictive than those obtained previously. 179 ) 

r -)■ £P° (P° = TT°,rj,r]') 

Early Belle results on the search for the r decays into a lepton and a neutral 
pseudoscalar (vt°, 77, 7/), were based on a data sample of 3.6 x 10 8 t + t~ pairs; the 
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Fig. 33. Mtu-AE distributions for (a) r~ — > e~e e~ and (b) t~ — > /x~/i + /i~ modes. J The 
black dots and shaded boxes show the data and signal MC, respectively. The ellipse is the 
signal region. The region formed by the two parallel lines, excluding the signal ellipse region, is 
the side-band region used to evaluate the expected number of backgrounds in the signal region. 



Table XVI. Summary of the efficiency (Eff.), the expected number of BG events {N^q), and the 
upper limit on the branching fraction (UL) at 90% C.L. for r~ — > £~ £' + £"~ . 



Mode 


Eff.(%) 




UL (10" s ) 


Mode 


Eff.(%) 




UL (KT S ) 


e~e + e~ 


6.0 


0.21 ±0.15 


2.7 


e~ fj, + 


6.1 


0.10 ±0.04 


2.7 


e~e + fi~ 


9.3 


0.04 ± 0.04 


1.8 


fi~ e + fj,~ 


10.1 


0.02 ±0.02 


1.7 


e~ )i + e~ 


11.5 


0.01 ±0.01 


1.5 




7.6 


0.13 ±0.06 


2.1 



resulting upper limits were in the range (0.8-2.4) x KT 7 at 90% C.L.. 180 ) Recently, 
we have updated the results with a data set two times larger. By studying the 
backgrounds in detail, we obtain on average about 1.5 times higher efficiency than in 
our previous study while maintaining a background level in the signal region of less 
than one event in all modes. The results are summarized in Table XVII. A single 
event is found in r — > erj(—> 77) while no events are observed in other modes. The 
obtained 90% C.L. ULs on the branching fraction are in the range (2.2 — 4.4) x 10~ 8 . 



Table XVII. Summary of the efficiency (Eff.), the expected number of BG events (N^q), and the 
upper limit on the branching fraction (UL) for r — > £P°, where (comb.) means the combined 
result from subdecay modes. 



Mode 


Eff.(%) 




UL (Hr s ) 


Mode 


Eff.(%) 


Nbg 


UL (HT S ) 


w(-> 77) 

fJ,T](—¥ 7r7T7r°) 


8.2 
6.9 


0.63 ±0.37 
0.23 ±0.23 


3.6 
8.6 


er?(-> 77) 
ei](— > 7T7r7r ) 


7.0 
6.3 


0.66 ± 0.38 
0.69 ±0.40 


8.2 
8.1 


/iry(comb.) 






2.3 


er?(comb.) 






4.4 


fjbf] {— > TVTTTl) 

w'(-> 7P°) 


8.1 
6.2 


u - uu -o.oo 
0.59 ±0.41 


10.0 
6.6 


er]' {_—¥ ixirq) 
er/(-> 7p°) 


7.3 
7.5 


0.63 ±0.45 
0.29 ± 0.29 


9.4 

6.8 


/1J7' (comb.) 






3.8 


er/(comb.) 






3.6 


M 7T U 


4.2 


0.64 ±0.32 


2.7 


evr 


4.7 


0.89 ± 0.40 


2.2 



In total, Belle has completed searches for 46 lepton-flavor-violating r decay 
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Fig. 34. Current 90% C.L. upper limits for the branching fraction of r LFV decays obtained 
in the CLEO, BaBar, and Belle experiments. Red, blue, and black circles show Belle, 
BaBar, and CLEO results, respectively. 



modes using nearly the entire data sample of 1000 fb _1 except for the ongoing anal- 
ysis of the £7 modes. No evidence for LFV decays has been observed in any of the 
modes studied and we set 90% C.L. ULs on the branching fractions at the O(10~ 8 ) 
level. The current status of r LFV searches in P-factory experiments and in the 
CLEO experiment is summarized in Fig. 34. The sensitivity for LFV searches has 
been improved by two orders of magnitude in comparison with the CLEO results. 
This is due to the effective background rejection as well as the increase in the size 
of the data sample. In the near future, SuperKEKB/Belle II at KEK will reach a 
sensitivity at the O(10 _9 )-O(10 -10 ) level and explore a wider region of parameters 
in various NP scenarios. 

6.1.2. CP-violating r decays 

To date CPV has been observed only in the K and B meson systems. In the 
SM, all observed CPV effects can be explained by a single irreducible complex phase 
in the CKM quark mixing matrix. It is important to look for other CP-violating 
effects where SM CPV is not expected in order to find NP. One such system is the 
r lepton. In hadronic r decays, no CP-violating effects from the SM are expected 
except for cases in which the decay products contain K$ mesons. In other words, the 
r decay is an ideal place to look for other CP- violating effects that could originate 
from new physics scenarios, such as the minimal supersymmetric model 181 ) or from 
multi-Higgs-doublet models. 182 ) 

If there is a CP-violating NP amplitude in a r decay, interference between the 
SM and NP amplitudes should occur. Even in this case, as was emphasized by 
J.H. Kiihn and E. Mirkes, 183 ) one cannot observe the CPV effects as a difference 
of the total decay rates between r~ and r + , but instead one needs to measure the 
difference between the decay-angular distributions of the hadronic system for r~ and 
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r + . The analysis of the decay-angular distribution is therefore crucial. 

We searched for CP violation in — > Kgir^Ur using a 699 fb _1 data sample. 184 * ) 
In order to search for CP violation in the angular distribution, we define the CP 
asymmetry observable as the difference between the mean value of the product of 
the decay angles cos /3 cos (j) in the K^tt^ system for r~ and r + : 



iCP 



(cos /3 • cos tp) T - — (cos (3 ■ cos ip) T -t 



where j3 (ip) is the angle between the direction of the (r) and the direction of 
the e + e~ CM system measured in the K^-k^ rest frame. 



We obtain 3.2 x 10 5 — >■ K%k ± v t candidates. 



The K^tt^ invariant mass 



distribution shown in Fig. 35 clearly indicates that, in addition to the if* (890) 
resonance, other resonant contributions are also needed to explain the full spectrum 
(see Sect. 6. 2. 3 for more details). The measured CP asymmetry A CP is shown in 
Fig. 36 as a function of Kgir^ invariant mass after correcting for known detector 
effects. The result indicates that there is no CP asymmetry at the 1% level. Then 
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Fig. 35. Invariant mass spectrum of the K%ir 
system inr-> K%-k^v t candidates. 184 ' 



W (GeV/c 2 ) 

Fig. 36. Measured CP-violating asymmetry 
A CF as a function of the Kgit invariant 
mass W after subtraction of background 
(black squares). 184 - 1 The inverted red tri- 
angles show the expected asymmetry when 
Q(Vs) = 0.1. Note that the previous CLEO 
result 185 ' corresponds to Ss(?7s) < 0.19. 

we obtain the upper limit for the CP-violating scalar coupling constant t/ s 184 ^ to be 

|Im(r/ s )| < (0.012 - 0.026), 

at the 90% C.L. Our study achieved ten times higher sensitivity than the previous 
CLEO results shown by the inverted red triangles in Fig. 36. 

6.1.3. Tau electric dipole moment 

If an elementary particle has a non-zero electric dipole moment (EDM), this 
is a clear indication of violation of the T-reversal symmetry and thus violation of 
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CP invariance according to the CPT theorem. The current limit for the r EDM 
(d T ) is several orders of magnitude less restrictive than that for the electron, muon, 
or neutron. Measurement of d T is difficult because of r's short lifetime. However, 
improvements in sensitivity are interesting both theoretically and experimentally. 
As explained below, one can measure the r EDM by using the correlation of decay 
product momenta in the process e + e~ — > t + t~ . 

The matrix element for the process e + e~ — > t + t~ , is given by the sum of the 
SM term Mg M , the EDM term \d T \ 2 M. 2 d2 , and the interference between them: 

M 2 = M 2 SM + Re(d T )Ml c + ^(dr)M 2 m + \d T \ 2 M 2 2 , 

where Re(d r ) ( Im(d T ) ) is the real (imaginary) part of the EDM. These interference 
terms -M-L/lm contain the following combination of spin-momentum correlations: 



M 2 Rc oc 



(S + xS-)-k, (S+xS-)-p, 

M 2 lm *(S + -S-)-k, (S+-S-)-p, 
where S± is a spin vector, and k and p are the unit vectors of the r" 



and 



momenta in the CM system, respectively. These terms are CP-odd since they change 
sign under a CP transformation. 

One could evaluate the value of the matrix elements if the values of S± and k 
could be measured on an event-by-event basis from the r-decay products. Although 
one cannot know them completely due to missing neutrinos from r decays, one can 
obtain the most probable values of S± and k by calculating approximate averages 
from measurements of the momenta of r decay products. In the analysis, we employ 



the method of 

In this method, the observables Or c and 

Oim are 



optimal 



observables. 186 ) 



0Re = 



ML 
M 2 SM ' 



o 



Im 



M 2 lm 
M 2 SM 



evaluated using the most probable values of 
S± and k. The means of 0R e , Cim are pro- 
portional to the EDM value and have max- 
imal sensitivity. The relation between the 
mean values and the EDM d T is shown in 
Fig. 37 for the r + r~ — > (itv T )(pv T ) mode. 

We carried out the EDM analysis with 
a 29.5 fb _1 data sample collected by the 
Belle detector. 187 ) In order to obtain 
the maximal sensitivity, we measured the 
EDM in 8 modes, t + t~ — > (ev e v T )([iv^v T ), 
(ev e v T )('Ku T ), (^v^v t ){ttv t ), (ev e v T )(pp T ), 
(/j,^u T )(pu T ), (irv T )(irv T ), (-KU T )(pv T ), and 
(pv T )(pv T ). The values of EDM obtained 





1 




> 


0.8 


a 


0.6 




0.4 




0.2 












o 


-0.2 






-0.4 




-0.6 




-0.8 




-1 

























• 


Re(d T 


) 








o 


lm(d T ; 


















































































-4 


2 'c 
Re/lm(cy 


(10" 16 


\ 

ecm) 



Fig. 37. Correlation between the r EDM 
and the optimal observable obtained 
by MC simulation for t + t~ — > 
{i\i> T ){pv T )} &7 ^ Black dots and circles 
indicate the relations for the real and 
imaginary parts, respectively. 
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from the mean values of the optimal ob- 
servables are shown in Fig. 38. All results 
are consistent with zero EDM. 

We obtain mean values for Re(d T ) and Im(d T ) by taking the weighted mean of 
8 modes to be 

Re(d T ) = (1.15 ± 1.70) x l(T 17 ecm, Im(d r ) = (-0.83 ± 0.86) x 10~ 17 ecm. 
The 95% C.L. intervals are 

-2.2 xl0~ 17 < Re(d T ) < 4.5 x 10" 17 ecm, -2.5xl0" 17 < lm(d T ) < 0.8 x 10~ 17 ecm. 
These limits are ten times more restrictive than previous experiments. 




Fig. 38. Results on the r EDM for 8 modes and the weighted mean for the (a) real and (b) imaginary 
parts. 

6.2. Precision measurements 
6.2.1. Tau lepton mass 

A precise measurement of the r lepton mass is very important to test electroweak 
theory and lepton universality, since the decay width is proportional to the r lepton 
mass to the fifth power. For a long time the world average for the tau mass was 
dominated by a single precise measurement carried out at the e + e~ threshold by the 
BES experiment in 1996. 

Belle measured 188 ^ the r mass by using a pseudomass method and showed that 
a precision similar to that obtained in the threshold region can be obtained with 
completely different systematic errors. Figure 39 shows the pseudomass distribution 
obtained by Belle, where a few hundred thousand — > -k^-k^-k+Vt- events are used. 
The pseudomass is defined by 

M min = VMJ + 2(£ bcam -E X )(E X -P X ), 

where M X ,P X ,E X are the mass, absolute momentum, and energy of the 3ir system, 
respectively. We obtain 

m T = (1776.61 ± 0.13(stat) ± 0.35(syst)MeV/c 2 
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ARGUS 1992 

1 776.30 ± 2.40 ± 1 .40 

BES 1996 
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Fig. 39. Pseudomass distribution M m - m for 
-)■ 371-*^ candidates. 188) 



Fig. 40. Summary of the r mass measure- 
ments. 



for the r mass and 

\m T + — m T - \/m T < 2.8 x 10~ 4 at 90% C.L., 

the most stringent limit for the relative mass difference between positive and negative 
r leptons. 

Measurements with a similar precision were subsequently carried out by BaBar 
with the pseudomass method 189 ) and by KEDR with the threshold-scan method. 19 °) 
The current status of r mass measurements is summarized in Fig. 40. 

6.2.2. Spectral function in — > tt^tt ^ decay 

Among the decay channels of the r lepton, r ± — > 7r ± 7r°^ T has the largest branch- 
ing fraction. From the conservation of vector current (CVC), the 7r _ 7r° mass spec- 
trum can be related to the cross section for the process e + e~ —> tt + it~ and thus can 
be used to improve the theoretical error on the anomalous magnetic moment of the 
muon a M = (g M - 2)/2. 

Using a sample of 5430000 — > 7r ± 7r°z^ T , Belle measures the branching fraction 
and the 7T7r° mass spectrum, 191 ) which are important for obtaining the theoretical 
value of — 2. After unfolding using the singular-value-decomposition method, 58 ' 
the 7T7r° mass spectrum obtained is shown in Fig. 41, where the shapes for p(770), 
^(1450), and />(1700) resonances and their interference pattern are measured very 
precisely. Figure 42 is the pion form factor in the p(770) region obtained from the 
mass spectra in Fig. 41. The measured branching fraction is 

0(7-+ -> tt+ttV) = (25.24 ± 0.04(stat) ± 0.40(syst))%. 

It is known that there is a significant difference in the value of obtained 
from e + e~ and r data. A lengthy discussion is ongoing about a possible source 
of this difference. 192 ^ Belle r data are in very good agreement with the recent 
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measurement of the e + e~ —> ir + Tr~ cross section from initial-state radiation (ISR) 
data by BaBar. 193 ) In addition, it has been pointed by F. Jegerlehner 194 )' 195 ) that 
7 — p interference, which is present only in e + e~ and does not contribute in the r 
decay, plays an important role. After taking this interference effect into account, 
the discrepancy between the e + e~ and r data disappears; i.e. the hadronic term of 
(g^ - 2) from the e + e~ data is a)f d [e] = 690.8(4.7) x KT 10 , while including the r 
data it becomes a£ ad [e, r] = 691.0(4.7) x 1CT 10 . 194 ) Note that without the p - 7 
interference correction, a^ ad [e, r] was a^ ad [e, r] = 696.6(4.7) x 10~ 10 . 

The resulting difference between theory and experiment for is greater than 
3d, which strengthens the difference further. Recently there have been efforts to 
evaluate a^ ad in lattice QCD. 196 **' 197 ) The reported values scatter in the range from 
a had = 641 x 1Q -io to 748 x 10 -io with an error f (30 _ 6 4) x 10 -io The error ig 

one order of magnitude larger than that obtained so far from e + e~ and/or r data. 




(M «) 



Fig. 41. Unfolded it ± -k mass spectrum for 
Solid circles are the data 
and the solid curve is a fit. The error bars 
include both statistical and systematic un- 
certainties. 
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Fig. 42. Pion form factor |F^(s)| 2 in the 
/o(770) region extracted from the mass 
spectra in Fig. 41. 



191) 



6.2.3. Observation of decays with three kaons 

Using a data sample of 401 fb _1 corresponding to 3.6 x 10 8 t + t~ pairs, Belle 
reported the first observation of decays with three charged kaons in the final state. 198 ^ 
We select events in which a K + K~ pair comes from the </> meson and, after taking 
into account a serious peaking background from r~ —> §-k~v t , report the branching 
fraction, B(t~ — > 4>K~v T ) = (4.05 ± 0.25 ± 0.26) x 10~ 5 . In addition, we observe 
r~ — > cj)iT~v T and — > (f>7r~ tt° v t decays, which is a serious peaking background for 
the three kaon process. Later BaBar confirmed the existence of this decay with a 
branching fraction consistent with ours. 199 ) 
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6.2.4. Study of r" -> K s k~v t 

A data sample of 351 fb _1 has been used to study the Ksis~v T final state. 20 °) As 
a result of the analysis, 53110 lepton-tagged signal events have been selected. The 
measured branching fraction, B(t~ -)■ K s n~v T ) = (0.404 ± 0.002 ± 0.013)%, is the 
most precise of all the published measurements. Although the Belle result is consis- 
tent with the other results within errors, the central value is somewhat lower than all 
of them. 
An analysis of the Kstt~ invariant mass 
spectrum shown in Fig. 43 reveals the dom- 
inant contribution from if* (892) ~ with ad- 
ditional contributions of higher states at 
1400 MeV. A satisfactory fit is obtained 
only if the existence of a broad scalar state, 
Kq(800), is assumed. For the first time the 
K* (892) ~" mass and width have been mea- 
sured in r decay: M = (895.47 ± 0.20 ± 
0.44 ± 0.59) MeV, T = (46.2 ± 0.6 ± 1.0 ± 
0.7) MeV, where the third uncertainty is 
from the model. The K*(892)~ mass is 
significantly higher than the world aver- 
age value based on various hadronic exper- 
iments and is much closer to the world av- 
erage for the neutral if* (892). 

6.2.5. Measurement of hadronic r decays 
with an r] meson 

Using a data sample of 490 fb" 1 we 
have studied hadronic r decay modes in- 
volving an T) meson. Candidate r\ mesons 
are reconstructed from their decays into 77 
and 7r + 7r~7r°. 201 ) Table XVIII lists the measured branching fractions or the up- 
per limits. In all cases the number of observed events is significantly higher and 
the results are more precise than previous measurements by CLEO 202 ^ 204 ) and 
ALEPH. 205 ) For the K~ ~r\r\v T decay mode, our result is the first measurement. For 
7r^-K°r]v T the invariant mass spectrum and the branching fraction are consistent with 
a prediction based on the conserved vector current (CVC) theorem. 206 ) 

6.2.6. Decays with three hadrons in the final state 

With a data sample of 666 fb" 1 Belle has also studied various decay modes of the 
r lepton with three hadrons in the final state. 207 ) The results of this analysis for the 
branching fractions of various three-prong modes are listed in Table XIX together 
with recent results from BaBar. 1 ") Note that, for the 7r~ir + ir~i , T and K~n + 'K~v T 
modes, the branching fractions listed do not include any K° contribution, while the 
result for K~ K + K~v T includes 4>K~v T . 

In Fig. 44, our results are compared with the previous measurements. For 




1 1.2 1.4 

Vs, GeV/c 2 



Fig. 43. The Ksir mass distribution. Points 
with error bars are data while the his- 
togram shows the fitted result for the 
spectrum expected in a model incorpo- 
rating only K*(892). The background 
is already subtracted. 
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Table XVIII. The branching fractions of various decay modes with an r\ meson. The upper limits 
are at the 90% C.L. 

Decay mode B 

K~riu T ,lQ- i 1.58 ±0.05 ±0.09 

7r _ 7rVv, 10~ 3 1.35 ±0.03 ±0.07 

K--k°tiu t , lQ- b 4.6 ±1.1 ±0.4 

K%it-r)v T ,lQ->' 4.4 ±0.7 ±0.3 

K*-^ T ,1Q- A 1.34 ±0.12 ±0.09 

K~ K° s r)v T , 10~ 6 <4.5 

A'g7r"7r°^ T ,10" 5 < 2.5 

K-rrqv T ,lQ- & < 3.0 

ir~rrqu T , 10~ 6 < 7.4 

(A-7rV T ) nonrcs , IP' 5 < 3.5 



Table XIX. Comparison of the branching fractions of three hadron decay modes from Belle and 
BaBar. 

Decay mode BaBar Belle 

-K-ir + Tv-is r , % 8.83 ±0.01 ±0.13 8.42 ± 0.00to.25 

K-n + 7r-u T , % 0.273 ± 0.002 ± 0.009 0.330 ± O.OOl+g;"? 

K-K+-r-v t , % 0.1346 ±0.0010 ± 0.0036 0.155 ± O.OOl^'oos 

K-K + K-v T , IP' 5 1.58 ± 0.13 ± 0.12 3.29 ± 0.17^;^ 



all modes studied, the precision of the branching fraction measurements for both 
BaBar 1 ") and Belle is significantly better than previous results. The accuracy of 
our results is comparable to that of BaBar, but the central values show striking 
differences in all channels other than 7r~ 7r + 7r~ v t . For this mode, our result is lAa 
lower than that of BaBar, while for the other modes the branching fractions obtained 
by Belle are higher by 3.0a, 3.0a, and 5.4a than those of BaBar for the K~ir + it~v T , 
K~ K + ir~v T , and K~ K + K~v T modes, respectively. 

6.2.7. Summary of precise measurements 

These measurements, as well as additional measurements of missing modes, are 
very important for obtaining separately the inclusive branching fractions of vector, 
axial-vector, and strange decay modes and corresponding spectral functions. 

For the rare decay modes with branching fractions of less than 10~ 2 , there is a 
significant improvement compared to the previous experiments. 

§7. D° mixing and CPV 

The neutral D meson system is one of the four flavored neutral particle-antiparticle 
systems that can exhibit oscillations. Particle-antiparticle mixing causes an initial 
(at time t = 0) pure D° meson state to evolve in time to a linear combination of D° 
and D° states: 

\D°(t)) = 



\ D o )cosh ( i ^+l ri 



+ ^|DVnh^7prt 



x e 



(7-1) 
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Fig. 44. Summary of the branching fraction measurements for three-prong r decays. 



where the two parameters that describe the D° - D° mixing, x and y, are defined as 
the mass and width difference of the two mass eigenstates \D\^) = p\D°) ± q\D°): 

mi -™ 2 r 1 -r 2 a + r 2 

X = — — y = ^2I^' F = ^— ' (7 ' 2) 

and r is the mean decay width. The coefficients p and q are complex, satisfying the 
normalization condition \p\ 2 + \q\ 2 = 1. The time-dependent decay rates for D° — > f 
(favored) and D° — > f (suppressed) decays are given by: 

r(D°(t) -> /) = \A f \ 2 e~ rt P + ^ /|2 cosh(yrt) - K[A/] smh( y rt) 

+— cos ( xrt ) + 9 I A /] sm(xrt)^J , (7-3) 
1 + |A^| 2 

^ cosh(yrt) - ^[Xj 1 ] sinh(yrt) 



r(D\t) f) = \A 7 \ 2 







Q 


\- rt ( 






P 





|A- , 

^ cos(xrt) - ^[Xj 1 } sm(xrt) , (7-4) 

where Af = --J- and A-f = -^r-. The time evolution of neutral D meson decays is 
J p A/ J p Aj J 

exponential with lifetime t d o = modulated by the mixing parameters x and y 

(see the expressions above). Time-dependent measurements of D° and D° decays 

thus enable us to measure the mixing parameters x and y. Since the dependence on x 
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Fig. 45. Short distance (left) and long distance (right) contributions to D°-D° mixing in the SM. 

and y depends on the final state, different decay modes exhibit different sensitivities 
to the parameters x and y. 

Out of the four flavored neutral meson systems, the neutral D meson system 
is the only one in which down-type quarks are involved in the box diagram loop 
(see Fig. 45). The neutral pion is its own antiparticle and the top quark decays 
before it forms a hadron and therefore cannot oscillate. Hence studies of charm 
mixing offer a unique probe of NP via flavor changing neutral currents in the down- 
type quark sector. In the SM, mixing in the neutral D meson system can proceed 
through a double weak boson exchange (short distance contributions) represented 
by box diagrams, or through intermediate states that are accessible to both D° and 
D° (long distance effects), as represented in Fig. 45. The potentially large long 
distance contributions are non-perturbative and therefore difficult to estimate, so 
the predictions for the mixing parameters x and y within the SM span several orders 
of magnitude between 10 -8 and io -2 . 208 )' 209 ) Due to large uncertainties in the SM 
mixing predictions, it is difficult to identify NP contributions (a clear hint would be 
if x is found to be much larger than y); however, measurements can still provide 
useful and competitive constraints on many NP models, as will be discussed later. 

The study of CP violation in decays of charmed hadrons also holds the poten- 
tial for uncovering NP. In the SM, direct CP violation can occur in singly Cabibbo 
suppressed (SCS; c — > dud, c — > sus) decays, but not in Cabibbo favored (CF; 
c —7- sud) or doubly Cabibbo suppressed (DCS; c — > dus) decays. This is due to 
the fact that the final state particles in SCS decays contain at least one quark- 
antiquark pair of the same flavor, which makes a contribution from penguin-type 
or box amplitudes induced by virtual 6-quarks possible in addition to the tree am- 
plitudes. However, the contribution of these second order amplitudes is strongly 
suppressed by the small combination of CKM matrix elements V c i,V* b . Therefore, 
in processes involving charmed hadrons, mainly the first two generations of quarks 
are involved. From the Wolfenstein parameterization of the CKM matrix 13 ) one can 
see that the elements related to the first two generations of the quarks are nearly 
real. Of course, by using the unitarity of the matrix one can still parameterize and 
estimate the imaginary part of those elements. For example, examining the phase 
difference between the decays D°(D°) — > K + K~, one finds that it is 2 arg(V^ s T^* s ), 
which can be expressed using the unitarity and the Wolfenstein parameterization 
as 2A 2 \ 4 r] « 1CT 3 . Hence the expected CPV asymmetries in the charm sector are 
of the order of 10 -3 , which is small compared to the measured CP asymmetries in 
the bottom sector. Recently, with the experimental precision reaching the per mille 
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level, some authors 210 - 1 have argued that the asymmetries could be much larger than 
naively expected. Nevertheless, at the current level of experimental sensitivity, the 
measurements of the CPV in the charm sector are mainly a search for a significant 
effect, which would point to so-far unknown NP processes. 

7.1. Experimental techniques of D°-D° mixing and CP symmetry violation 

Often, the flavor of initially produced neutral D mesons needs to be tagged 
in order to identify D°-D° transitions and CP violation. The flavor is tagged by 
requiring that neutral D mesons originate from D* + — > D°ir + decays, where the 
charge of the pion accompanying D tags the flavor of the neutral D meson at 
production. Another common property of the measurements described below is 
the selection of D meson candidates based on the CM momentum, typically p* > 
2.5 GeV/c for data taken at the T(4S) resonance. This requirement completely 
removes charmed mesons arising from possibly CP- violating B meson decays that 
have a displaced production vertex. Hence the Belle charm samples consist entirely 
of e + e~ — > cc continuum data. 

The most precise constraints on the mixing parameters x and y are obtained 
using the time dependence of D° decays. In time-dependent measurements, the D° 
decay time is determined according to t = rri£)o(L-p D o) /\pd°\ 2 > where L is the vector 
joining the D 0, s production and decay vertices, and p D o and m D o are its momentum 
and nominal mass. The reconstructed tracks of D° decay products are refitted to 
a common vertex to determine the D° decay point, and then the D 0, s production 
point is determined from the kinematic fit of the D momentum vector with the 
beam spot profile. The decay-time uncertainty at is evaluated event-by-event from 
the covariance matrices of the production and decay vertices. Typically, for decays 
with two charged tracks in the final state, (at) ~ t d o/2. Candidates with badly 
reconstructed decay time (with large at) are excluded from the analysis. 

The mixing parameters are extracted by performing a fit to the decay-time 
distribution using the following PDF: 

/+oo 
r sig (t'; x, y)R sig (t - t')dt' + P hkg (t), (7-5) 
-oo 

where the signal contribution is a convolution of the (final state dependent) time- 
dependent decay rate (-T S ig) and the detector resolution function (i? s ig)- To reduce the 
systematic uncertainties related to the parameterization of the resolution function, 
kinematically similar decays (from high statistics control samples) are usually used 
to determine the resolution function parameters directly from data. The background 
Vbkgit) is parameterized using an exponential (to describe the background candidates 
originating from mis-reconstructed charm decays) and a 5 function (to describe ran- 
dom combinations of final state particles), each convolved with its corresponding 
resolution function. The parameters of the background PDF are determined using 
events populating the sideband region in the invariant mass of D° candidates. 

Experimental determinations of CPV can be divided into those using the decay 
time distribution of certain decays to determine the unknown parameters and those 
using the decay time- integrated methods. The unknown CPV parameters often 
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follow from the parameterization below: 

sin</> . (7-6) 



A f 



A 



2 2 



P 



1^1 

pA f 



A^ D / is the asymmetry from CPV in decays to a final state /, Am ^ is from 
CPV in the mixing, and sin^> / is a manifestation of CPV in the interference 
between decays with and without mixing. 

While in charged D meson processes only the CPV in decays is present, neutral 
charmed mesons may include contributions from all three types of violation. 

All measurements are performed blindly, i.e. the selection criteria are determined 
using samples of simulated events or data events that are statistically independent 
from those used to perform the measurement, in order to avoid possible biases. 

7.2. Time- dependent measurements of D°-D° mixing and CP violation 
7.2.1. Decays to CP eigenstates 

Belle found the first evidence for D°-D° mixing 211 ) in a data sample of 540 fb _1 
using the ratios of lifetimes extracted from a sample of D° mesons produced through 
the process D* + — > D°tt + , which decay to K~ir + , K~K + , or ir~ir + . The time- 
dependent decay rates of the CF mode, K~ir + , and the SCS modes h~h + (h = K or 
7r) are obtained from the time-dependent decay rates given in the previous section: 

r(D°(t) -> K-7r+,D°(t) -> K + tt~) oc e-^D" (7-7) 
r(D°{t),D°(t) -> h + hT) oc e - {1+ycp ^ T D^ (7-8) 

where we assume that x,y < 1 and \A f /A f \ = 1 (\A f /A f \ < 1) for D° meson 
decays to h~h + (K~ir + ). For D° — > h + h~ decays, linear terms in xt and yt are the 
first terms in the Taylor expansion of the exponential function above. The lifetime 
difference between decays to the CP eigenstates h~h + and CP-mixed state K~ir + , 
VCPi is defined as 

VCP = TyT7r± - 1 = y cos 4> - \a m x sin <j). (7-9) 
i~h+h- 2 

The lifetimes tk-k and Thh are the effective lifetimes extracted from samples of D° 
mesons decaying to the CP mixed final state K~ir + , and CP even final states 
K~K + and ■k~it + . If \q/p\ = 1 and 4> = aig(q/p) = 0(ir), CP symmetry in mixing 
and interference between mixing and decay is conserved, and hence the parameter 
yep corresponds to the mixing parameter y. In these time-dependent measurements 
of neutral D mesons decaying to CP eigenstates, indirect CP violation is also probed 
by comparing the lifetimes of D° and D° mesons decaying to CP eigenstates: 

D° D° 

A r = T ff " ~ ^ ~ = -A M y cos </> - x sin cf>. (7-10) 
T h+h- + T h+h- 

By performing a simultaneous fit to the decay-time distributions of around 0.15 (1.2) 
million reconstructed tagged D° decays with purity above 90% to h~h + (K~ir + ), 
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Fig. 46. Results of the simultaneous fit to decay-time distributions of (a) D° — s- K + K , (b) 
D° -> R--K+, and (c) D° — > tv + tt decays. The cross-hatched area represents background 
contributions, the shape of which was fitted using D invariant mass sideband events, (d) Ratio 
of decay-time distributions between D° — > K + K~ , -k + -k~ , and D° — > K~tt + decays. The solid 
line is a fit to the data points. 



Belle found y C p = (1.13 ± 0.32 ± 0.25)% and A r = (0.01 ± 0.30 ± 0.15)%. Figure 46 
shows the decay-time distributions with fit results superimposed as well as the decay- 
time dependent ratio of D° decays to CP-even eigenstates K~K + and tt~it + to the 
CP mixed final state K~tt + , as measured by Belle. 211 - > In case of yep = 0, this ratio 
would be constant, which is inconsistent with Belle's data at 3.2cr. No evidence for 
indirect CP violation is found. 

7.2.2. Decays to hadronic wrong sign decays 

Belle also performed a search for neutral D meson mixing and CP violation in 
a time-dependent study of DCS D° — > K + tt~ decays 212 - 1 based on 400 fb _1 of data. 
These decays (also referred to as wrong sign decays) can proceed both through 
mixing followed by a CF decay, D° — > D° — > K + ir~, or directly through a DCS 
decay such as D° — > K + it~ . To distinguish the two processes, an analysis of the 
decay time distribution is performed. The most general form (e.g. allowing for direct 
CP violation in DCS decays, mixing and interference between mixing and decay) for 
the time-dependent decay rates of the two-body wrong sign decays D° — > K + tt~ or 
D° — > K~n + to second order in x and y is given by: 
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Fig. 47. (Left) The decay-time distribution for wrong sign (WS) events. Superimposed on the data 
(points with error bars) are projections of the decay-time fit when no CPV is assumed. The 
mixing and interference terms are shown at the 95% C.L. upper limits. (Right) 95% C.L. regions 
for (x' 2 ,y). The point is the best fit result assuming CP conservation. The dotted (dashed) 
curve is the statistical (statistical and systematic) contour for no CPV . The solid curve is the 
statistical and systematic contour in the CPV-allowed case. 
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(7-11) 



where Rd is the ratio of DCS to CF decay rates, and the parameters x' and y' are 
rotated mixing parameters, which are rotated by an unknown strong phase difference 
between the DCS and CF amplitudes, 6k-k- x' = x cos Sk-k + ysinJ^Tr and y' = 
y cos 5 Kn — xs'm 5 Kit- The three terms in the time-dependent decay rates of wrong 
sign decays are due to the DCS amplitude, the interference between the DCS and 
CF amplitudes, and the CF amplitude, respectively. In addition to the wrong sign 
decays, the Cabibbo favored (or right sign) D° — > K~ir + decays are reconstructed 
in order to determine the resolution function parameters, as well as the distribution 
of wrong sign signal events in D° invariant mass and mass difference distributions, 
which are fitted to extract the number of correctly reconstructed wrong sign decays. 

From a fit to the decay-time distribution of around 4 x 10 3 signal wrong sign 
decays (and with purity around 50%) Belle found x' 2 = (O.I8I0J3) x 10~ 3 and 
y' = (0.6^3'g) x 10~ 3 assuming no CP violation (setting Ad = Am = 4> = in Eq. 
7-11). The errors in x' 2 and y' include both statistical and systematic uncertainties. 
A projection of this fit superimposed on the data is shown in Fig. 47 and the non- 
mixing point (x' 2 = y' = 0) is found to be excluded at 95% C.L. In a second fit, 
CP- violating parameters are allowed to float and no evidence for either direct or 
indirect CPV is found. Belle obtains the following 95% C.L. intervals for CP- 
violating parameters: A D £ (-76, 107) x 10~ 3 and A M £ (-995, 1000) x 10" 3 . 
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7.2.3. Self-conjugated three-body decays 

Several intermediate resonances can contribute to a hadronic three-body de- 
cay of a neutral D meson. For example, D° — > Kg-K + Tr~ decays can proceed via 
D° ->• K*-tt+ (CF amplitude), D° ->• K^p° (SCS amplitude and CP eigenstate), 
D° — > K* + tt~ (DCS amplitude), and others. In the isobar model, the instantaneous 
amplitudes for D° and -D decays to the three-body final state / are parameterized 
as a sum of Breit-Wigner resonances and a constant nonresonant term (in the case 
of no direct CP violation, e.g., there is no difference between amplitudes and phases 
in D° and D° decays): 

A f {s + ,S-) = Y,are i4 ' r A r (s+,s-)+a NR e i *™, (7-12) 

r 

A f {s + ,s-) = ^a r e i ^A(s-,s+) + a NR e i * NR , (7-13) 

r 

where y/s± is the invariant mass of a pair of final state particles (e.g. Kgir^, and 
the sum runs over possible intermediate resonances r. The time-dependent decay 
rate for D° decays is thus given by (the corresponding expression for D° decays is 
obtained by multiplying the equation below by \p/q\ 2 )- 

^^>aH I (^,..)|VJf«) + | A (^..)|» e -JC-) 



+29l[Ai(s + , s_)^2(«+, cos 



t 

x— 

T 



+2Q[^i(s+,s_)^(s + ,s_)]sin ( x- )e~r, (7-14) 



T 



where Ai t 2{s+, s~) = \ (^Af(s+, s_) ± |^l/(s+, S-)^ • Different regions in the s+ — s_ 

plane (also called the Dalitz plot) exhibit different forms of time dependence, as can 
be seen from the above decay rate; therefore, the time-dependent Dalitz plot analysis 
of neutral D meson decays to a self-conjugated three-body final state enables us to 
measure the x and y parameters simultaneously. In the case where the analysis is 
performed separately for D° and D° samples, indirect CP violation can be probed 
by measuring the amplitude and phase of q/p. This method of measuring the mixing 
parameters x and y was pioneered by CLEO in D° — > Kgir + ir~ decays, 213 ) and was 
applied by Belle to D° -»• K° s tt + tt- decays 214 -* using a data sample of 540 fb 1 in 
which around 530 000 signal events are reconstructed with a purity of around 95%. 

The decay amplitude is not a priori known and has to be extracted from the 
data. This is done by first performing a time-integrated Dalitz plot analysis in 
which a model for the decay amplitude (A(s+, s_)) that describes the observed de- 
cay kinematics best is obtained. Belle finds that a good description of the Dalitz 
plot is obtained when 18 quasi-two-body resonances and a nonresonant term are 
used (see Fig. 48 for the results of the fit). Once the decay amplitude compo- 
sition is determined, a time-dependent Dalitz analysis is performed to determine 
the mixing parameters. In a fit with conserved CP symmetry (\q/p\ = 1 and 
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Fig. 48. Projections of the Dalitz distribution (points with error bars) and the fit result (curve) 
for D° — > KgTr + ir~ decays. 214 ^ Here, m\ corresponds to m 2 (KgTT ) for D° decays and to 
m 2 {K° S TT*) for D° decays. 

= 0) the mixing parameters are found to be: x = (0.80 ± 0.29±g;g?±g;J2)% and 
y = (0.33 ± 0.24l£<gl^)%, excluding the non-mixing point with 95% C.L. The 
errors quoted are the statistical, systematic error arising from experimental sources 
(e.g. modeling of background, resolution function, etc.) and the systematic error 
arising from the decay model (determined by using alternative models with different 
parameterizations, excluding resonances with small contributions, etc.). In a fit al- 
lowing for CPV, the \q/p\ and 4> parameters are found to be consistent with no CP 
violation: \q/p\ = 0.86l°J°l°;°6 ± 0.08 and <j> = (-14^1^)°. 

Large fractions of D° — > K®K + K~ decays proceed via D° — > Kg(f> (CP-odd) 
and D° — > Kgao(980) (CP-even) decays. Belle took advantage of this fact and 
performed a measurement of the yep mixing parameter integrated over the Dalitz 
plot using an untagged sample of D° — > KgK + K~ decays. 215 ^ We measure the 
effective lifetimes of D° mesons, ton, OFF) in two different regions of K + K~ invariant 
mass (at the <p peak (ON) and in <p sidebands (OFF)), which are given by ton,off = 
(1 + (1 — 2/oN,OFF)yc7p) r D°) where /on,off is the CP-even fraction in the ON or 
OFF region calculated using the decay model obtained by BaBar. 216 ) The obtained 
value of yep = (0.11 ± 0.61 ± 0.52)% is consistent with yep obtained in D° — > hh 
decays. 

7.3. World average and constraints on new physics models 

Various measurements of D°-D° mixing performed in different decay modes can 
be combined to obtain the world average values of x and y. The Charm subgroup 
of the Heavy Flavor Averaging Group has done this by performing a global x 2 fit 
from measurements of relevant observables 28 ^ performed by the Belle*) , BaBar, CDF, 
LHCb, CLEO-c, Focus, and FNAL E791 experiments. The world average values are 
found to be 

x = (0.63 ±lg$)%, (7-15) 
y = (0.75 ± 0.12)%, (7-16) 

*' In the world average fit, HFAG also includes Belle's time-integrated measurement of the 
mixing rate Rm = (x 2 + y 2 )/ in D° — > K + £ue 217 ^ which is not described in detail in this report. 
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and 

+0.18 



(0.88 i^), (7-17) 
(-10.1±ljj;{j)°. (7-18) 



The non-mixing point, (x,y) = (0,0), is excluded at the 10.2 standard deviation 
level while the \q/p\ and 4> values are consistent with conservation of CP symmetry 
in mixing and interference between mixing and decay. 

Golowich et al. 218 -* studied the implications of existing D°-D° measurements on 
many NP models. In many scenarios they found strong constraints that surpass those 
from other search techniques and provide an important test of flavor changing neutral 
currents in the up-quark sector. One simple extension to the SM that they studied 
is the addition of a fourth family of fermions. The obtained constraint on the CKM 
mixing parameters V c b>V* b , (b' is the down-quark of the fourth generation) is an order 
of magnitude more restrictive than those obtained from unitarity considerations of 
the CKM matrix. 

7.4. Time-integrated measurements of CPV in charm 

In the time-integrated measurements one usually determines the asymmetry of 
the partial decay widths, 



(/ _ r(D^f)-r(D^f) 
r(D^f) + r(D^f) 



A t cp = - u ;^ n y_ (7 . 19) 



The measured asymmetry 

Af _ N(D+ -> f) - N(D- -> f) 

^ N ( D + /) + N(D~ -> /) ' ( ' 

where N denotes the number of detected decays, receives a contribution from sev- 
eral non-CP violating sources, the detector-induced asymmetries due to a possible 
asymmetry in the acceptance of positively and negatively charged pions and kaons, 
or the different acceptances for neutral kaons and their antiparticles. In addition, the 
physical forward-backward asymmetry in the process e + e~ — > cc affects the mea- 
sured asymmetry, as we will see in the following. All these effects must be carefully 
determined using control data samples in order to achieve an accuracy at the level 
of 0(1O~ 3 ), the level of uncertainty of Aqp measurements in various final states 
reached by the Belle experiment. The existing MC simulation tools cannot be used 
for corrections at this level of accuracy. 

Currently the best sensitivity on Aqp at Belle has been achieved in the decays 
D+ -»• Tr+ifO. This decay mode is a mixture of CF (D+ -> tt + K°) and DCS 
(D + — > tt + K°) decay. If NP processes with unknown CP- violating phases would 
contribute, the CPV in the decays may be significantly different from zero. The 
measured asymmetry in these decays can be written as 



(7-21) 
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where A^f? is the physical CPV asymmetry, the detector- induced asymmetry 
between the tt + and ir~ reconstruction efficiencies, and Afb the contribution of the 
forward-backward asymmetry. The latter is an odd function of the D meson polar 
angle in the CM cos#* (see e.g. Ref. 219 )), while the first term is independent of any 
kinematic variables. The detector-induced asymmetry depends on the momentum 
and the polar angle of the charged track in the laboratory frame. In bins of these 
variables the measured asymmetry can be corrected for A^ + using samples of D° — > 
K~it + it and D + — > K~tt + it + decays. The measured asymmetries for these decays 
are 

Ag* = A FB + Af ~ + A^ + A$ 

A*™ = A FB + Af " + J$ , (7-22) 

assuming negligible CP violation in the Cabibbo favored D meson decays and the 
universality of the forward-backward asymmetry for different types of charmed 
mesons *\ By inspecting Eqs. (7-22) one finds that in the difference of the measured 
asymmetries in D° — > K~ir + ir° and D + — > K~it + it + some of the detector-induced 
asymmetries and the forward-backward contribution cancel and hence one can de- 
termine Al + . In turn, A^ + is then used to correct A^ 7r+ in bins of the charged pion 
momentum and polar angle, and to extract A^, p n . 22 °) However, in the D + — > h + K® 
decay modes, one needs additional corrections due to the presence of a neutral kaon 
in the final state. In such D + meson decay modes either a K° or K° is produced, 
which interact differently in the detector material. However, in the final state a Kg 
is reconstructed, and hence this affects the value of the asymmetry. A separate ded- 
icated study 221 ) was performed and the appropriate correction factor applied to the 
asymmetry. Furthermore, because of the CP violation in the neutral kaon system, 
the asymmetry expected in this final state with K s is A Ks = (—0.332 ± 0.006)%. 
The Belle result is given in Table XX and is in good agreement with the expectation 
due to CP violation in the neutral kaon system. 

Belle searched extensively for non-zero time- integrated CP asymmetries in a 
number of other decay modes and achieved the best sensitivity in many of these. 
The results (see Table XX) are consistent with no CPV at levels varying from 
0(Kr 2 ) to 0(KT 3 ). 

7.5. Conclusions 

With the world's largest sample of recorded charmed hadron decays Belle has 
experimentally observed mixing phenomena in the last remaining neutral meson 
system, D°. The mixing parameters in this system are nowadays becoming a pre- 
cision measurement, with world average values 28 - 1 of x = (0.63 ± l[j 20)% anc ^ ^ = 
(0.75 ± 0.12)%. Further measurements and advances in theoretical predictions are 
required to determine whether the observed values are consistent with the SM or re- 
ceive contributions from NP. Furthermore, an extensive search for CPV in the charm 

*-* Within the SM only Cabibbo suppressed decays of charmed mesons have two possible am- 
plitudes with different weak and strong phases — the tree and the penguin amplitude — which is a 
necessary condition for non-zero CPV in decays. 
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Table XX. Measured time-integrated CPV asymmetries in the D meson system. 



Decay mode 


C (fb" 1 ) 


A C p (%) 


Comment 


Ref. 


D u - 






-0.28 ±0.19 ±0.10 






D° - 


>K° s ri 


791 


±0.54 ±0.51 ±0.16 




222) 


D°- 


>K%rf 




±0.98 ±0.67 ±0.14 






D°- 




540 


±0.43 ±0.52 ±0.12 




223) 


D°- 


> K+K~ 




-0.43 ±0.30 ±0.11 






D°- 




532 


±0.43 ± 1.30 




224) 




+ K+n-n" 


281 


-0.6 ±5.3 




225) 


D°- 






-1.8 ±4.4 






D+ - 




977 


-0.363 ± 0.094 ± 0.067 


signif. asymmetry due to K% 


220) 


D+ - 


->■ + 


955 


±0.51 ±0.28 ±0.05 


universality of Afb in Df 
and D + decays to v + 4> tested 


226) 


D+ - 


rjiT + 


791 


±1.74 ±1.13 ±0.19 


D+ -+ K+ V v> 


227) 


D+ - 


i + 

-> 7] 7T 


791 


-0.12 ±1.12 ±0.17 


also observed 




D+ - 


-^K%K+ 


673 


-0.16 ±0.58 ±0.25 




228) 


Df- 




673 


±5.45 ±2.50 ±0.33 




228) 


Df- 






±0.12 ±0.36 ±0.22 




228) 



sector was carried out. The measurement methods that were developed allowed for 
the observation of a significant CPV asymmetry in decay modes with a neutral kaon 
in the final state and sensitivities to possible time- integrated CP asymmetries at the 
per mille level in a variety of decay modes. No significant indirect CP violation has 
been observed so far. 

§8. B physics at the Y(5S) 

The T(10860), (M = 10876 ± 11 MeV/c 2 , r = 55 ± 28 MeV), 15 ) is generally 
interpreted as the T(5S), the fourth excitation of the vector bound state of bb, and is 
just above B*B* threshold. The Belle experiment collected a total of 121.4 fb -1 at 
the T(10860) peak energy and a total of 27.6 fb -1 at off-peak CM energies nearby, 
between 10.683 and 11.021 GeV. The on-resonance data sample corresponds to 37 
million resonance events and includes 7.1 million B s events. These data were ana- 
lyzed to pursue investigations of B s meson properties, hadronization to B q and B s 
events (q is a u- or <i-quark), energy dependence of various types of events, and pos- 
sible new bottomonia and bottomonium-like states. Published on-peak results are 
based on two subsets, 1.86 fb -1 and 23.6 fb -1 (including the 1.86 fb^ 1 ), as well as 
the full set of 121.4 fb -1 , which will be referred to as sets 2FB, 24FB, and 121FB, 
respectively. 

The e + e _ — > T(10860) is an excellent venue for studying several aspects of 
B s decay; given clean, efficiently triggered events with precisely known CM energy, 
collected by a well-understood detector, the Belle experiment has been uniquely po- 
sitioned to measure absolute branching fractions, access modes that include photons 
in the final state, and do comparative studies of B and B s mesons with minimal 
systematic uncertainties. 
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8.1. Bs masses: method of full reconstruction 

At the energy of the T(10860), three types of B s events are allowed: B S B S , B*B*, 

and B S B* (and B S B*) events. Each is an exclusive 2-body decay, so the energy of 

(*) 

the daughter Bg ' in the collision CM frame is fully constrained. The method of 
"full reconstruction," where all decay products are detected and measured, was used 
with great success for B q at the T(4S). The reconstruction of B s in B S B S events 
is analogous: each B s carries energy equal to the beam energy (in the collision 
CM system), so upon reconstructing a candidate, the quantity AE accumulates 
at AE = GeV and M^ c at the true B s mass, m^. In the decay B* — > B s j, 
the photon carries away essentially all of the released energy, which is equal to the 
mass difference, 5M ps 50 MeV/c 2 . In a B*B S event, the B s (B*) carries energy 
~ ^bcam - SMc 2 /2 (~ £ b eam + SMc 2 /2). The daughter B s from B* ->• B s j carries 
energy ~ E\ 2lCSLm — 5Mc 2 /2. Thus, for both of these B s 7 s, one can expect reconstructed 
decays to accumulate around AE = —SMc 2 /2 and M^ c = m,B s +SM/2. Carrying the 
process another step further, both i? s 's in B*B* events accumulate at AE = —5Mc 2 
and Mb c = ms s + 6M = ms*- Given the Belle detector's momentum resolution, 
these three event types accumulate in well-separated regions of AE and M^ c , as 
shown in Fig. 49 for B s — > D~ir + candidates, signal MC simulations, and data. 229 -* 



J 5.42 - 
0) 

a 5.4 - 

^5.38 r 
5.36 r 
5.34 r 
5.32 r 
5*JT 



0.05 

A E (GeV) 



> 5.42 p.- 
O 5.4 ^ 
^5.38 re 
5.36 j- ' 
5.34 H' 



0.05 

A E (GeV) 



Fig. 49. Illustration of the full reconstruction method, B s — > D s -k. Distributions in AE and Mbc of 
candidates, (left) Monte Carlo simulation, (right) data, 24FB set. Also shown are signal regions 
for B*B* (upper signal box), B'*B a (middle box), and B s B a (lower box) events. 229 ' 

As can be seen from Fig. 49, the B*B* events dominate in the data. As explained 
above, the M^ c distribution peaks at mg s (with very minor corrections). The B* — B s 
mass difference is found from the mean AE of the candidates; the B s candidate 
mass reconstructed as M£ c = y/{E£ clim + (AE)) 2 — (j>* and ) 2 accumulates at the B s 
mass. The modes B s -»• D+^-{D s -»• K+R-)^, K*°(^ K + K-)R-, K s (^ 
7r + ir~)K~} were reconstructed for this measurement. B*B* candidates are selected 
by requiring —0.08 < AE < —0.02 GeV. From the 24FB data set we measure 229 ) 

m Bt = 5416.4 ± 0.4 ± 0.5 MeV/c 2 
m Bs = 5364.4 ± 1.3 ± 0.7 MeV/c 2 

8.2. Event composition at the T(10860) peak 

To study production and decay rates of B s , their abundance and properties in 
X(10860) events are needed. This evaluation proceeds in three steps. First, we 
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measure the hadronic 66 cross section. 230 ) We then find the fraction of 66 events 
containing B s . 230 ^ Finally, we measure the relative rates to the three possible event 
types. 229 ) 

8.2.1. a(e + e~ -»• bb) 

As is the case at the T(4S), bb events at the T(10860) (where "66" includes 
both resonance and 66 continuum events, which are indistinguishable) are readily 
distinguished statistically from the continuum of lighter quarks e + e~ — >■ qq (q = 
u, d, s, c) via their distribution in R2, the ratio of the second and zeroth Fox- Wolfram 
moments, 25 ^ a measure of "jettiness" that tends to be lower for the more isotropic 
66 events. The R2 distribution for the 2FB data and a scaled continuum sample 
are shown in Fig. 50. We found o e+e -^ bl = (3.01 ± 0.02 ± 0.16) x 10 2 pb, 230 ) 
which constitutes ~10% of the total hadronic cross section. This value was averaged 
with the corresponding result from CLEO 231 ) to obtain the PDG average of 07, = 
(3.02 ±0.14) x 10 2 pb. 15 ) 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Fig. 50. Distribution in R2, (histogram) data 
set 2FB, at the T(10860) and (points) con- 
tinuum below T(4S), scaled. 




Fig. 51. Distribution of D s in x = Pd s / ■sj E^ eajn — m\ 
(points) T(10860) and (histogram) scaled 
continuum. 230 ' 



8.2.2. a{e + e- -> B s B s )/a{e + e~ -»• 66) 

The fraction (f s ) of 66 events that hadronize to B s (B^B^) may be determined 
through measurement of the inclusive rate fi(T(10860) — > D S X) = Br- The B s 
decays predominantly via the spectator mechanism, as do the lighter B mesons, and 
as such we can assume a direct correspondence between B — > DX and B s — > D S X 
for a large fraction of the decays. Based on our understanding of the mechanisms of 
B decay and the measured branching fractions B{B — > DX) and B{B — > D S X), a 
reasonable estimate may be made: 232 ) B(B S -»• D S X) = (92 ± 11)%. The inclusive 
rate of T(5S)— > D S X is an average over B s , Bj, and B u , weighted by abundance: 

B ( r ( ss)->n,x ) =hB(B ^ DtX) 
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Hr-f,) Bini V D - X) <<»> 

where f s is the fraction of B s and we assume that B^ and B u are produced equally 
and that non-B production is negligible. The distributions of D s in normalized mo- 
mentum x = Pd 3 /\J -E'beam ~~ m L> s f° r Y(5S) and scaled continuum data are shown 
in Fig. 51. The measured value, B(T(5S) -»• D s X)/2 = (22.6 ± 1.2 ± 2.8)% for the 
2FB data set, is fed into Eq. (8-1) and solved to obtain 230 ) f s = (16.4 ± 1.4 ± 4.1)%, 
which corresponds to (4.95 ± 1.31) x 10 4 B s events/fb _1 . The same analysis may 
be performed for D° to obtain an independent value of / s , albeit with larger uncer- 
tainties; B(B S ->• D°X) < ^(.Bq ->• The results are combined to obtain 230 ) 
/ s = (18.0 ± 1.3 ± 3.2)%. The Belle result is averaged with the corresponding CLEO 
result 231 ) to obtain the PDG average 15 ) 

f 9 = (i9.5t 3 : °)%. 

The same method applied to the 121FB set yields 

f 8 = (17.1 ± 3.0)%. 

8.2.3. B*B* : B*B S : B S B S 

As described in Sect. 8.1, reconstructed B s signals from the three event types 
are well separated in AE and M^ c . These three modes account for 100% of B s 
events, so the fraction comprised by each is derived from a simultaneous fit to AE 
and Mb c that yields all three signals. For this measurement we use B s — > D~ir + , 
the mode with the greatest statistical significance. To date, statistically significant 
signals have been observed in the B*B* and B*B S + B S B* channels in the 24FB 
data set, from which we obtain 229 ) 

_ q(e+e- -> B*B*) mn i+3.8, n9W 
fB t B t = ^T^ W) = (90.1_, ±0.2)% 

_ a(e+e-^B* 3 B a + B a B* 3j , :! 

fB'.B, = — w w = (7.3_ 3 ± 0.1)% 

<r(e+e- -)• B s 'B s ') 



The value f B *B* = (87.0 ± 1.7)% from the 121FB data set (unpublished) 233 ) is used 
in evaluating branching fractions from the 121FB set. 

8.2.4. B(*)BW(7r)(7r) 

The well-tuned methods of B reconstruction at the T(4S) (Sect. 4.4.2) have been 
applied to study the more complicated assortment of B events at the T(108 60). 234 ) 
The following final states that include non-strange B mesons are energetically al- 
lowed: B { * ] B { *\ B q B ( *\, B q B qirn. The relative rates can improve our understand- 
ing of hadronization dynamics. Neutral and charged -B's are reconstructed in the 
following modes and submodes: B + ->■ J/ipK + , D°ir + ; B° -)• J/ipK*°, D~ir + ; 
J/ip -»• e+e", /t/V - ; K*° -»• K + tt~; D° -»• D~ -»• ^+^~7r~. 

As with the fully reconstructed B s , the signal events populate the (AE,M^ C ) plane 
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in clusters depending on the type of event. Figure 52(a) shows the projections in Mb c 
of the distributions for the various event types. The distribution of candidates in 
data, after background subtraction, are shown in Fig. 52(b). While the distributions 
for events containing additional pions overlap each other, it is clear from the data 
that their contribution is relatively small and that the majority of the rate is due to 
two-body events, B^B^*\ It is also noted that there is an accumulation of events 
in the region of high M^ c , where BBtttt events would accumulate, according to the 
MC simulation. The fractions of bb events fragmenting to BB, B*B, and B*B* are 
measured to be (5.5±J;g ± 0.4)%, (13.7 ± 1.3 ± 1.1)%, and (37.51^ ± 3.0)%, respec- 
tively. The events where M^ c is above the two-body limit are grouped together as 
"large M bc " and found to comprise (17.5±i;| ± 1.3)%. 




M bc (GeV/c z ) 



(b) 



5.35 



5.4 



M bc (GeV/O 



Fig. 52. (a) Distributions in M hc for reconstructed B° -¥ D~n + , for BB, BB* + B*B, B* B* , 
and BB-KiT channels (cross-hatched histograms, left to right) and for the three-body channels 
BB* -k + B*B-k (plain histogram), BB-k (dotted), and B*B* -k (dashed). The distributions are 
normalized to unity, (b) Mbc distribution in data after background subtraction. The sum of the 
five studied B decays (points with error bars) and results of the fit (histogram) used to extract 
the two-body channel fractions are shown. 



Events containing one or more additional pions may be identified by pairing 
reconstructed -B's with additional charged pions in the event and examining the 
residual, or missing, event energy and momentum, E m [ ss and P m i ss , which by infer- 
ence are carried by the opposing B^ and up to one additional pion. From these 
we reconstruct AE m i SS and Mb c , m iss- Projections onto AE m i ss + Mb c , m iss — mB for 
various simulated event types are shown in Fig. 53(a). The corresponding distri- 
bution in data, with the fit result, is shown in Fig. 53(b). The fractions of bb 
events hadronizing to three-body modes BBtt, BB*tt, and B*B*ir are found to 
be (0.0 ± 1.2 ± 0.3)%, (7.3 +Jj ± 0.8)%, and (1.0 t\i ± 0.4)%, respectively. Para- 
doxically, no evidence for BBtttt was observed, so this channel does not account 
for the remaining (9.2 +28 ± 1-0)% of the "large Mb c " contribution observed in 
T(10860)— > BX. The residual is quantitatively consistent with initial state radiation, 
e + e~ — > e + e~7, e + e~ — > bb, where about half the bb form the T(4S) resonance. 234 ) 
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Fig. 53. (a) The AE m i BB + M bc ,miss — ms distribution normalized per reconstructed B meson for 
MC simulated B+ -»■ J/ipK+ decays in the (peaks from left to right) BBtt + , BB* Tv+ + B*Bn + , 
B* B* 7r + , and BBtttt channels, (b) The AEmiss + Mb c ,miss — tub data distribution for right- 
sign B~/°n + combinations for all five studied B modes. The curve shows the result of the 

fit 234 ' 



8.3. B s decays 

To a large degree, the general properties of the B s meson parallel those of the 
non-strange B mesons. Like its lighter cousins, the B s is expected to decay predomi- 
nantly by a spectator process, where the lighter valence quark has no role in the weak 
interaction, and its spectator-dominated properties such as lifetime are expected to 
be similar. This is serendipitous at the .B-factory, allowing many of the techniques 
developed for analysis of i?'s at the T(4S) to be applied to B s at the T(5S). Further- 
more, close correspondences between the hadronic final states in spectator decays of 
B s and B^ allow for sensitive tests of quark-hadron duality and of hadronic models 
that may reduce theoretical uncertainties limiting precision CKM tests in B physics. 

Being electrically neutral, the B s experiences mixing and may thus address ques- 
tions of interest regarding CP violation and roles for physics beyond the Standard 
Model. Notably, B s experiences a much higher rate of mixing than Bj, and very 
little CP violation in the SM. 

All branching fractions described in this section are listed in Table XXI. The 
branching fractions measured in the 24FB set are evaluated using /b*b* = (90.1^4 q)%, 
f s = (19.5±^)%, and er e+e -^ b i = 0.302±0.014 nb (a weighted average from 230 )' 231 )). 
The results based on the 121FB set have used NgW^w = (7.1 ± 1.3) x 10 6 = 
£ x °e+e-^bb x /- and Jb*b* = (87.1 ± 3.0)%. 
8.3.1. Modes with single D s 

The decays B s — > D^~h + , where h is a light non-strange meson, proceed domi- 
nantly via a CKM-favored spectator process. The D s are reconstructed in the modes 
K + K-)ir~, K*°(^ K + K-)K~, and K s {-^ it + it-)K- and p ± in t^tt . As 
described in Sect. 8.1, the signal is extracted by fitting the distributions in AE and 
Mbc (and decay angles, for B s — > D* s ~ p + ). Shown in Fig. 54(left) is the projection 
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into Mbc for B s — > Djir + candidates in the 24FB set. The branching fraction for 
modes other than D~ir + are obtained using only the B*B* sample and the value of 
fB*B* measured with D~ ~ir + . 229 )> 235 ) 

The distribution of the angle between the B s momentum and the beam axis in 
the CM frame, 9* B , is of theoretical interest 236 ) and is presented in Fig. 54(right) 
for the signal events in the B*B* region. A fit of the distribution to 1 + a cos 2 0* Ba 
returns x 2 /n.d.f. = 8.74/8 with a = -0.59+°;}g. We naively expect a = -0.27 by 
summing over all the possible polarization states. 

For B s — > D*~p + , a pseudoscalar decay to two vectors, the distributions in the 
helicity angles D *- and 9 p + depend on the relative contribution from the different 
helicity states, which depends on the detailed hadronization mechanism for the de- 
cay; for example, the factorization hypothesis predicts that longitudinal polarization 
dominates: Jl ~ 88%. 237 ' A four-dimensional fit yields (7.4a) signal events 

and fr - 1 n^ +a08+0 ' 03 235 ) 
ana jl — i - uo -o.io-o.04- 



Table XXI. Branching fractions with statistical and systematic uncertainties. A third uncertainty, 
due to f 3 , is quoted where it is separated from other systematics. The data set analyzed is 
identified in the rightmost column. 



Mode 


23(1CT 4 ) 


Data set 


Single- D s modes 






D-7T + 


36.7l 3 3 ji 4 4 J±4.9 


24FB 




24± 5 4 ± 3 ± 4 


24FB 




85±il ± 11 ± 13 


24FB 


DTp + 


118^20 ± 17± 18 


24FB 


cess modes 






J/lpT) 


5.10 ± 0.50 ±0.251^9 


121FB 


J/tpr/f 
J/V/o(980) 
J/V/o(1370) 


3.71 ±0.61 ±0.18±£ : g 

1 1C+0.31 +0.15+0.26 
llo -0. 19 -0.17 -0.18 
n oa+0. 11+0. 03+0. 08 
u -°^-0. 14 -0.02 -0.05 


121FB 
121FB 
121FB 


D*+D* s - 


200 ± 30 ± 50 


121FB 


D*+D- + c.c. 


180 ± 20 ± 40 


121FB 


DfDj 


581 1 1 ± 13 


121FB 


£)(*) + £)(*)" gum 


430 ± 40 tfo ± 90 


121FB 


hh modes 






K + R- 


0.38±g;Jg ± 0.05 ± 0.05 


24FB 


K°K° 


< 0.66 (90% C.L.) 


24FB 


K~7V + + c.c. 


< 0.26 (90% C.L.) 


24FB 




< 0.12 (90% C.L.) 


24FB 


<t>1 
77 


+0.18+0.1LJ 
u - d '-0.15-().ll 

< 870 (90% C.L.) 


24FB 
24FB 



8.3.2. Flavor- neutral channels 

An interesting characteristic of the B s stems from the fact that it experiences 
an appreciable rate to the flavor-neutral combination cess, via a tree-level CKM- 
favored process. The massiveness of the participating quarks and proximity to mass 
threshold argue for the applicability of predictions at the limit m, b ^ — > oo with 
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Fig. 54. (Left) Mb c distribution of B„ — Y D 3 n + candidates with AE in the B*B* signal region 
[-80,-17] McV, 24FB data set. 229) The different fitted components are shown with dashed 
curves for the signal, dotted curves for the B s — > D*~-k + background, and dash-dotted curves for 
the continuum. (Right) Fitted distribution of the cosine of the angle between the B s momentum 
and the beam axis in the CM frame for the T(5S) — >• B^Bl signal. 



(m b — 2m c ) — > and N c (number of colors) —> oo, where the cess final states are 
CP-even and the D^DJ and D* s + D* s ~ modes (along with DfD~) saturate the 
width difference Ar^ p between the two CP-eigenstates. 238 - 1 This parameter equals 
Ar s / coS(f) s , where Ar s is the decay width difference between the mass eigenstates, 
and (f> s is the CP-violating phase in B s -B s mixing. 239 ) Thus the summed branching 
fraction B(B® — > Dg ) gives a constraint in the Ar s -(j) s parameter space. 

Both parameters can receive contributions from NP; see, e.g., 

Refe> 240),24i) AssuminK 
negligible CP violation (</> s ~0), the branching fraction is related to Z\P S via 

Ar s /r s = 2B/(l-B). (8-2) 

The quantity of interest, the summed branching fraction B = B{B S — > D), D s , 
is more easily measured in e + e~ — > T(5S) than at a hadron machine because the 
decay D* — > can be fully reconstructed. 

The final Belle result is based on the 121FB set. 242 ) It includes the first measure- 
ment of the fraction of longitudinal polarization (_/x) of B® — > D* + D*~ . The final 
states reconstructed consist of D+Dj, D*+D£+D*-D+ (= D*±DJ), and D* + D*-, 
where D* + D+j, Df 0vr+, K° S K+, K*°K + , <Pp + , K° S K* + , and K*°K* + , 
K Q S ^vr+vr-, K*°^-K+n~, K*+^K^ir + , (jj^K+R-, p+^vr+vr , and vr°^77. 243 ) 

Events containing candidates satisfying 5.25 GeV/c 2 < M bc < 5.45 GeV/c 2 and 
—0.15 GeV < AE < 0.10 GeV are selected. Approximately half of the selected 
events have multiple B® — > candidates. These typically arise from pho- 

tons produced via ir° — > 77 that are wrongly assigned as D* daughters. For these 
events we select the candidate that minimizes a x 2 constructed from the recon- 
structed Df and (if present) D* s + masses. 

Signal yields are measured by performing a two-dimensional unbinned maximum- 
likelihood fit to the M hc ~AE distributions. The combinatorial effects of analyzing 
multiple multi-body decays present a particular challenge in this analysis. The signal 
PDFs have three components: correctly reconstructed (CR) decays; "wrong combi- 
nation" (WC) decays in which a non-signal track or 7 is included in place of a true 
daughter track or 7; and "cross-feed" (CF) decays in which a D^D'f {D* + D*~) is 
reconstructed as a DfD~ (D+D~ or D* s ), or a DfD~ (D* D^) is reconstructed 
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as a DgDf or D* s + D*~ (D* + D*~). In the former case, the 7 from D* + -¥ D+7 
is lost and AE is shifted down by 100-150 MeV; this is called "CF-down." In the 
latter case, an extraneous 7 is included and AE is shifted up by a similar amount; 
this is called "CF-up." In both cases M bc remains almost unchanged. The small 
contributions from B S B S and B S B* events are fixed relative to B*B* according to 
our measurement on B® — > D~tt + decays. 233 ) The fitted signal yields from B*B* 
only are used to determine the branching fractions. 

The projections of the fit are shown in Fig. 55. The branching fraction for chan- 
nel % is calculated as B i = Y i /{e l MC ■ N g _ • f B ^, • 2), where Y i is the fitted CR 

yield, and £ l MC is the MC signal efficiency with intermediate branching fractions 15 ) 
included. The statistical significance is calculated as y ' —2 ln(£o/£max)> where Co 
and £ max are the values of the likelihood function when the signal yield Y i is fixed to 
zero and when it is floated, respectively. We include systematic uncertainties (dis- 
cussed below) in the significance by smearing the likelihood function by a Gaussian 
having a width equal to the total systematic error related to the signal yield. 




M bc (GeV/c 2 ) M bc (GeV/c 2 ) M bc (GeV/c 2 ) 



Fig. 55. M bc projections and fit results, (left) B° -> Df D~ , (center) B° — > D^Df, (right) B° ->■ 
D* + D*~ . The red dashed curves show CR+WC signal; the blue and purple solid curves show 
CF; the gray solid curves show background; and the black solid curves show the total. 

Inserting the total B from Table XXI into Eq. 8-2 gives 

Ar s /r s = 0.090 ± 0.009 ± 0.023 , (8-3) 

where the first error is statistical and the second is systematic. This result has 
precision similar to that of other recent measurements. 244 -*' 245 ) The central value is 
consistent with, but lower than, the theoretical prediction; 240 ) the difference may be 
due to the unknown CP-odd component in B®^D* + D*~, and contributions from 
three-body final states. With more data these unknowns can be measured. The 
former is estimated to be only 6% for analogous B°^D* + D*- decays, 246 ) but the 
latter can be significant: Ref. 247 ) calculates Ar(B s -> D^K^)/ T s = 0.064 ± 
0.047. This calculation predicts Ar s /T s from L>i* )+ D^" alone to be 0.102 ± 0.030, 
which agrees well with our result. 

To measure f L , we perform an unbinned ML fit to the helicity angles 9 1 and 9 2 , 
which are the angles between the daughter 7 momentum and the opposite of the B s 



Physics achievements from the Belle experiment 



83 



momentum in the D* + and D* rest frames, respectively. The angular distribution 
is (|^+| 2 + \A_\ 2 ) (cos 2 6>i + l) (cos 2 (9 2 + l) + |^4 | 2 4sin 2 (9j_ sin 2 (9 2 , where A + , A_, 
and A are the three polarization amplitudes in the helicity basis. The fraction 
f L = \A \ 2 /(\A \ 2 + \A + \ 2 + \A_\ 2 ). We obtain 242 ) 

f L = 0.06 l°i 8 7 ±0.03, (8-4) 

where the first error is statistical and the second is systematic. Reconstruction of B s 
decays to well-defined CP final states are of interest for studies of CP violation. In 
the SM, mixing- mediated CP violation occurs in neutral mesons due to the complex 
argument of the product of CKM matrix elements participating in the mixing "box 
diagram." For B s the relevant product is Kt 2 V^g, which is real, so no significant 
asymmetry is expected. Searches for CP asymmetry in decays of B s thus present 
an opportunity to reveal NP. Such measurements will require the reconstruction of 
a sizable sample of CP-defined final states. 

The decays B s — > J/ipn^ (CP = +1) proceed by the same process as B — > 
J/ipK°, so the branching fractions may be estimated based on the measured branch- 
ing fractions, 15 ) B{B° d -)• J/ipK°) = 8.71 x 10" 4 : B{B S -> J/iprj) ps 3.5 x 10~ 4 , and 
B{B S -> J/iprf) ps 4.9 x 10" 4 . The decays are reconstructed in the following modes: 
J/tj) — > e + e~, r\ — > 77, ir + ir~ir°; n' — > 7]-k + it~ , p°7. The signals are extracted 

via a 2-dimensional fit in AE and Mb c . Projections in M^ c are shown in Fig. 56. 




Fig. 56. Projections in M bc , based on the 121FB data set at T(10860): (left) B s -¥ J/i>v(lj), 
(right) B s — > J/ipr/(Tv + Tv~Tv°). Solid curves show projections of fit results. Backgrounds are 
represented by the blue dotted curves. Two small bumps around 5.37 and 5.39 GeV/c 2 are 
contributions from BgBg and B*B° production channels, due to the overlap of the AE signal 
regions. 

The same b — > ccs process can also produce the decay B® — > J/*0/o(98O), another 
promising channel for CP studies, with the clear advantage of being an all-charged 
final state with no angular analysis required because of the J p = + quantum 
numbers of the /o(980). The mode was reconstructed as B s — > J / iJjtt + n~ , {J/ip — > 
e + e~}, analyzing the 121FB set. The fit to data include the /o(980) and 
another resonance in the irn mass spectrum at ~ 1.4 GeV 2 , fx (Fig. 57). The fx 
mass, measured at 1.405±0.015±{j;{j$ GeV/c 2 , is consistent with that of the /o(1370). 
The nonresonant yield is consistent with zero. 
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Fig. 57. Pion pair mass distribution for B s — > J/ipir + -K candidates in the 121FB set, for 
-79.7 MeV < AE < -19.7 MeV. The solid line represents the total PDF. The dash-dotted 
curve represents the total background, the dashed curve shows other J/ip background, and the 
dotted curve the nonresonant component. 

We have also searched for the 2-body CP-eigenstate modes B s — > K + K~, 
K°K°, and 7r _ 7r + , as well as the flavored mode B s — > K~ir + , in the 24FB data 
set. 248 ) The findings for K+K~ and K°K° were the first absolute branching frac- 
tion and first reported limit, respectively. 

8.3.3. Radiative decays 

Radiative penguin decays, which produce a photon via a one-loop Feynman 
diagram, are a promising venue to search for physics beyond the SM because particles 
at mass scales not yet directly accessible at accelerators can contribute to such loop 
effects. The B s —> <frf mode is a radiative process described within the SM by a b — > 
s"f penguin diagram; it is the counterpart of the B — > K*(892)j decay. In the SM, the 
B s — > 07 branching fraction has been computed with 30 % uncertainty to be about 
40 x lO" 6 . 249 )' 250 ) This channel was first observed at Belle, with <f> reconstructed 
in the mode K + i^~ 251 ) . For photon selection, major sources of background in 
the signal region included ir° — > 77 and rj — > 77 as well as calorimeter hits that 
were out of time with the beam crossing. Based on the 24FB set, we reported 
B(B S — > (pj) = (57l]J(stat)l]; 2 (syst)) x 10~ 6 , which is in agreement with both 
the SM predictions and with extrapolations from measured B + — > i^*(892) + 7 and 
B° — > K*(892)°7 decay branching fractions. 251 ) 

8.3.4. Modes suppressed in the Standard Model 

The B s — > 77 mode is described in the SM by a penguin annihilation diagram 
(Fig. 58), and its branching fraction has been calculated to be in the range (0.5- 
1.0)xI0~ 6 . 252) - 254) Belle has searched for this mode in the 24FB data set. 251) No 
significant signal was observed (Fig. 59), and a 90%C.L. upper limit oiB{B s — > 77) < 
8.7 x 10 -6 was obtained. This limit significantly improves on the previously reported 
one and is only an order of magnitude larger than the SM prediction, providing the 
possibility of observing this decay at a future Super B-factory. 255 )' 256 ) 

8.4. Measurement o/sin2</>i 

The method of full B reconstruction, used to study the assortment of B events at 
the T(5S'), 234 ) has been applied to a novel tag to measure sin 20i. 25? ) Three-body fi- 
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nal states bM°{-> B (j)}B^-tt + (+c.c.) are identified through full reconstruction 
of a neutral B in a C-P-eigenstate and a charged pion. The event residue, consisting 
of a charged B and up to two photons, is characterized through "missing mass," 
calculated through energy and momentum conservation: 

£miss = £beam ~ E B o w ; p miss = -p B 0^, MM(B°Tt) = M miss = yj E^ - p 2 ^. 

The missing mass distributions are well separated for BBtttt, BBtt, BB*tt, and 
B*B*tt events, as can be seen in Fig. 60(left). The sign of the charged pion tags the 
initial flavor of the neutral B and enables a time-independent measurement of CP 
asymmetry, which is related to sin20i as: 

= N BBn - - N BBn+ = Sx + A 
BBn ~ N BBlT - + N BB7T+ l + 

where S = —r]cpsin2(j)i (rjcp is the CP-eigenvalue of the B° mode), x = Am/T, 
and A, a measure of direct CP violation, is zero in the SM. 

Neutral B 7 s are reconstructed in the following modes and submodes: B° — > 
J/ipKs; J/ip — > e + e~ , ■ Figure 60 shows the distributions in M m ; ss for (center) 

B°ir + and (right) B°n~ combinations, respectively, where the fits yield a total of 
21.5 ±6.8 events. The asymmetry is found to be A BBn = 0.28 ±0.28, corresponding 
to sin2(/>i = 0.57 ± 0.58 ± 0.06. This result establishes a new time-independent 
method of measuring sin20i. The value is consistent with measurements in T(4S) 
data. 



§9. New resonances 



9.1. Charmonium physics 

In e + e~ collisions at CM energies near ^/s ~ 10.58 GeV, there are a number 
of ways to produce final states that contain a cc quark pair. These include: i) 
P-meson decays, in which b — > ccs is a favored transition; ii) 77 fusion, which is 
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Fig. 60. Distributions in M m i ss of tagged B°n candidates for (left) simulated BBir, BB'ty, and 
B*B*ir events and 121.4 fbr 1 of data, (center) B°ir + and (right) B°ir~ . 



proportional to the square of the quark charge and, thus, favors production of cc 
and uu pairs over ss and dd pairs; hi) near-threshold s-channel cc production via 
initial-state radiation; and iv) cc associated production with J/tp mesons in e + e~ 
annihilation, which Belle found to be the dominant mechanism for J/tp productions 
in e + e~ annihilation near y/s = 10 GeV. Belle exploited all four of these processes to 
make a series of interesting discoveries related to the spectroscopy and interactions 
of cc charmonium mesons. 

9.1.1. First observation of the T) C (2S) 

Prior to 2002, the only "positive" observation of the rj c (2S), the first radial ex- 
citation of the charmonium ground state meson, the rj c , was a peak in the 7 energy 
spectrum from exclusive ip(2S) — > 7X decays reported by the Crystal Ball Experi- 
ment. 258 ^ However, this result was somewhat suspicious since the hyperfine tp(2S)- 
rj c (2S) mass splitting inferred from the measured mass value, AMhf s (2S) = 92 ± 
5 MeV, is substantially higher than the theoretical expectation of AM^ ory (2S) ~ 
58 ± 8 MeV; see, e.g., Ref. 259 ) In 2002, Belle reported the observation of a higher- 
mass r) c (2S) candidate in the rj c (2S) — > KsK^n 11 mass distribution produced via 
the B Kr) c (2S), r? c (25) -> Ks^n^ decay chain (see Fig. 61(left)). 260 ) Belle 
subsequently observed a signal at the same mass in the J/ip recoil mass spectrum 
for inclusive e + e~ — > J/ipX processes, 261 ^ shown in the right-hand panel of Fig. 61. 

The original Belle rj c (2S) signal has since been confirmed by a number of reports, 
including a higher statistics Belle study of B — > KKsK^tt^ decays. 262 ) The current 
PDG world- average hyperfine splitting value, AM^ G (2S) = 49 ±4 MeV, 15 ) is close 
to theoretical expectations and inconsistent with the Crystal Ball result, which is 
now generally thought to have been incorrect. 

9.1.2. TheX(3872) 

The X(3872) was first observed by Belle 263 ) as a small narrow peak in the 
it + it~ J/tp invariant mass spectrum from B — > Ktt + tt~ J/tp decays shown in the left- 
most panel of Fig. 62. It was subsequently confirmed by CDF, DO, and BaBar. 264 ) 
Other X(3872) decay modes that have been identified include the radiative de- 
cay, X(3872) — > 7«//^, 265 )' 266 ) which establishes its charge conjugation parity as 
C = +1, subthreshold decays to uiJtp, 267 ^ and the decay to open charm, X(3872) — > 
D*°D° . 268 ). 269 ) The Belle signals for X(3872) -> jj/tp are shown in the right panel 
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Fig. 61. Left: The KsK+tt^ mass distribution from B KKsK*^ decays. 260 ' The large peak 
on the left is the rj c ; the smaller peaks on the right are the J/psi (around 3.1 GeV) and r/ c (2S) 
signals. Right: The J/ip recoil mass spectrum in inclusive e + e~ — > J/ipX processes. 261 ' A fit 
with rjc, XcO, and r) c (2S) contributions is shown as a solid curve. The dashed curve in the figure 
corresponds to the case where the contributions of the J/ip, Xci, Xc2, and ip(2S) are set at their 
90% C.L. upper limit values. The dotted curve is the background function. 

of Fig. 62. Angular correlation studies by CDF 270 ) and Belle 271 ) indicate a preferred 
quantum number assignment of J PC = 1 ++ , although 2 h cannot be ruled out. 
The only available 1 ++ cc charmonium assignment for the X(3872) is the x'd- How- 
ever, the 3872 MeV mass value is significantly lower than the expected x'ci mass of 
3905 MeV, a value that is pegged to the measured 3929 ± 6 MeV mass of its J = 2 
multiplet partner, the x' C 2i which was discovered by Belle in 2006 (see below). A 
x' c i mass of 3872 MeV would imply that the mass splitting for the radially excited 
Xcj(2P) multiplet is larger than that for the % c j(lP) multiplet, contrary to expec- 
tations from potential models and lattice QCD. 272 ) There are similar problems for 
the J PC = 2~ + assignment, for which the only available cc level is the r/ C 2, the 1 L>2 
state. In this case, the 3872 MeV mass value is too high compared to the expected 
value of 3837 MeV, an expectation that is tightly constrained by the measured mass 
of its 3 L>i multiplet partner, the well established ^(3770). 

The lack of a natural charmonium assignment and the close proximity of the A(3872) 
mass, 3871.68 ± 0.17 MeV, 15 ) to the D*°D° mass threshold, 3871.94 ± 0.35 MeV, 15 ) 
has led to speculations that the V(3872) is a loosely bound D*°D° molecule-like 
structure; see, e.g., Ref., 273 ) although other interpretations have been proposed; see, 
e.g., Refs. 274 ) 

9.1.3. TheV(3940) 

The V(3940) was first observed by Belle as the near-threshold peak in the uJ/if) 
invariant mass distribution in B — > KujJ/tf) decays, 275 ) as shown in the left panel 
of Fig. 63. This observation was subsequently confirmed by BaBar. 276 ) The Belle 
experiment reported a similar peak in the near-threshold oj J/Y> mass distribution pro- 
duced in the two-photon process 77 — > ojJ/iJj 277 ^ (see Fig. 63 (right)). Although the 
mass of the V(3940) is well above the open-charm threshold, decays to DZ) 278 )> 279 ) 
and D*Z) 268 ) have not been seen; in the latter case, a 90% C.L. upper limit of 
B(V(3940) -»• D*D) < 1.4£(Y(3940) -»• uJ/ip) has been established. This limit 
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Fig. 62. Left: AM = M{ty+tt- £ + T) - M(£ + £~) distributions for B -> K-k^-k" J/ip, J/ip 

decays for a) data and b) inclusive B -> JipX MC. 263) The peak near AM ~ 0.6 GeV is due to 
tp' — > -k + -k~ J/ip decays; the peak at AM ~ 0.75 GeV in the data, which does not show up in 
the MC, is due to X(3872) -5- tt+tt" J/ip. Right: M(^J/ip) distributions for a) B+ -> K+-yJ/ip 
and b) B° K°jJ/i} decays. 265 ' 



and the rate of production in two-photon processes, implies that the partial width 
to u)J/ip is large, namely 1^(1^(3940) — > uJ/tp) > 1 MeV, which is very large for 
charmonium. 




M(d/f) (MeV) M(oxJ/f) (MeV) 



Fig. 63. Left: The points with error bars show the M(ujJ/ip) distribution for B — > KtoJ/ip decays. 
The curve in a) shows results of a fit to a phase-space-like threshold function. The curve in 
b) shows the results of a fit with a Breit-Wigner resonance function included. 275 ' Right: The 
u}J/tp invariant mass distributions for the two-photon fusion process 77 — > uiJ/t/j. 277 ^ The bold 
solid curve shows results of a fit including a resonance (thinner solid curve) and the dot-dashed 
curve shows a fit to a phase-space-only distribution; the histogram shows J/-i/> sideband data. 



Belle's 77 — > y(3940) — > uJ/ip observation was confirmed by BaBar, which also 
included results of an angular analysis that favors a J FC = ++ quantum number 
assignment. 280 ^ The only available ++ cc assignment is the x' c o> f° r which the mass 
value is somewhat high, but, perhaps, acceptable. The x' c o ~^ D*D decay mode is 
forbidden by parity, but x' c o ~^ DD is allowed and expected to be a strongly favored 
mode, 281 ) so the lack of any prominent signal for it is a mystery. 278 ) 
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9.1.4. The Z(3930) candidate for the x' C 2 charmonium state 

The left panel of Fig. 64 shows the DD invariant mass distribution for the 
process 77 — > DD measured by Belle, 282 ) where a strong peak near 3930 MeV is 
evident. The right panel shows the | cos#*| distribution for events in the ±20 MeV 
mass interval centered at 3930 MeV, where 0* is the CM angle between the D meson 
direction and the beamline. Small values of | cos#*| are favored, which is consistent 
with expectations for a J = 2 resonance (shown in the figure as a solid curve). 
The mass, angular distribution, and the strong decay to DD are all consistent with 
expectations for the x' C 2i *- e -> the radially excited 2 3 P2 charmonium state. 




4.1 4 . 2 4 . 3 ' ' ' ' ' * i * i 
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Fig. 64. Left: Invariant mass distributions for DD pairs produced via the 77 — > DD two-photon 
process. The curves show fits to the data with (solid) and without a resonance term (dashed). 282 ' 
Right: The yield of events with 3.91 < M(DD) < 3.95 GeV versus [cos#*[. The curves are 
expectations for J — 2 (solid) and J — (dashed); the histogram shows the M(DD) sideband 
yield. 282 ' 



9.1.5. TheX(3940) 

Belle discovered a third meson state with mass near 3940 MeV, the X(3940), 
produced in association with a J /ip in e + e~ annihilation. The left panel of Fig. 65 
shows the distribution of masses recoiling from the J/tp in inclusive e + e~ —> J/ip X 
reactions. 283 ^ With a partial reconstruction technique, Belle was able to isolate 
samples of exclusive e + e~ — > J/ipDD and J/ipD*D events. The DD and D*D 
invariant mass distributions for these samples are shown in the right panels of Fig. 65. 
There is no sign of the V(3940) in the DD events, but there is a distinct signal for 
V(3940) D*D. 

To address the question of whether or not the V(3940) is the same state as the 
F(3940), a search 283 ) was made for e + e — > J/^uJ/tp. No signal for V(3940) — > 
uJ/1/) was seen and a 90% C.L. lower limit £(X(3940) D*D) > 1.7B(X(3940) 
ujJ/ip) was established, which is inconsistent with the corresponding upper limit 
for the y(3940) discussed above. This implies that the y(3940), produced in B 
decays and decaying to ojJ /ip, and the X(3940), produced in association with a J/ip 
and decaying to DD*, are distinct states. The only cc assignment available for the 
X(3940) is the r/ c (35), for which decays to D*D are expected to be dominant and 
decays to DD are forbidden by parity. However, the 3 S*i triplet partner state of the 
^(35) is the well established -0(4040), with a mass of 4040 ± 4 MeV. 15 ) Assigning 
the V(3940) as the n c (3S) would mean AM hfs (3S) = 98±8 MeV, i.e., twice as large 
as AMhf s {2S) (see above) and in strong disagreement with theoretical expectations. 
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Fig. 65. Left: The distribution of masses recoiling from the J/tfi in inclusive e + e — s- J/ipX reac- 
tions. 283 ' The solid curve shows the result of a fit that includes r/c, XcO, r/c(2S), and X(3940) 
resonance terms as well as a smooth background function that has a step at the DD threshold 
(dotted curve). Right: The a) DD and b) D* D invariant mass distributions from exclusive 
e+e' ->■ J/ipD^D annihilation. 283 ' The curves are fits that include possible resonance terms 
and the histograms are backgrounds determined from the D-meson sidebands. The dashed 
curves show: a) the 90% C.L. upper limit on the signal; b) the background function. 



9.1.6. Anomalous J = 1 states seen in initial-state-radiation processes 

In 2005, BaBar reported the discovery of a striking ir + tt~ J / peak near 4260 MeV 
in the initial-state-radiation process e + e~ — > 7j sr 7r + 7r~ «//?/>. 284 ) This observation was 
subsequently confirmed by CLEO 285 ) and Belle. 286 ) The cross section for e + e~ — > 
it + it~ J / ij) from the Belle paper is shown in the left panel of Fig. 66, where a promi- 
nent signal for the Y (4260) with a peak cross section of ~ 70 pb is evident. Curiously, 
the total cross section for e+e" annihilation into open charmed mesons shows no 
sign of a peak at 4260 MeV; the total cross section for open charm at the Y(4260) 
peak is about 3 pb, 287 ) which, taken together with the measured natural width 
r tot [y(4260)] = 95 ± 14 MeV, implies a 90% C.L. lower limit on the partial width 
T(Y(4260) -> 7T+7T- J/ip) > 1.6 MeV. 288) This is much larger than values that are 
typical for 1 charmonium states (e.g., ^(-0(3770) — > TT + ir' J/ip) = 53 ± 8 keV). 

BaBar also reported a similar peak in the TT + ir~ip(2S) cross section at 4325 MeV. 289 ) 
With higher statistics, Belle confirmed this (now called the Y(4360)), and found a 
second, higher mass peak, the V(46 60) 290 ^ (see the right panel of Fig. 66). Here too, 
there are no evident accompanying structures in the open charm cross sections near 
these masses. Another peculiar feature is that, with the currently available statistics, 
there are no signs of the Y(4260) in the 7r + 7r _ ^'(2S') channel or of the Y(4360) or 
V(4660) in the ir + Tr~J/ip channel. 

9.1.7. The electrically charged Z~ charmonium-like meson candidates 

In 2008, Belle reported peaks in the ip'ir~ and XciK~ invariant mass distribu- 
tions inB^ iP'ir-K (Fig. 67 (left) 291 )' 292 ) and B -»• Xci^~K (Fig. 67 (right), 293 ) 
respectively. If these peaks are meson resonances, they would necessarily have a 
minimal quark content of ccdu and be unmistakably exotic. Although in both cases 
the peaks have greater than 5a statistical significance, the experimental situation 
remains uncertain since none of these peaks have yet been confirmed by other ex- 
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periments. Analyses by BaBar of the same channels neither confirm nor contradict 
the Belle claims. 294 ) 
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Fig. 67. Left: The M^, n _ projection of the Dalitz plot with the K* bands removed is shown as 
data points. 292 ' The histograms show the corresponding projections of the Dalitz-plot fits with 
(red solid) and without (blue dotted) a Z~ — s- ^V - resonance term. The dashed histogram is 
the background. Right: The data points show the M Xcl7r - projection of the Dalitz plot with 
the K* bands removed. The histograms show the corresponding projections of the fits with 
(red solid) and without (blue dotted) two Z~~ x.ci n ~ resonance terms, the dotted histograms 
represent the contribution of the two Xci 7 !" - resonances. 



9.1.8. Studies of open charmed hadron pair-production via initial-state-radiation 

The observation of the 1"(4260) motivated a Belle program of measurements 
of exclusive e + e~ cross sections for charmed hadron pairs near threshold. Belle 
presented the first measurements of exclusive cross sections for the production of 
charmed-hadron pairs in electron-positron annihilation in the vicinity of the thresh- 
old for open-charm production performed at CM energies near the T(4f>) resonance 
using the initial-state-radiation process. The continuous energy spectrum of this 
radiation allows investigating the production of charmonium with quantum numbers 
J PC = 1 over the whole energy range. The electromagnetic suppression of hard 
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photon radiation is compensated by an enormous integrated luminosity collected at 
the .B-factories, and selection criteria specific for the ISR processes provide high effi- 
ciency at considerable suppression of the background. Taken together, these factors 
resulted in measurements that are competitive in precision with the CLEOc 285 ) and 
BESII 295 -* experimental data in which charmed-hadron cross sections were measured 
using e + e~ energy scans with and without electromagnetic suppression. 

The exclusive e+e" cross sections to DD (D = D° or D + ), D + D*~, D*+D*-, 
D D~7t + , and D°D*-tt + final state using ISR, 296 ^ 2 ") shown in Fig. 68, have no 
evident peaks that can be associated with any of the above-mentioned Y states, 
contrary to expectations for conventional J PC = 1 charmonium states with such 
large masses and total widths. 




Fig. 68. Exclusive cross sections for charmed-hadron pair production measured in Belle. ' ' 
Left: DD, D = D° or D+ (upper); D + D*~ (middle); D* + D*~ (lower). Right: D°D~n + 
(upper); D°D*~tt + (middle); A£A~ (lower). The vertical dashed lines indicate the mass values 
of the established charmonium states: V(4040), ^(4160), and t/>(4415). 



In 2008, the Belle collaboration reported the observation of a significant near- 
threshold peak, called the X(4630), in the e + e~ — > A+A~ exclusive cross section 
shown in the lower right-hand panel of Fig. 68. 300 ) It remains unclear whether or not 
this observed peak is a resonance. In particular, peaks near the baryon-antibaryon 
pair mass threshold are observed in many processes, including three-body baryon 
decays of B-mesons. 128 ) The mass and width of the X(4630) peak determined under 
the assumption that the X(4630) is due to a resonance are M = (4634 ± 10) MeV/c 2 
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and r = (92 ± 40) MeV. These values agree within errors with the mass and the 
full width of the y(4660) peak seen in the y(4660) -> vr+vr-^' decay channel 290 ) 
as mentioned above. Such a coincidence (including quantum numbers) may not 
be accidental, although the possibility that the X(4630) and Y(4660) peaks have 
different origins cannot be ruled out. Among possible conventional interpretations, 
it has been suggested that the A (4630) is the ip{5S) or ip(6S) 1 charmonium 
state, 301 - 1 or a threshold effect caused by the presence of the ip(SD) state with mass 
slightly below the A~ threshold. 

9.1.9. Summary 

This section has highlighted only a fraction of the charmonium and charmonium- 
related results from Belle. In addition to the observations described above, Belle 
reported a number of other observations related to charmonium. A near-threshold 
peak was found in the D*D* e + e~ — > J/ipD*D* annihilation process. 302 ) A Belle 
search for the Y(4140) — a <pJ/ip resonance reported by CDF 303 ) — in the (f>J/ip mass 
distribution produced via the 77 — > <f)J/ip two-photon process found no evidence for 
the F(4140) but, instead, uncovered a 3.2a significant peak at higher mass that was 
dubbed the X (4350). 304 ) Belle cross section measurements of exclusive processes of 
the type e + e~ — > J / iprjc 261 ^ and e + e _ — > J/tpD^ l)(*) 305 ) found order-of-magnitude 
disagreements with NRQCD predictions 306 ) and have had a profound impact on 
subsequent developments in the theory; see, e.g., Ref. 307 ) A recent study of the 7Xci 
mass distribution in the l?-meson decay process B — > K'jXci found strong evidence 
for the long-sought-for tp C 2, the 3 Z?2 charmonium state. 308 ) 

In the original physics program planned for Belle outlined in the Belle Letter 
of Intent, 309 ) no mention was made of charmonium physics or searches for non- 
conventional, multi-quark meson states. Somewhat unexpectedly, thanks in part to 
the huge data samples provided by the KEKB collider, Belle turned out to be a 
powerful instrument for both conventional charmonium physics, and for uncovering 
a new class of charmonium-like states that have yet to be understood. 310 ) 

9.2. Bottomonium(-like) states 

As described in the previous section, most of the new charmonium states dis- 
covered in recent years at the ^-factories do not seem to have a simple cc structure. 
Although the masses of these states are above the corresponding thresholds for decay 
into a pair of open charm mesons, they decay readily into J/tp or tp(2S) and pions, 
which is unusual for cc states. In addition, their masses and decay modes are not 
in agreement with the predictions of potential models, which, in general, describe 
cc states very well. For these reasons, some of these charmonium-like states are 
probably more complex than simple quark-antiquark states and are candidates for 
exotic objects such as hybrid, molecular, or tetraquark states. Recently, Belle has 
made a series of exciting discoveries of new states in the bottomonium sector using 
its unique data sample taken around the T(5S) resonance. 

Bottomonium refers to bound states of bb quarks and is considered an excellent 
laboratory to study QCD at low energy. The spin-singlet states hb(nP) and rj^nS) 
alone provide information concerning the spin-spin (or hyperfine) interaction in bot- 
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Fig. 69. The inclusive M m j as spectrum with the combinatorial background subtracted (points with 
error bars) and the signal component of the fit function overlaid (smooth curve). The verti- 
cal lines indicate boundaries of the fit regions. The expected high-statistics reference signals, 
T(5S) — > Y(nS)Tr + iT~ where n — 1, 2, 3, as well as the newly observed hb(lP) and hb(2P) states 
are seen. 



tomonium. Measurements of the h^nP) masses would provide unique access to the 
P-wave hyperfine splitting, AMhf s {nP) = < M(n 3 Pj) > — M(n Px), the difference 
between the spin- weighted average mass of the P-wave triplet states (xbj(nP) or 
n 3 Pj) and that of corresponding hb(nP), or r^Pi. These splittings are predicted to 
be close to zero. Recently, CLEO observed the process e + e~ —> h c (lP)Tr + ir~ at a 
rate comparable to that for e + e~ —> J/^tt + tt^ in data taken at the ■0(4160) reso- 
nance. Such a large rate was unexpected because the production of h c (lP) requires 
a c-quark spin-flip, while production of J/ip does not. Belle previously reported 
anomalously high rates for e + e~ —> Y(nS)7r + ir~ (n = 1, 2, 3) at energies near the 
T(5S) mass. 311 ) If the T(nS) signals are attributed entirely to T(5S) decays, the 
measured partial decay widths P[T(55) — > T(n5)7r + 7r _ ] ~ 0.5 MeV are about two 
orders of magnitude larger than typical widths for di-pion transitions among the four 
lower T(nS) states. Using the large data sample collected at energies near the T(5S) 
resonance and motivated by the suggestive CLEO result, Belle decided to investigate 
the missing hb(mP) singlet bottomonium states. 312 ) 

We do not expect the hb(mP) states to have a large dominant exclusive decay 
mode, which would allow their reconstruction with high efficiency Instead, they 
are reconstructed inclusively using the missing mass (recoil mass) of the 7r + 7r~ pair. 
The 7r + 7r _ missing mass is defined as M^ iss = (Pr(ss) — P-k+tt-) 2 -, where Pr(5S) is 
the 4-momentum of the T(5S) determined from the beam momenta and P w + W - is 
the 4-momentum of the tt + tt~ system. The tt + tt~ transitions between Y(nS) states 
provide high-statistics reference signals as shown in Fig. 69. The hb(nP) states 
are also very clearly, and for the first time, observed here. The measured masses 
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of the h b (\P) and h b (2P), M = (9898.3 ± l.lll;?) MeV/c 2 and M = (10259.8 ± 
0.6lj;'o) MeV/c 2 respectively, correspond to hyperfine splittings that are consistent 
with zero. The processes T(5S) — > h b (mP)ir + k~ , which require a heavy-quark 
spin flip, are then found to have rates that are comparable to those for the heavy- 
quark spin conserving transitions T(5S) — > Y(nS)ir + ir~ . These observations differ 
from a priori theoretical expectations and strongly suggest that exotic mechanisms 
contribute to T(5S) decays. 

To understand the T(nS) and h b (mP) production mechanism at the T(5S) res- 
onance, it is necessary to study in detail the resonant structure of the T(5S) — > 
T(nS)ir + iT- and T(5S) -> h b {mP)TT+ tt' transitions. 313 ) In the case of T(5S) -> 
Y(nS)ir + ir~ , the T(nS) is reconstructed in the channel and one examines 

the 7r ± T(n5) mass spectra. This is illustrated for the T(2S) case in Fig. 70. Two 
charged bottomonium-like resonances, the Z&(10610) and ^(10650), are observed 
(Table XXII). A similar structure is found (Fig. 71) for the h b (mP)7r + 7r~ de- 
cay, where this time the appropriate observable is M m i ss (-7r T ), the missing mass 
of the opposite sign pion as the decays are reconstructed inclusively using the 
missing mass of the it + it~ pair. Production of the Z b s saturates the T(5S) — > 
h b (mP)7r + Tr~ transitions and accounts for the high inclusive h b (mS) production 
rate. All channels yield consistent results and weighted averages over all five chan- 
nels give M = 10607.2 ± 2.0 MeV/c 2 , r = 18.4 ± 2.4 MeV for the Z 6 (10610) and 
M = 10652.2 ± 1.5 MeV/c 2 , T = 11.5 ±2.2 MeV for the Z b (10650), where statistical 
and systematic errors are added in quadrature. The Z&(10610) production rate is 
similar to that of the ^(10650) for each of the five decay channels. Analyses of 
charged pion angular distributions favor the J p = 1 + spin-parity assignment for 
both the Z 6 (10610) and Z 6 (10650). 




10.4 10.45 10.5 10.55 10.6 10.65 10.7 10.75 
M(Y(2S)7i) max , (Gev/e?) 

Fig. 70. Comparison of fit results (open his- 
tograms) with experimental data (points 
with error bars) for events in the Y(2S) 
signal regions. M(T(25)7r) max is the max- 
imum invariant mass of the two Y(2S)ir 
combinations. The hatched histogram 
shows the background component. 
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Fig. 71. Comparison of fit results for events in 
the h b (lP) signal region. The Z 6 (10610) 
and Z b (10650) are clearly observed in both 
cases; the result of the fit is represented by 
the histogram. 
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Table XXII. Comparison of results on Z(,(10610) and .£[,(10650) parameters (mass and width in 
MeV, relative normalization and phase in degrees) obtained from T(5S) —¥ T(nS)-K + ir~ (n = 
1,2,3) and T(5S) — > hb(mP)-K + -K~ (m = 1,2) analyses. 



Final state 




r(2S)-K + -K- 




h b (lP)-K + n- 


h b {2P)n+n- 


M[Z 6 (10610)] 
r[Z 6 (10610)] 
M[Z 6 (10650)] 
r[Z 6 (10650)] 
Rel. norm. 
Rel. phase 


10611 ±4± 3 
22.3 ± 7.7±l;g 
10657 ± 6 ± 3 
16.3 ± 9-8±| g 
0.57 ± 0.2llg;S 
58 ± 43+4 


10609 ± 2 ± 3 
24.2 ± 3.1± 2 3 ° 
10651 ± 2 ± 3 
13.3±3.3l3;2 

0.86 ± O.lllo °o 
-13 ± 13±^ 7 


10608 ± 2 ± 3 
17.6 ±3.0 ±3.0 
10652 ± 1 ± 2 
8.4 ±2.0 ±2.0 
0.96±0.14lg^ 
-9 ± 191 


10605 ± 2 + i 

u 4 +4.5+2.1 

10654 ± 3la 
on q+5 4+2.1 

1.39 ± 0.37l° ?5 
107+44+3 

i8 '-57-12 


10599+.3+.4 
-, 0+10+9 

10 -8 -7 

106541^13 
19 ± 71" 

1 fi +0.6+0.4 
± - u -0.4-0.6 
101 +65 +74 
±ol -105-109 



These states defy a standard bottomonium assignment. In principle, a bottomo- 
nium particle's electric charge is zero; therefore, the minimal quark content of the 
Zf,(10610) and ^(10650) is a four-quark combination. Theoretical interpretations 
of these hidden-bottom meson resonances were proposed immediately after their ob- 
servation. The proximity (within a few MeV/c 2 ) of the measured masses of these 
unexpected new states to the open beauty thresholds, BB* (10604.6 MeV/c 2 ) and 
B*B* (10650.2 MeV/c 2 ), suggests a "molecular" nature of these new states, which 
can in turn explain most of their observed properties. In the case of a molecule, it 
would be natural to expect that Z^(10610) and (10650) states to decay respectively 
to BB* and B*B* final states at substantial rates. Recently, Belle reported prelim- 
inary results on the analysis of three-body T(5S) — > BB*ir {B + B*°ir~ , B~B*°ir + , 
B°B*-ir + and B°B* + tt-) and T(5S) -»• B*B*tt (B*+B*°tt- and B*-B*°ir + ) in- 
cluding an observation of the T(5S) -»• ^(10610)7^ -»• [BB*]^^ and T(5S) -»• 
Z^(10650)tt t — > [S*5*] ± 7r =F decays as intermediate channels. Evidence (with a sig- 
nificance of 4.9a ) for a neutral Z° (10610) decaying to T(2S)ir° has been also obtained 
by Belle in a Dalitz plot analysis of T(5S) — > T(2S)tx°tt° using their full T(5S) data 
sample. 314 ^ Its measured mass, M(Z£(10610)) = 10609lg ±6 MeV/c 2 , is consistent 
with the mass of the corresponding charged state, the (10610). 

The Zb states have also been interpreted as cusps at the B*B and B*B* thresh- 
olds and as tetraquark states. 

After observing that the decay T(5S) — > hb(nP)ir + ~tt~ proceeds via the Zb in- 
termediate resonances, Belle 315 - 1 exploited this information to look for the 7/&(l,2<S) 
resonances in the processes e + e~ — > hb{nP)ir + ir~ , hb(nP) — > r)b(mS)j. Here only the 
7T + , 7r~, and 7 are reconstructed and the requirement 10.59 GeV/c 2 < M m i ss (7r ± ) < 
10.67 GeV/c 2 helps to reduce the background significantly. The M m i ss (7r + 7r~) spec- 
tra are fitted for different M^ ss (it + it "7) bins to measure the hb(nP) yield. The 

hb(nP) yield peaks at M^ s (tt + tt~ 7) values corresponding to m Vb ^ mS ^ due to the 
hb(nP) — >■ r]b{mS)^ transitions (Fig. 72). 

The hb{\P) —7- i]b(lS)^ and hb{2P) — > rjb(lS)^ transitions are observed for the 
first time and first evidence for the rjb(2S) is obtained using the hb(2P) — > %(2S*)7 
transition. The mass and width parameters of the %(1«S*) and %(2S I ) are measured to 
be m Vb{ls) = (9402.4 ±1.5 ±1.8) MeV/c 2 , m %(25) = (9999.0 ± 3.5l 2 ;^) MeV/c 2 , and 
r Vb{ls) = (10.81^1^) MeV. Our value of the mass is about 11 MeV higher 
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Fig. 72. The h b (lP) yield versus M^ s (7r + 7r"7) (a), and h b (2P) yield versus M^ s (7r + 7r _ 7) in the 
r)b(lS) region (b) and in the r)b(2S) region (c). The solid (dashed) histogram presents the fit 
result (background component of the fit function). 



than the previous world average and the hyperfine splittings are 57.9 =L 2.3 MeV 
and 24.3^45 MeV for the 15 and 2S states, respectively, consistent with theoretical 
predictions. 

9.3. Others 

In addition to cc and bb states, Belle has also studied charmed mesons and 
baryons. They are copiously produced at KEKB either directly in e + e~ collisions or 
as products of B meson decays. At the 10.53 GeV CM energy, the cross section for 
prompt cc pair production exceeds that of bb, assuring large samples of ground and 
excited charmed states hadronizing from the produced cc quarks. Charm hadrons are 
usually studied inclusively, however such an approach often suffers from large back- 
ground. Charm production in B decays is governed by the Cabibbo-favoured b — > c 
transition. The restricted kinematics of T(4S) — > BB production enables selection of 
clean B samples. The fixed spin of the parent B constrains possible quantum num- 
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bers of daughter particles, simplifying spin-parity determinations. However, charmed 
states with high spin and highly excited charm states are suppressed in B decays. 

9.3.1. Charmed mesons 

The spectra of quark- antiquark systems are predicted using potential models, 
which attempt to model QCD features by describing the interquark potential. 316 ) 
Charmed mesons, having cu, cd or cs quark content, are heavy-light systems for 
which the models employ Heavy Quark Symmetry (HQS). In the limit of an infinitely 
heavy-quark mass, heavy-light mesons become similar to a hydrogen atom; which 
gives many theoretical simplifications. However, since the c quark mass is finite, 
HQS is only an approximate symmetry. An important consequence of its breaking 
is the -D( s ) — D* s ^ splitting. The orbitally excited P-wave multiplet (L = 1), denoted 

D**y is expected to consist of a broad J p = (0 + , 1 + ) doublet having total light-quark 

angular momentum j q = \ and a narrow (1 + , 2 + ) doublet with j q = \- 

Before the advent of the P-Factories, in addition to the ground state and 
D* s ^ mesons, only the narrow DJ\ doublets were established: (D-y (2420), D\ (2460)) 
and (D s i(2536), -D* 2 (2573)); the broad ones remained missing. Discovery of two 
narrow and unexpected states, the L>* (2317) + and _D s i(2460) + , began a renais- 
sance in charm spectroscopy. 317 ) They were found in the D^ir and D* + ir° final 
states, respectively, and were produced inclusively in the cc continuum. Spectra of 
the AM(D { * ) tt g ) = Mp^vr ) - M{D { * ] ) mass difference measured by Belle are 
shown in Fig. 73; prominent peaks at AM(D s ir°) ps 350 MeV/c 2 and AM(D*tt°) ps 
350 MeV/c 2 , are the L>* (2317) and D s i(2460), respectively. Their masses and upper 

limits on their widths, measured from fits to the Z\M(Z)i*V°), are summarized in 
Table XXIII. 

Observation of radiative D s ^y (Fig. 73) and dipion D s it + it~ decays of the D s i (2460) 
ruled out a J p = 1 assignment. For the D* (2317) no decay channel was found 
apart from the discovery mode. Such a decay pattern was consistent with spin-parity 
assignments of + and 1 + for the -D* (2317) and -D s i(2460) respectively, as expected 
for the P-wave cs doublet with j q = \. However, the measured masses were much 
lower than predicted by potential models and, thus, decays to D^K, expected to 
be dominant, were not permitted kinematically. Instead, the dominant decays into 




0.1 0.2 0.3 0.4 0.5 0.6 0.1 0.2 0.3 0.4 0.5 0.6 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

M(Ds it ) - M(Ds) (GeV/c 2 ) M(Ds* 7t°) - M(Ds*) (GeV/c 2 ) M(Dsy)-M(Ds) (GeV/c 2 ) 

Fig. 73. Distributions of AM(D s n°) (loft), AM(D* s n°) (middle), and AM(D s j) (right). His- 
tograms show data from the D s and/or it sideband regions. 
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Table XXIII. Parameters of D* (2317) and D s i(2460), compared with PDG parameters of j q = § 
states. 





D* s * 


Decay modes 


Mass (MeV/c 2 ) 


Width (MeV/c 2 ) 


o + (i) 


do (2317) 


D S 7V 


2317.2 ±0.5 ±0.9 


< 4.6 




L> s i(2460) 


D*7v,Dsj 


2456.5 ±1.3 ±1.3 


< 5.5 




D s i(2536) 


D*K 


2535.3 ± 0.2 


< 2.3 


2+(f) 


D: 2 (2573) 


DK 


2572.6 ± 0.9 


20 ±5 



isospin-violating modes resulted in very small widths. All this triggered exotic inter- 
pretations of these mesons as DK molecules, multiquark states, mixtures of P-wave 
cs meson with a csqq tetraquark, or chiral partners of Di*\ 318 ^ 

To clarify the nature of these states, Belle searched for them in exclusive B — > 
DD s j decays, where D s j denotes any excited charmed-strange meson. 319 ) These 
reactions proceed via b — > cW + — > ccs, which is the dominant cs production 
mechanism in B decays; here D** with j q = | are expected to be more readily 
produced than j q 
D s i(2460) in B - 



| states. Thus, one expected to observe the -D* (2317) and 



DD 



S J) 

states studied were D^ir 



if they were the missing cs doublet. The D s j final 
D { *\ and D { * ] 7T + 7r . Figure 74 shows distributions 
of D s j invariant mass for B candidates satisfying AE and Mj, c signal region re- 
quirements, and for the channels with significant signals found: D* (2317) — > D s ir°, 
D s i(2460) D*tt° and L> s i(2460) D s j. The L> s i(2460) helicity angle distribu- 
tion for the _D s i(2460) — > D s ^f mode (Fig. 74) showed that the data were consistent 
with the J = 1 hypothesis. Study of the D s j production rates in B — > DD s j decays 
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Fig. 74. Left: M(D sJ ) distributions for D 3 j final states: D s n° (top), D*tt° (middle), D a y (bot- 
tom). Hatched histograms show the AE sidebands. Right: the D s i(2460) — > D s ~f helicity 
distribution. Data points are compared with MC predictions for J = 1 (solid line) and J = 2 
(dashed) assignments. 



seems to support the interpretation of the D* (2317) and D s i(2460) as the orbitally 
excited cs j q = \ doublet. Although some of the models managed to reproduce the 
low masses of these states, 32 °) our understanding of cs spectroscopy still seems to be 
incomplete. 
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On the other hand, the corresponding j q = ^ doublet in the cu spectrum, dis- 
covered by Belle about the same time as the narrow D** states, has properties that 
perfectly match potential model predictions. The D** mesons, expected to decay 
dominantly into D^n final states, were studied at Belle in a full Dalitz plot analysis 
of B + — > Z)M~7r + 7r + decays. 321 ) To distinguish between the two identical final-state 
pions, D(*^~tt + combinations having minimal and maximal mass values were used 
as the Dalitz plot variables. The M 2 (D ( -*K) m i n vs. M 2 (D^7r)max plots for B can- 
didates within the AE-Mb c signal region, are shown in Fig. 75. Non-uniformly dis- 
tributed events indicate intermediate resonances emerging in the M 2 (Dir) m i n spec- 
trum. The fitted resonance contributions to the M(D^~7r + ) m i n projection are 
shown in Fig. 75. The Dir system was found to be composed of a tensor D^ and 
broad scalar state while the D*n system consists of a narrow axial D±, a tensor 
L>2, as well as a broad axial D l . The two broad states, observed for the first time, 
were consistent with the j q = \ P-wave cu doublet. The measured parameters of 
the D**° states are summarized in Table XXIV: differences between the D** and 
D** properties are striking. Belle also performed a similar analysis for the D** +, s 
produced in B° -»• Z)M°7r+7r-. 322 ) 




Fig. 75. Dalitz distributions for B~ — > D + tt~tt~ (first from the left) and B~ — > D* + tt~tt~ (third) 
signal region candidates. The corresponding M(D + ir~) min and M(D* + -K~) m in projections, 
with background subtracted, are shown as the second and fourth plots, respectively. Hatched 
histograms show the fitted resonance contributions. The open histogram is the coherent sum of 
all contributions. 



Table XXIV. Parameters of the D** mesons. 
J^(j q ) D" Decay modes Mass (MeV/c 2 ) Width (MeV/c 2 ) 

+ (|) Do (2400) Thv 2308 ± 17 ± 15 ± 28 276 ± 21 ± 18 ± 60 

l+(i) Dj(2420) D*n 2427 ± 26 ± 20 ± 15 384t \f ± 24 ± 70 

1+(|) Di(2420) D*n 2421.4 ± 1.5 ± 0.4 ± 0.8 23.7 ± 2.7 ± 0.2 ± 4.0 

2+(|) D2(2460) D M 7v 2461.6 ± 2.1 ± 0.5 ± 3.3 45.6 ± 4.4 ± 6.5 ± 1.6 



Studies performed by Belle allowed one to investigate important implications of 
HQS breaking. Theory predicts that the two 1 + mesons, with j q = ^ and j q = |, 
decay into D*tt in an S and D wave, respectively. Due to the finite c-quark mass, 
the observed (physical) 1 + states can be a mixture of such pure states and, thus, the 
resulting D 1 and D\ amplitudes are superpositions of S- and D-wave amplitudes. 
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The corresponding mixing angle was measured to be nonzero. 321 ) Similarly, mixing 
between the two cs axial states can be expected. An angular analysis performed for 
the .D s i(2536) + — > D* + K$ mode showed that, contrary to the HQS prediction of a 
pure D-wave decay, the 5-wave decay dominates. 323 ) 

Potential models also predict multiplets of higher orbital and radial excitations 
of charmed mesons. The first example of such a cs meson, the Z?* 1 (2700) + , was 
observed in the D°K + final state produced in doubly-charmed B + -4 D°D°K + 
decays. 278 ^ Its mass was measured to be 2708 ± 9l}g MeV/c 2 , while its width is 
108 ± 23+3i MeV/c 2 . The D^(2700) + 's spin-parity of 1~ was established from a 
study of its helicity angle. The M(D°K + ) spectrum together with the measured 
intermediate resonance contributions, as well as the 1)^(2700) helicity distribution 
are shown in Fig. 76. Observation of the 1^(2700) —> D*K decay with a rate 
comparable to that for DK, suggests that the -0^(2700) is a D* radial excitation. 324 ) 




2.4 2.6 2.8 3 3.2 3.4 -1 -0.5 0.5 1 
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Fig. 76. Left: Background-subtracted M(D° K + ) distribution for B + — ¥ D°D°K + with a contribu- 
tion from Ds 1 (2700) + (blue histogram), reflections from i/>(3770) (green) and t/'(4160) (yellow) 
decaying to D°D°, non-resonant contributions (brown and red). Right: D^i(27QQ) helicity 
distribution compared to predictions for J — (green), 1 (red) and, 2 (blue) spin assignments. 



9.3.2. Charmed baryons 

Charmed baryons provide a laboratory for the study of the dynamics of a light 
diquark in the environment of a heavy quark and allow one to test many theoretical 
predictions. 325 - 1 For the charmed baryons with cud, cdd or cuu quark content, the 
only states known before the start of the i?-Factories were the and I7 c (2455) 0,+ ' ++ 
ground states with J p = ^ , the | spin excitation i7 c (2520), as well as four A c 
excitations observed by CLEO in the A c irir final state. Two states, the 7l c (2595) 
and the yl c (2625), were identified as orbitally excited states, while the interpretation 
of yl c (2765) and 7l c (2880) remained unknown. Except for the A c , quantum num- 
bers of charmed baryons were not measured but, instead, either assigned based on 
model predictions or unknown. Since the predicted spectra are rich and dense, J p 
assignment for a given state is difficult and requires an experimental determination. 

The first such measurement, performed for the yl c (2880), is an excellent ex- 
ample of a comprehensive study of baryon properties. 326 ) Figure 77 shows the 
A~^tc + it~ invariant mass, with the A^ — > pK~ir + mode reconstructed. In addi- 
tion to the yl c (2880) signal, there are also peaks associated with the yl c (2765), 
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Fig. 77. Left: M(/l+7r + 7r~) distribution with /I+tt* within E c (2455)°' ++ signal (black) and side- 
band (red) regions. Middle: Helicity distribution of -A c (2880) — > X' c (2455)7r with fit results for 
the J = 5 (dotted), | (dashed) and § (solid) hypotheses. Right: /l c (2880) yield as a function 



as well as the /l c (2940) found by BaBar in the D°p final state. 327 ' The param- 
eters of the narrow baryons, obtained from a fit to the M(A c tt + tt~) distribution 
are: M Ac{2880) = 2881.2 ± 0.2 ± 0.4 MeV/c 2 , r Ac{2880) = 5.8 ± 0.7 ± 1.1 MeV/c 2 , 

M Ac(2940 ) = 2938.0 ± l-3±f;2 MeV/c 2 , r Ac(2940 ) = 13+* ± 27 MeV/c 2 . The mea- 
sured /l c (2880) helicity distribution (see Fig. 77) is consistent with the spin | hy- 
pothesis. The quark model predicts the lowest | and | + A c spin excitations at 
about 2900 MeV/c 2 , in agreement with the /L c (2880) mass. Distribution of the 
yl c (2880) yield as a function of the M{A^ir ± ), shown in Fig. 77, indicates contribu- 
tions from the I7 C (2455) and I7 C (2520). The measured yl c (2880) partial width ratio, 




n I i i I i i i I i i i I i i I i i i I i i i I i i I i i i I i i i I 
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M(A+ir) - M(A+), GeV/c 2 

Fig. 78. M(A c tt) — M(A C ) distributions for the A£ signal window (points) and scaled sidebands 
(red histogram). Insets show background subtracted distributions for the X' c (2800). The peaks 
at 0.43 GeV/c 2 are cross-feeds from yl c (2880) -> E c (2A55)tv where the pion from the i7 c (2455) 
decay is missing. 

Belle also studied excited charmed baryons decaying to A c ir final states. Fig- 
ure 78 shows distributions of the AM(A c ir) = M(A c ir) — M(A C ) mass differences 
for the A~£ir~, Afn , and /l + vr + combinations. 329 - ) Peaks near 0.51 GeV/c 2 were 
attributed to new baryons forming an isotriplet denoted as Z , c (2800) 0,+ ' ++ . The 
measured I7 C (2800) mass splittings relative to the A c and I7 C (2800) widths are sum- 
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Table XXV. Parameters of charmed baryons discovered by Belle. The 170(2800) masses were mea- 
sured with respect to the A c mass of 2286.46 ± 0.14 MeV/c 2 . 



Name 


Quark content 


Decay mode 


Mass (MeV/c 2 ) 


Width (MeV/c z ) 


I7 o (2800) u 
£ c (2800)+ 
i; c (2800) ++ 


cdd 
cud 
cuu 


A+TY- 

A+tt° 
AU+ 


515.4™*- +M Ac 
505.4±»;|±iS;* + M Ac 
5U.5 + _litH +M Ac 


fil + 18+22 
ul -13-13 
R9 +37+52 
oz -23-38 
7 r + 18+12 
' J -13-ll 


£ c (2980) + 


csu 


A+K~it + 


2978.5 ± 2.1 ± 2.0 


43.5 ± 7.5 ± 7.0 


£■0(2980)° 


csd 


A+K° S 7V + 


2977.1 ±8.8 ±3.5 


43.5 (fixed) 


^(3077)+ 


csu 


A+K~n+ 


3076.7 ±0.9 ±0.5 


6.2 ± 1.2 ±0.8 


~ c (3077)° 


csd 


AtK%K+ 


3082.8 ± 1.8 ± 1.5 


5.2 ±3.1 ±1.8 



marized in Table XXV. These new states could be members of the U C 2 triplet with 
J p = I with total angular momentum of the light diquark equal to two, which 
are expected to have AM(A c n) 0.5 GeV/c 2 and a width of 15 MeV/c 2 . Mixing 
of the S C 2 with other states predicted to lie nearby, could be a reason for the wider 
observed state. 

For charmed-strange baryons formed from csd or csu quarks, in addition to the 
ground states Sc ^°' + and the | + spin excitation £ , c (2645) 0,+ , there were also two can- 
didates for P-wave excitations, the S" c (2790) and ^(2815), observed in the S' c tt and 
£' c (2645)7r final states, respectively. Belle studied S c states decaying into A c Ktt, 330 ^ 
which requires the c and s quarks in the initial states to be carried away by different 
daughter particles. Two peaks in the M{A^ K~tt + ) spectrum shown in Fig. 79, were 
attributed to the new excited baryons denoted S cx (2980) + and S cx (3077) + . Their 
neutral isospin partners were found in A+K®ir~. The S' cx (2980) and S cx (3077) 
parameters, obtained from fit to the M(A c Kn) distribution, are summarized in Ta- 
ble XXV. The spin parity assignments for these baryons remain unknown. 
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Fig. 79. Distributions of M(A+K~tt + ) (left) and M(AtK%~) (right) with the fit curves overlaid. 

§10. Two-photon physics 

An e + e~ collider is also a 77 collider. Through measurements of two-photon 
collision processes, we can study hadron spectroscopy. Two-photon physics at Belle 
includes searches for new resonances, tests of perturbative QCD, and measurements 
of photon-meson couplings and form factors. In this section, we report our investi- 
gations of scalar resonances and QCD tests in meson-pair production processes from 
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two-photon collisions in the energy range between 1 GeV and 3 GeV, respectively. 
New resonances produced in two-photon processes are discussed in Sect. 9. 

10.1. Hadron physics and QCD 

The Feynman diagram for the two-photon process, 77 — > X at an e + e~ collider is 
shown in Fig. 80, where the reaction is regarded as a collision of two photons, each of 
which is emitted from one of the initial e + e~ beams, i.e. e + e~ — > e + e~jj — > e + e~ X. 
The CM energy of the two-photon collision system covers a wide range, and hence 77 
reactions can be measured over a continuous and broad energy range. Usually, two- 
photon measurements are performed by exclusively reconstructing the final-state 
particle system X in order to determine the collision energy of the two photons 
(W = Mx) for each event. 

Meson resonance formation processes are explored in measurements in the low 
energy region (W < 3 GeV). Since two or more overlapping resonances are often 
produced we extract each component from a partial-wave analysis, which takes into 
account interference. It is known that some light-quark scalar mesons, such as the 
/o(980) and ao(980), cannot easily be explained in a qq constituent model. The 
two-photon decay width r~~ of these light quark mesons, which is measured by two- 
photon processes, is the most important parameter that provides information on the 
internal structure of such mesons. 

In the higher energy region (W ^ 3 GeV), we study the properties of char- 
monia and search for new hadronic states with even charge-conjugation C. Since 
the contributions from resonances are relatively small in this region, we can test 
QCD by measuring the differential cross section of meson-pair production processes, 
77 — > MM' , which is calculated theoretically in a model with quark-pair produc- 
tion 77 — > qq followed by quark hadronization. The hadronization part is described 
by several different models based on perturbative and non-perturbative QCD. The 
Belle data sample has been used to perform such QCD tests with by far the highest 
statistics to date. 

10.2. Principles of a two-photon process measurement at an e + e~~ collider 

In a two-photon process at an e + e~ collider, photons emitted from the beam 
particles are always virtual, and the four-momentum squared (q 2 , which is the same 
as the invariant mass squared of the photon) is always negative. The virtuality of 
the photon Q 2 (= —q 2 ) is well approximated as Q 2 = AE^E' sin 2 |, where Eb is the 
CM beam energy, and E' and 6 are the recoil energy and the scattering angle of the 




Fig. 80. A two-photon collision diagram for the process e + e — > e + e X. 
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beam particle, respectively, when 9 is not very close to zero. However, the emission 
angle of the photon has a strong peak near 9 ~ (Q 2 ~ 0). When the 9 angles of 
both photons are small and the recoil e~ and e + are not detected, the reaction is 
regarded to a good approximation as a collision of two real photons (we call this case 
a "zero-tag event"). 

In a zero-tag event, the transverse momentum component (pt) of the final-state 
system X tends to be balanced, that is, close to zero. Requiring pt balance and a 
much smaller detected energy compared with that for the e + e~ beams, we can easily 
separate the two-photon signal process from background e + e~ annihilation processes. 
However, if W is greater than about 0.5-y/s for the e + e~ beams, measurement of two- 
photon processes becomes difficult due to the large background from annihilation 
processes and/or the small statistics of the signal. In the S-factory energy range, 
measurements up to W ^$ 4.5 GeV are feasible. 

The measured cross section a(e + e~ — > e + e~X) can be translated into a two- 
photon collision cross section using the relation : 

/dL 
cj(77 -»• X;W)-^dW, 

where dL^/dW is the two-photon luminosity function calculated in QED as a prob- 
ability density distribution for the CM energy of the two-photon systems, which are 
emitted from the incident e + e~. The two-photon cross section depends very weakly 
(logarithmically) on the e + e~ beam energy. 

10.3. Single meson formation process 

If only one meson is produced in a collision of two real photons, the quantum 
numbers of the meson, C, and spin-parity (J ) are restricted to be (even) ±+ or 
(odd, J / 1) ++ - The production of J = 1 mesons is forbidden. Thus, two-photon 
production is complementary to e + e~ annihilation processes where only 1 mesons 
are produced directly. 

In these processes, the production cross section of a meson is proportional to its 
two-photon decay width .T 77 via the relation: 

a(W) = 8n(2J + 1) final state) 

where Mr and are the mass and total width of the meson, and B is the branching 
fraction. 

In the zero-tag mode, we measure only the final-state particles from the decay 
of a produced meson. This significantly reduces backgrounds compared to the case 
of 77 inclusive meson production. The two-photon decay width of neutral mesons 
is a direct and sensitive probe of their internal structure, as mentioned above. In 
addition, detailed analyses of final states are useful to study the branching fractions 
and decay structures. 

10.4. Production of light-quark mesons 

Meson production through two-photon processes had been studied in the past at 
PEP, PETRA, TRISTAN, and LEP (see, e.g., the compilation in Ref. 331 )). However, 
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the more than three orders of magnitude larger statistics available at a U-factory 
compared to past experiments have qualitatively improved the analyses, allowing 
detailed studies of resonances that were impossible in the past. 

Figure 81 shows an example of the large Belle two-photon data statistics; here 
we give the integrated cross section (| cos 6**| < 0.6) for 77 — > it + it~ as a function 
of W, where 9* is the angle of the produced particle relative to one of the incident 
photons in the CM system of the two photons and W is the total CM energy. 332 )' 333 - 1 
This analysis used an early Belle data sample with an integrated luminosity of only 
85 fb _1 , (~ 9% of the full data). The Belle data have negligibly small error bars and 
a structure due to the /o(980) is clearly visible near W ~ 1 GeV, as shown in the 
inset. 
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Fig. 81. Integrated cross section (| cos#*| < 0.6) of 77 — > tt + tt . The inset shows an enlarged view 
of the Belle data near the /o(980) peak. A fit with a resonance parameterization is superimposed. 



10.5. Measurements of pseudoscalar-meson-pair production at Belle 

Belle has performed a study of 77 —> P1P2, where P1P2 are 7r + 7r - , 332 )~ 334 ) 
K+K-, and i-T^o^Hse) ^0^0^37), 338) ^0 ; 339) and w _340) The angular coverage 

for charged-particle-pair production is restricted to the range | cos6**| < 0.6 due to 
the limitations of the charged track triggers. On the other hand, for 7r°7r°, A, 

and 

rjrj, we can extend the angular range to | cos#*| < 0.8 or even to | cos#*| < 1.0 (full 
angular coverage) owing to the wider coverage of the calorimeter trigger for multi- 
photon final-state events. It should be noted that a wider angular coverage plays an 
essential role in separating partial waves. 

A study of resonance production in two-photon collisions gives several resonance 
(R) parameters: its mass, its total width, and r^B(R — > P\P2)- The latter is 
difficult to obtain otherwise. 
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10.6. Differential cross sections and partial wave amplitudes 

Partial waves with even angular momenta contribute to the cross section of 
77 — > P\P 2 . Up to G waves may be considered at low energy (W 2 GeV). 341 ) The 
differential cross section can then be written as 

^(77 -> P 1 P 2 ) = \SY ° + D Y 2 ° + Go y 4 °| 2 + \D 2 Y 2 2 + G 2 Y 2 \ 2 , (10-1) 

where S is the J = partial wave, Do and Go (D 2 and G 2 ) are the helicity-zero 
(-two) components of partial waves for J = 2 and 4, respectively, and Yj l s are 
spherical harmonics; partial waves determine the energy (W) dependence, while 
spherical harmonics govern the angular dependence. Because spherical harmonics 
are not independent of each other, the partial waves cannot be determined by only 
fitting the differential cross section. 
If we write 

ihl -(77 P1P2) = S 2 \Y °\ 2 + D 2 \Y 2 Q \ 2 + D\ \Y 2 2 \ 2 + G 2 \Y»\ 2 + G 2 \Y 2 \ 2 , 



47rd| cos 8* 

^ (10-2) 

we can determine the "hat amplitudes" S 2 , D 2 , D 2 , Gg, and G 2 by fitting differential 
cross sections in each H^-bin, because the |Yj"| 2 s are independent of each other. 
Spectra of hat amplitudes can give useful information on partial waves even though 
they contain terms arising from the interference of partial waves (S, Do, D 2 , Go, 
and G 2 ). 337 ) 

In order to obtain information on possible resonances, we have to parameterize 
the partial waves and then fit differential cross sections according to Eq. (10-1). Such 
analyses allow measurement of the two-photon widths of some mesons, including the 
/o(980) and a (980). 

The existence of the low- lying scalar nonet (/o(500) (or a), if* (800) (or re), 
/o(980), and ao(980)) is a long-standing puzzle, yet these scalar mesons are thought 
to play the role of a "Higgs boson in QCD", by spontaneously breaking the chiral 
symmetry of the QCD vacuum. 342 ) 

The measured two-photon widths of /o(980) and ao(980) are small (although 
the results have large systematic errors) compared to those of the /2(1270) and 
02(1320), as listed in Table XXVI. This pattern of widths supports a picture in 
which the low-lying scalar mesons are made of color-triplet diquark pair. 342 -* 

A more satisfactory way to derive information on partial waves is to do partial 
wave analyses utilizing hadron data of the past and fully taking into account theo- 
retical constraints; 343 ^ we eagerly await such analyses using the high-statistics data 
from Belle. 

10.7. QCD in the higher energy region 

In the higher energy region (W > 3 GeV) where resonance contributions are 
small, QCD can be studied by measuring exclusive two-body hadron production. It 
is believed that QCD gives reliable predictions at sufficiently high energy but the 
applicable energy is not known. Belle can measure two-photon processes up to W of 
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Table XXVI. Two-photon width (xB) 



Meson 


T 77 (xS) (eV) 


Ref. 


/o(980) 
do (980) 


T 77 = 286 ± 
r„B(w°) = 1281™ 


337) 
339) 


/a(1270) 


r 77 = 3030 ± 350 


15) 


a 2 (1320) 


r 77 = iooo ± 60 


15) 



4.5 GeV. S.J. Brodsky and G.R. Farrar predicted 

^ = s 2 -^f(9*) (10-3) 

for the hadron-pair production in a two-photon process at sufficiently high energy, 
using the Mandelstam variables s(= W 2 ) and t. 344 ) n c is the total number of elemen- 
tary particles involved in the initial and final states, eight for baryon-pair production 
(.". a~W~ 10 ) and six for meson-pair production (.-. cr~VF~ 6 ). S.J. Brodsky and 
G.P. Lepage (BL) also calculated the differential cross section for meson-pair produc- 
tion. 345 ) Their calculation was based on perturbative QCD where the perturbatively 
calculable 77 — > qq part is convoluted with the quark distribution amplitude. They 
obtained 



d|cos0*| s I sin 4 6»* 



, 2(eie 2 )(ei-e 2 ) 2 ^ 

H ^r^* 9{ - u 

sin y* 



+2(eie 2 )V(n}> (10-4) 



where Fm(s) is the electromagnetic form factor for a meson M, ej is the charge 
of a constituent quark, and 3 is a function that depends on the quark distribution 
function. For charged meson-pair processes this calculation predicts da jd cos 8* ~ 
sin- 4 r, and da{K + n-) / da{K+ K~) = (f K /U) A - The first term in Eq. (10-4), which 
is dominant for charged meson pair processes, does not depend on g because the de- 
pendence on the quark distribution function is absorbed into Fm- This prediction 
was improved by taking into account the effect of the s quark and modifying distri- 
bution functions. 346 ) Predictions for neutral meson-pair processes are not straight- 
forward, since the terms that include g are dominant. 

On the other hand, a non- perturbative calculation in the handbag model 347 ) 
(DKV) factorizes the non-perturbative hadronization part and gives the differential 
cross section 

da 8na 2 1 , „ . , l9 ,„„ „. 

\R M m(s)\ 2 . (10-5) 



d\ cos 6*\ ~ s sin 4 0*' MMV 

Although this model cannot predict absolute values for the cross sections, it gives a 
relation between annihilation form factors Rmm( s ) i n different processes. 

The Belle experiment has measured cross sections for 7r + 7r - , 334 ) tt°tt , 338 ^ 7771-°, 339 ) 
W) 340) x+if~, 334 ) KOjfO, 336 ) andpp 348 ) production in two-photon production. Be- 
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fore the Belle experiment no data were available to test these models due to limited 
statistics and poor particle identification capabilities. 

The angular distribution measurements are summarized in Table XXVII. The 
W~ n dependence of the cross section and ratios of cross sections are listed in Ta- 
ble XXVIII. The measured angular dependences agree with sin™ 4 0* except for the 
nn process. We obtained larger n values than the BL prediction of six, and in the 
neutral meson-pair process the value is close to the BC prediction of ten, which may 
be due to a significant higher order contribution in this energy region. 346 ) 

The ratios of cross sections asymptotically approach a constant as energy in- 
creases, but no model can systematically reproduce all the measured values. 

For baryon-pair processes, the measured n value is larger than the perturbative 
QCD prediction of ten, but decreases as W incr eases. 348 ) The angular distribution 
above 2.5 GeV agrees qualitatively with the perturbative QCD prediction but has a 
steeper rise. 

Table XXVII. Comparison between measured angular distribution and perturbative QCD predic- 
tion of sin -4 8* 



Mode 


sin- 4 (9* 


W (GeV) 


| cos 8* | 


Ref. 


7T + 7T~ 


Well matched 


3.0-4.1 


< 0.6 


334) 


K+K~ 


Well matched 


3.0-4.1 


< 0.6 


334) 


K° S K° S 


Matched 


2.4-3.3 


< 0.6 


336) 


jtV 


Better agreement with sin™ 4 8* ± bcosd* 
Approaches sin™ 4 8* above 3.1 GeV 


2.4 - 4.1 


< 0.8 


337) 


r,n° 


Good agreement above 2.7 GeV 


3.1-4.1 


< 0.8 


339) 


VV 


Poor agreement 
Close to sin™ 6 8* above 3.0 GeV 


2.4 - 3.3 


< 0.9 


340) 



Table XXVIII. Energy dependence of the measured cross section, (n value in ao oc VK™ n ) and ratios 
of ao between different processes, ao is the cross section integrated over the sensitive angular 
region. An SU(3) octet (mixture of octet and singlet with mixing angle —18°) is assumed for 
the r\ meson. 7?/ is the ratio of decay constants squared, f%/f%o- 



Process 


n or ao ratio 


W(GeV) 


| cos 6>*| 


BL 345) 


BC 346) 


DKV 347J 


7T + Tr™ 


7.9 ±0.4 ±1.5 


3.0-4.1 


< 0.6 


6 


6 




K+R- 


7.3 ± 0.3 ± 1.5 


3.0-4.1 


< 0.6 


6 


6 




K%K% 


10.5 ±0.6 ±0.5 


2.4-4.0 


< 0.6 


6 


10 




jtV 


8.0 ±0.5 ±0.4 


3.1-4.1 


< 0.8 


6 


10 




V n° 


10.5 ± 1.2 ± 0.5 


3.1-4.1 


< 0.8 


6 


10 




VV 


7.8 ±0.6 ±0.4 


2.4 - 3.3 


< 0.8 


6 


10 




PP 


12.41";* 


3.2-4.0 


< 0.6 


10 






K+K-/n+n- 


0.89 ±0.04 ±0.15 


3.0-4.1 


< 0.6 


2.3 


1.06 




K%K%/K+K- 
TrVVTr+TT™ 


~ 0.13 to ~ 0.01 


2.4-4.0 


< 0.6 




0.005 


0.08 


0.32 ±0.03 ±0.06 


3.1-4.1 


< 0.6 




0.04-0.07 


0.5 


rjri/-K TV 


0.48 ± 0.05 ± 0.04 


3.1-4.0 


< 0.8 


0.2472/(0.467?/) 
0.36 ~/ (0.627?^) 






0.37 ±0.02 ±0.03 


2.4-3.3 


< 0.8 







110 



J. Brodzicka et al. 



10.8. Summary and outlook 

Two-photon processes can be a background when studying CP violation, the 
main theme at a -B-factory, as well as for other physics topics. However, a de- 
tailed study of two-photon data can contribute much to the understanding of hadron 
physics in its own right as described above. The topics that can be addressed are 
divided into four categories: the search for and study of new or exotic particles, the 
production and decay structure of charmonia, the nature of light-quark resonances, 
and tests of perturbative QCD. The overwhelming statistics available at a i?-factory 
has opened a new era in two-photon physics. 

Our study so far has mostly been limited to collisions of two real photons; a vast 
unstudied region remains open for future investigation, in which one or both of the 
photons are virtual, i.e. the study of single and double tagged two-photon physics. 

§11. Summary 

The Belle experiment at KEKB is described in Sect. 2. Belle accomplished its 
main mission, which was the verification of Kobayashi and Maskawa's bold proposal 
that a single irreducible complex phase can explain all matter-antimatter asymme- 
tries (CP-violating phenomena). 

As discussed in detail in Sect. 3, Belle's observation of large time-dependent 
CP asymmetries in modes such as B — > J/ipKg (together with similar results from 
BaBar) in 2001 demonstrated that the KM proposal was correct and laid the foun- 
dation for their 2008 Nobel Prize in Physics. In addition, the results provided a 
theoretically clean measurement of one of the interior angles of the unitarity trian- 
gle, <p\ (or (3). After the accumulation of the one ab _1 data set, the measurements of 
CP asymmetries involving 0i became precision results and important calibrations 
for new physics studies. 

To check the consistency of the SM of particle physics, it is also necessary to 
measure the other two interior angles of the unitarity triangle, 02 ( or &) and 03 (or 
7). Although theoretical plans for the determination of these angles were proposed 
at the start of the P-factories, the final and most precise results were obtained by 
new methods that were not originally anticipated; e.g. for 03, the best sensitivity 
was obtained from Dalitz analysis of B — > DK, D — > Kstt + it~ decays. 

The development of the methods for determination of the length of the sides of 
the unitarity triangle also followed a somewhat unexpected path that was determined 
by the convergence of high statistics i?-factory data and theoretical insight. The 
results and methods used for \V c i,\ and \V u i,\ determination are described in Sect. 4. 

The results for the sides and interior angles of the unitarity triangle are consis- 
tent. However, reasonably large new physics contributions, of order 10% the size of 
the SM amplitude, are still allowed. In parallel with the work on fixing the weak 
interaction parameters of the unitarity triangle, Belle also completed a decade of 
studies and publications on rare decays, as described in Sect. 5. 

In rare decays for which the SM amplitude contribution is highly suppressed, 
the effects of NP could be clear and dramatic. Belle established the existence of a 
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number of highly suppressed processes including b —> d'y and b — > sl + l~ . In addition, 
as the data sample has increased, there have been a number of intriguing hints of NP 
in various channels, e.g. exclusive hadronic b — > s CP-violating modes, B — > tv, and 
B —7- K*£ + £~ , but so far there is no compelling evidence of NP at the current level 
of sensitivity in Belle. Exploration of NP will require the luminosity of SuperKEKB 
and Belle II. 

A i?-factory is also a high energy tau-charm factory and has the largest samples 
of r leptons and reconstructed charm. The results on r lepton physics are described 
in Sect. 6. Searches for lepton-flavor-violating (LFV) decays and CP violation in 
the r sector have reached an interesting sensitivity at Belle but so far no NP signals 
have been found. The foundation for Belle II explorations of this sector has been 
established. The results on charm are discussed in Sect. 7. The highlights include 
two classes of unexpected and unanticipated results: the discovery of D-D mixing 
and the existence of a large number of new charmonium-like resonances (Sect. 9). 
The latter was completely unexpected by the theoretical community and was guided 
by Belle data. 

Belle is also the world's leading two-photon facility. The results in this domain 
of physics are discussed in Sect. 10. Finally, KEKB's capabilities to operate in a 
range of center of mass energies allowed Belle to record a number of unique large 
data sets at the T(1S), T(2S), and T(55) resonances. The T(55) data were used, 
as expected, to study some properties and decays of B s mesons (Sect. 8). However, 
theorists did not anticipate that these data could be used to discover a series of 
peculiar bottomonium-like resonances or find the missing bottomonia states such 
as the 776 (25), hb(lP), and hb(2P). These discoveries in hadron spectroscopy are 
described in Sect. 9. 

In addition to establishing the KM model, measuring weak interaction parame- 
ters, and observing suppressed SM processes, the analysis of Belle data was marked 
by a series of unexpected discoveries driven by data. At the next stage in Belle II 
at SuperKEKB, the focus will shift to NP exploration. However, it is likely that the 
large increase in luminosity will also lead to unanticipated results and discoveries. 
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